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Introduction 



In any protocol aimed at manipulate or transmit information, symbols are encoded in states of some physical 
system such as a polarized photon or an atom. If these systems are allowed to evolve according to the laws of 
quantum mechanics, novel kinds of information processing become possible. These include quantum cryptog- 
raphy, teleportation, exponential speedup of certain computations and high-precision measurements. In a way, 
quantum mechanics allows for information processing that could not be performed classically. 

Indeed, in the last decade, we have witnessed a dramàtic development of quantum information theory, mostly 
motivated by the perspectives of quantum-enhanced communication, measurement and computation systems. Most 
of the concepts of quantum information were initially developed for discrete quantum variables, in particular quan- 
tum bits, which have become the symbol of the recently born quantum information technology. More recently, 
much attention has been devoted to investigating the use of continuous variable (CV) systems in quantum infor- 
mation processing. Continuous-spectrum quantum variables may be easier to manipulate than quantum bits in 
order to perform various quantum information processes. This is the case of Gaussian state of light, e.g. squeezed- 
coherent beams and twin-beams, by means of linear optical circuits and homodyne detection. Using CV one may 
carry out quantum teleportation and quantum error correction. The concepts of quantum cloning and entanglement 
purification have also been extended to CV, and secure quantum communication protocols have been proposed. 
Furthermore, tests of quantum nonlocality using CV quantum states and measurements have been extensively 
analyzed. 

The key ingredient of quantum information is entanglement, which has been recognized as the essential re- 
source for quantum computing, teleportation, and cryptographic protocols. Recently, CV entanglement has been 
proved as a valuable tool also for improving optical resolution, spectroscopy, interferometry, tomography, and 
discrimination of quantum operations. 

A particularly useful class of CV states are the Gaussian states. These states can be characterized theoretically 
in a convenient way, and they can also be generated and manipulated experimentally in a variety of physical 
systems, ranging from light fields to atòmic ensembles. In a quantum information setting, entangled Gaussian 
states form the basis of proposals for teleportation, cryptography and cloning. 

In implementations of quantum information protocols one needs to share or transfer entanglement among dis- 
tant partners, and therefore to transmit entangled states along physical channels. As a matter of fact, the propaga- 
tion of entangled states and the influence of the environment unavoidably lead to degradation of entanglement, due 
to decoherence induced by losses and noise and by the conseqüent decreasing of purity. For Gaussian states and 
operations, separability thresholds can be analytically derived, and their influence on the quality of the information 
processing analyzed in details. 

In these notes we discuss various aspects of the use of Gaussian states in CV quantum information processing. 
We analyze in some details separability, nonlocality, evolution in noisy channels and measurements, as well as 
applications like teleportation, telecloning and state engineering performed using Gaussian states and Gaussian 
measurements. Bipartite and tripartite systems are studied in more details and special emphasis is placed on the 
phase-space analysis of Gaussian states and operations. 

In Chapter^we introduce bàsic concepts and notation used throughout the volume. In particular, Cartesian 
decompositions of mode operators and phase-space variables are analyzed, as well as bàsic properties of displace- 
ment and squeezing operators. Characteristic and Wigner functions are introduced, and the role of symplectic 
transformations in the description of Gaussian operations in the phase-space is emphasized. 

In Chapter[2]Gaussian states are introduced and their general properties are investigated. Normal forms for the 
covariance matrices are derived. In Chapter[3]we address the separability problem for Gaussian states and discuss 
necessary and sufficient conditions. 

In Chapter|4]we address the evolution of a n-mode Gaussian state in a noisy channel where both dissipation 
and noise, thermal as well as phase-sensitive ("squeezed") noise, are present. At first, we focus our attention on the 
evolution of a single mode of radiation. Then, we extend our analysis to the evolution of a ri-mode state, which will 
be treated as the evolution in a global channel made of n non interacting different channels. Evolution of purity 
and nonclassicality for single-mode states, as well as separability threshold for multipartite states are evaluated. 



vi 



In Chapter[5]we describe a set of relevant measurements that can be performed on continuous variable (CV) 
systems. These include both single-mode, as direct detection or homodyne detection, and two-mode (entangled) 
measurements as multiport homodyne or heterodyne detection. The use of conditional measurements to generate 
non Gaussian CV states is also discussed. 

Chapter[ó]is devoted to the issue of nonlocality for CV systems. Nonlocality tests based on CV measurements 
are reviewed and two-mode and three-mode nonlocality of Gaussian and non Gaussian states is analyzed. 

In Chapter0we deal with the transfer and the distribution of quantum information, i.e. of the information 
contained in a quantum state. At first, we address teleportation, i.e. the entanglement-assisted transmission of an 
unknown quantum state from a sender to a receiver that are spatially separated. Then, we address telecloning, 
i.e. the distribution of (approximated) copies of a quantum state exploiting multipartite entanglement which is 
shared among all the involved parties. Finally, in Chapter[8]we analyze the use of conditional measurements on 
entangled state of radiation to engineer quantum states, i.e. to produce, manipulate, and transmit nonclassical light. 
In particular, we focus our attention on realistic measurement schemes, feasible with current technology. 
Throughout this volume we use natural units and assume h = c = 1. 

Comments and suggestions are welcome. They should be addressed to 

matteo . parisS f isica . unimi . it 

Corrections, additions and updates to the text and the bibliography, as well as exercises and solutions will be 
published at 

http: / /qinf . fisica. unimi . it/ "paris/ QLect . html 



Chapter 1 

Preliminary notions 



In this Chapter we introduce bàsic concepts and notation used throughout the volume. In particular, Cartesian 
decompositions of mode operators and phase-space variables are analyzed, as well as bàsic properties of displace- 
ment and squeezing operators 1 1 1 . Characteristic and Wigner functions are introduced, and the role of symplectic 
transformations in the description of Gaussian operations in the phase-space is emphasized J2). 



1.1 Systems made of n bosons 

Let us consider a system made of n bosons described by the mode operators a/t, k = 1, . . . , n, with commutation 
relations [afc,4] = Su- The Hilbert space of the system TL = <8>5J =1 Tk is the tensor product of the infinite 
dimensional Fock spaces Tk of the n modes, each spanned by the number basis { \m) k}m&h by the eigenstates 
of the number operator a k a,k- The free Hamiltonian of the system (non interacting modes) is given by H = 
2fe=i(4 afc + Position- and momentum-like operators for each mode are defined through the Cartesian 
decomposition of the mode operators = K\(qk + ipk) with k,\ S R, i.e. 

1k = 77- (Ofc + 4) , Pk = (ttfc ~ 4) . (1.1) 

zKi ZlKi 

The corresponding commutation relations are given by 



kk,Pl] = 7r^ S kl ■ (1-2) 



Canonical position and momentum operator are obtained for ki — 2 _1 / 2 , while the quantum optical convention 
corresponds to the choice /íi = 1. Introducing the vector of operators R = (qi,p-\ q n , Pn) T , Ea. (II. 2t rewrites 

as 

[R k ,Ri] = ^n k i, (i.3) 

where íï^i are the elements of the symplectic matrix 

" = -=(-°i D- ^ 

k=l ^ ' 

By a different grouping of the operators as S = (qi, . . . , q n ,Pi, ■ ■ ■ ,Pn) T , commutation relations rewrite as 

i 

[Sk, Si] = ——^ Jki , (1.5) 
where Jki are the elements of the 2n x 2n symplectic antisymmetric matrix 

'"(i - . 1 ")- 

1„ being the n x n identity matrix. Both notations are extensively used in the literature, and will be employed in 
the present volume. 
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Analogously, for a quantum state of n bosons the covariance matrix is defined in the following ways 

°kx = [cr]ki = ^({R k ,Ri}) - (Ri)(R k ) , (L7a) 
V U = [V] kl = ^({S ki Si}) - (Si)(S k ) , (1.7b) 

where {A, B} = AB + BA denotes the anticommutator and (O) = O = Ti[g O] is the expectation value of 
the operator O, with g being the density matrix of the system. Uncertainty relations among canonical operators 
impose a constraint on the covariance matrix, corresponding to the inequalities 1 3 1 

<x + -^í2>0, V--^J>0. (1.8) 

Ineqs. follow from the uncertainty relations for the mode operators, and express, in a compact form, the 
positivity of the density matrix g. The vacuum state of n bosons is characterized by the covariance matrix er = 
V = (4k| )~ 1 l 2n , while for a state at thermal equilibrium, i.e. v = (££C =1 v k with 



Vk 



-P a W 1 / Ni \ m 

E(rrV) i m ^ fc < m '' (L9) 



we have 



1 



2 Diag(27V 1 + l,2iV 1 + l,...,...,2iV„ + l,2iV„ + l) , (1.10a) 

4/í -i 



i 



V v = — Diag {2Nx + 1, . . . , 2N n + 1, 2N X + 1, . . . , 2iV n + 1) , (1.10b) 



where Diag(ai, 0,2, . . .) denotes a diagonal matrix with elements a^, k — 1,2,... and Nk — (e 13 — 1) _1 is the 
average number of thermal quanta at the equilibrium in the k-th mode . 

The two vectors of operators R and S are related each other by a simple 2n x 2n permutation matrix S = PR, 
whose elements are given by 

{Sk,2i-i k <n 
x i< ' (LH) 
<Wfc,2i í < n 

Sk,i being the Krònecker delta. Correspondingly, the two forms of the covariance matrix, as well as the symplectic 
forms for the two orderings, are connected as 

v = P(tP t , j= pvlp t . 

The average number of quanta in a system of n bosons is given by Y^k=i ( a k ak )- l n terms of the Cartesian operators 
and covariance matrices we have 

n 2n 2n 

J2(a{a k ) = ~ \ = «?£ fa + #) - \ ■ d- 12 ) 

fc=i 1=1 i=\ 

Eqs. (O can be generalized to define the so-called quadrature operators of the field 



x k <t> 



J-Lke-^ + ale^) , (1.13) 



i.e. a genèric linear combination of the mode operators weighted by phase factors. Commutation relations read as 
follow s 

1 

[xk<t>, afiíi] = 7T^2 *kl sin(-0 - 4>) ■ (1-14) 

Position- and momentum-like operators are obtained for tp = and <p = ir/2, respectively. Eigenstates | x) ^ of the 
field quadrature x<p form a complete set Vçí>, i.e. J R dx\x) ( p ( p(x\ = I, and their expression in the number basis is 
given by 

H.--'-'^'"^'-"». CU» 

Hk(x) being the fc-th Hermite polynomials. 
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1.2 Matrix notations for bipartite systems 

For pure states in a bipartite Hilbert space H i ® H2 we will use the notation 1 4 1 

|C))=^c H |fe)i®|/) 2) ({C\ = J2 c *kii(k\^2{l\, (1.16) 

H M 

where cjy = [C]ki are the elements of the matrix C and \k) r is the Standard basis of H r , r = 1,2. Notice that a 
given matrix A also individuates a linear operator from Hi to 7^2, givenby A = cife; | A;) 1 2 | - In the following 
we will consider Hi and H2 both describing a bosonic mode. Thus we will refer only to (infinite) square matrices 
and omit the indices for bras and kets. We have the following identities 

A <g> B\C)) = \ ACB T )) , ((C\A<g)B = ((ACB T \, (1.17a) 
((A\B))=Tt[A^B], (1.17b) 

where (• • • ) T denotes transposition with respect to the Standard basis. Notice the ordering for the "bra" in J 1 . 1 7 al > . 
Proof is straightforward by explicit calculations. Notice that A® B = [A® <g> B), and therefore is enough 
to prové dl . 17ai for A <g> 1 and 1 <8> B respectively. Normalization of state |C)) implies Tr[C^C] = 1. Also useful 
are the following relations about partial traces 

Tr 2 [|A))«B|]=ABt, Tri [|A))((£?|] = A T B* , (1.18) 

where (•••)* denotes complex conjugation, and about partial transposition 

(\C))((C\) e = (C® \)E(&®1) , 

where E = J2ki ® 10(^1 i s th e swa P operator. Notice that AB^ and A T B* in dl.l8> should be meant as 
operators acting on Hi and H2 respectively. Finally, we just remind that (A T ) T = (A*)* = (A')^ = A, and thus 
A t = ( A Ty = ( A *y^ a t = ^ A t)* = (A*)l, and A* = (A f ) r = (A T )t. 



1.3 Symplectic transformations 

Let us first consider a classical system of n partides described by coordinates (q±, . . . , q n ) and conjugated momenta 
(pi, . . . , p n ). If H is the Hamiltonian of the system, the equation of motion are given by 

Úk = T-, P* = ~, (fc = l,...,n) (1.19) 
where i denotes time derivative. The Hamilton equations can be summarized as 

Rk = Qkl -gjj- > Sk = —Jki-gg-, (1-20) 

where R and S are vectors of coordinates ordered as the vectors of canonical operators in Section fl~T1 whereas 
fi and J are the symplectic matrices defined in Eq. HAl and Eq. dl .6I >. respectively. The transformations of 
coordinates R' = FR, S' — Q S are described by matrices 

and lead to 

R'k = Fks^stFlt TT77" ! ^A; — —QksJstQlt "TTTT" ■ (1-22) 

Equations of motion thus remain invariant iff 

FÍ2F T =rj, QJQ T = J, (1.23) 

which characterize symplectic transformations and, in turn, describe the canonical transformations of coordinates. 
Notice that the identity matrix and the symplectic forms themselves satisfies Eq. (I1.23> . 

Let us now focus our attention on a quantum system of n bosons, described by the mode operators R or S. A 
mode transformation R 1 — FR or S' — QS leaves the kinematics invariant if it preserves canonical commutation 
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relations (II .31 or (II. 51 . In turn, this means that the 2n x 2n matrices F and Q should satisfy the symplectic 
condition ljl.23|) . Since fi T = fT 1 = fí from {L23fl one has that Dct [F] 2 = l 1 and therefore F 1 exists. 
Moreover, it is straightforward to show that if F, F\ and F2 are symplectic then also F , F T and F1F2 are 
symplectic, with F^ 1 = SÏF T Q ~ . Analogue formulas are vàlid for the «7-ordering. Therefore, the set of 2n x 2n 
real matrices satisfying (11.231 form a group, the so-called symplectic group Sp(2n, R) with dimension n(2n +1). 
Together with phase-space translation, it forms the affine (inhomogeneous) symplectic group ISp(2n,R). If we 
write a 2n x 2n symplectic matrix in the block form 



A B 
C D 



(1.24) 



with A, B, C, and D n x n matrices, then the symplectic conditions rewrites as the following (equivalent) 
condi tions 

AD T - BC T = 1 ( A T D C T B = 1 

AB T = BA T , \ A T C = C T A . (1.25) 

CD T = DC T { BD T = D T B 

The matrices AB T , CD T , A T C, and B T D are symmetric and the inverse of the matrix F writes as follows 

(D T — R T \ 
_ c t a t j ■ (1-26) 

For a genèric real matrix the polar decomposition is given by F = TO where T is symmetric and O orthogonal. 
If F G Sp(2n, R) then also T,0 £ Sp(2n, K). A matrix O which is symplectic and orthogonal writes as 

(X Y\ XX T + YY T = t 

u -y_ Y x) ' XY T -YX T =0 ' (L - Z/> 

which implies that U = X + ïY is a unitary n x n complex matrix. The converse is also true, i.e. any unitary 
n x n complex matrix generates a symplectic matrix in Sp(2n, R) when written in real notation as in Eq. (I1.27> . 

A useful decomposition of a genèric symplectic transformation F £ Sp(2n, R) is the so-called Euler decom- 
position 

F = °{ D D-) '' W 

where O and O' are orthogonal and symplectic matrices, while D is a positive diagonal matrix. About the real 
symplectic group in quantum mechanics see Refs. ]6]|5), for details on the single mode case see Ref. Q. 



1.4 Linear and bilinear interactions of modes 

Interaction Hamiltonians that are linear and bilinear in the field modes play a major role in quantum information 
processing with continuous variables. On one hand, they can be realized experimentally by parametric processes 
in quantum optical 1 8 9 1 and condensate [ 10][n][21[21[n- 1 systems. On the other hand, they generate the whole 
set of symplectic transformations. According to the linearity of mode evolution, quantum optical implementations 
of these transformations is often referred to as quantum information processing with linear òptics. It should be 
noticed, however, that their realization necessarily involves parametric interactions in nonlinear media. The most 
general Hamiltonian of this type can be written as follows 

h = è 9Í 1} 4 + E -9Í? } a Ui + E 9m a Wi + h - c - ■ ( L29 ) 

fc=l k>l=l k,l=l 

Transformations induced by Hamiltonians in Eq. dl.29> correspond to unitary representation of the affine symplec- 
tic group ISp(2n, R), i.e. the so-called metaplectic representation. Although the group theoretical structure is not 
particularly relevant for our purposes algebraic methods will be extensively used. 

Hamiltonians of the form ( II .291 contain three main building blocks, which represents the generators of the 
corresponding unitary evolutions. The first block, containing terms of the form H oc a) + h.c, is linear in the 
field modes. The corresponding unitary transformations are the set of displacement operators. Their properties will 
be analyzed in details in Section [T.4.1l The second block, which contains terms of the form g^a'b+h.c, describes 

1 Actually Det [F] = +1 and never —1. This result may be obtained by showing that if e is an eigenvalue of a symplectic matrix, than also 
e _1 is an eigenvalue 0. 
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linear mixing of the modes, as the coupling realized for two modes of radiation in a beam splitter. The dynamics 
of such a passive device (the total number of quanta is conserved) will be described in Section fï.4.21 This block 
also contains terms of the form g^a^a, which describes the free evolution of the modes. In most cases these terms 
can be eliminated by choosing a suitable interaction picture. Finally, the third kind of interaction is represented by 
Hamiltonians of the form <?( 3 W 2 + h.c. and g( 3 >(vtf + h.c. which describe single-mode and two-mode squeezing 
respectively. Their dynamics, which corresponds to that of degenerate and nondegenerate parametric amplifier 
in quantum òptics, will be analyzed in Sections ll.4.3l and ll.4.4l respectivelv. Finally, in Section fï.4.5l we briefly 
analyze the multimode dynamics induced by a relevant subset of Hamiltonians in Eq. dl.29> . corresponding to the 
unitary representation of the group SU(p, q). 

Mode transformations imposed by Hamiltonians dl. 291 can be generally written as 

R — » FR + d R , S^QS + d s , (1.30) 

where the d's are real vectors and F, Q symplectic transformations. Changing the orderings we have d s = P d R 
and Q = PFP, P being the permutation matrix ( 11.1 11 . Covariance matrices evolve accordingly 

<t > F (t F T , V^QaQ T . (1.31) 

Remarkably, the converse is also true, i.e. any symplectic transformation of the form dl.30t is generated by a 
unitary transformation induced by Hamiltonians of the form dl.29> |6|. In this context, the physical implication 
of the Euler decomposition dl.28t is that every symplectic transformation may be implemented by means of two 
passive devices and by single mode squeezers | Í3j. 

As we will also see in Chapter|2]the set of transformations coming from Hamiltonians ( ll.29l l individuates the 
class of unitary Gaussian operations, i.e. unitàries that transform Gaussian states into Gaussian states. 



1.4.1 Displacement operator 

The displacement operator for n bosons is defined as 



£>(A) = (g) J D fc (A fc ) (1.32) 

k=l 

where À is the column vector À = (Ai, . . . , A„) T , A& € C, k = 1, . . . , n and D k (Xk) — exp{Afca£ — A^afc}, are 
single-mode displacement operators; notice the definition of the row vector A^ = (A^ , . . . , A* ). 

Introducing Cartesian coordinates as X k = «^(a^ + ib k ) we have D(\) = D(A) = D(K) where 

D{A) = exp{2ÍK 1 K 3 fi T ríA} , D(K) = exp {-2í Ki k 3 S t JK} , (1.33) 

and 



A = (ai,bi, . . . ,a„,b„) T , K = (ai, . . . , a„, bi, . . . , b„) T , (1.34) 

are vectors in R 2n (ki has been introduced in Section [T"Tl . Canonical coordinates corresponds to k\ = K3 = 
while a common choice in quantum òptics is K\ = 1, K3 = 1 /2. The two parameters are not independent on each 
other and should satisfy the constraints 2kik 3 — 1 (see also Section [T~5l . In the following, in order to simplify 
notations and to encompass both cases, we will use complex notation wherever is possible. 
Displacement operator takes its name after the action on the mode operators 

D\X)a k D(X) = a k + \ k (ft = 1, .. . ,n) . (1.35) 

The corresponding Cartesian expressions are given by 

D\A)RD(A) =R + A, D\K)SD{K) = S + K . (1.36) 

The set of displacement operators D(\) with À S C™ is complete, in the sense that any operators O on H can be 
written as 

0= / ^Tr[OD(A)]£>t(A), (1.37) 

where 



X [0](A)=Tr[OD(A)] 



(1.38) 
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is the so-called characteristic function of the operator O, which will be analyzed in more details in Section [T"5l 
Eq. (II . 37t is often referred to as Glauber formula [16|. The corresponding Cartesian expressions are straightfor- 
ward 

r K 2n d 2n \ 

0= 3 X [0](A)£>t(A), (1.39a) 

í K 2n r\ 2n 

°= 3 - X [0](K)D\K), (1.39b) 

with d 2n A = d 2n K = nLi da k dbk and 

X [0](A)=Ti[OD(A)], x[0](K) =Tr[0 D{K)\ . (1.40) 

For the single-mode displacement operator the following properties are immediate consequences of the defini- 
tion. Let A, Ai, A 2 £ C, then 

£)t(A) = £>(-A), D*(X) = D(\*) , D T (X) = D(-X*), (1.41) 
Tr[£>(A)] =ttÓ (2) (\) , (1.42) 
£>(Ai)X>(Aa) = D(X 1 + A 2 ) exp {±(AiA| - X{ A 2 )} . (1.43) 

The 2-dimensional complex ó-function in Eq. dl.42> is defined as 



f d 2 X f d 2 X 

S^ 2 Uz) = / — 5- exp {X* z — z*X\ = / — 5- exp {i(X*z + z* A)} 

JC T Ve ^ 



(1.44) 



Setting A = a + zb and using Eq. ( 11.431 1 we have 



D*(X)D{z)D(X) =L>(z + 2a)cxp{-2ib(a + 3fie[z])}, (1.45a) 

D\X)D{z)D{X) = D{z) exp{zA* - z*X} , (1.45b) 

D(X)D(z)D(X)=D(z + 2X), (1.45c) 

D T (X)D{z)D(X) = D(z + 2ib) exp{2ia(b + 3m[z])} , . (1.45d) 



Matrix elements in the Fock (number) basis are given by 



(n + d\D(a)\n) = J e^H 2 a d ^(|a| 2 ) , (1.46a) 



(np(a)|n + d) = ^ e"^ 2 (-a*) d L^al 2 ) , (1.46b) 

(n\D(a)\n) = L n {\a\ 2 ) , (1.46c) 

Ln(x) being Laguerre polynomials. 

The displacement operator is strictly connected with coherent states. For a single mode coherent states are 
defined as the eigenstates of the mode operator, i.e. a\a) — a\a), where a £ C is a complex number. The 
expansion in terms of Fock states reads as follows 

00 k 

\a)=e-^ 2 Y / ^= \k). (1.47) 

k=0 »*' 

Using Eq. dl .35i it can be shown that coherent states may be defined also as \a) = D(a)\0), i.e. the unitary 
evolution of the vacuum through the displacement operator. Properties of coherent states, e.g. overcompleteness 
and nonorthogonality, thus follows from that of displacement operator. The expansion dl.47> in the number state 
basis is recovered from the definition \a) = D(a)\0) by the normal ordering of the displacement 

D{a) = e Qa V^ H V Q * a , (1.48) 

and by explicit calculations. Multimode coherent states are defined accordingly as \a) = D(a)\0) where \a) 
denotes the product state ®k\ a k)- Coherent states are (equal) minimum uncertainty states, i.e. fulfill dl .81 with 
equality sign and, in addition, with uncertainties that are equal forposition- and momentum-like operators. In other 
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words, the covariance matrix of a coherent states coincides with that of the vacuum state dkk = Vkk = (4/tf) 1 , 
Vfc = l,...,n. 

The following formula connects displacement operator with functions of the number operator, 

d 2 z f 11 + v 2 



tt(1 — i/) L 2 1-1/ 



exp --- \z\ z }D(z), (1.49) 



f d 2 z 

/ D{z)OD{z) 

Jc n 


= n Tr[n o] , 


[ — D(z)ODUz) 
'c n 


= Tr[0] I, 


í d 2 z 

/ — D(z)OD*(z) 
'c n 


= T , 


r d 2 z 

1 D{z)OD T {z) 

IP 7T 


= 0* . 



with special cases 

| )(0| = / — exp {-i|z| 2 } Z>(z) , (-1)-*» = / p. D(z) . (1.50) 
Jc n Jc 2lT 

Proof is straightforward upon using the normal ordering fll .48l > for the displacement and expanding the exponentials 
before integration. 

From Eq. dl.37l >, for any operator O, we have 

Tr [O f D{zj\ = Tr [O L> 1 {z)\ , (1.51a) 
Tr[0* D(z)] = Tt[OD*(z)} , (1.51b) 
Tr [0 T D{z)] = Tr [O D T (z)} . (1.51c) 

Using Eqs. Jl .5 lci . dl.45d> and dl.50> . other single mode relations can be proved 

(1.52a) 
(1.52b) 
(1.52c) 
(1.52d) 

where IT = (— ) ata i.e. H = J2 p (— ) p \p)(p\ is the parity operator. Using the notation set outin Section [T"2l we 
introduce the two-mode states \D(z))) = D(z) (g>I\l)) with = \p) ® \p)- Then we have the completeness 
relation 

/ — \D(z)))((D(z)\ =I®I. (1.53) 

JC 7T 

Other two-mode relations 

• fi* 

(1.54a) 

(1.54b) 

(1.54c) 
(1.54d) 

where |J)) = J2 P (~) P \p) ® Ip)» (' ' ' ) e denotes partial transposition, and E and F are the swap operator and the 
parity-swap operator respectively, the latter being defined as 

F = Y J (-Y +q \p)(<l\®\<l)(p\- d-55) 

p,q 

The action of E and fona genèric two-mode state is given by 

B(Wi®|02) =b)i® W)2 (1.56) 

p{ mi <8 1^) 2 ) - (-r lai iv>i <8> (-r^ 2 iv> 2 . a.57) 



A d Z r./ \ 


)D*(z) 


= \t)){{i\ = Y J \p){q\®\p){q\, 


/ D(z) (8 

/C T 


>D T {z) 


= |J))((J|=^(-)^ |p>(g|®|p><íl 




2>D{z) 




/ — £>(*)■ 

JC T 


= F^(|J))((J|) e , 


f d2z n , ^ 

/ D(z)Ç. 

Jc v 


)D\z) 


= fs = (|l»((l|) e , 
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Notice that the operator associated to bipartite state |J)) is the parity operator defined above. Finally, notice that 
using properties of Hermite polynomials, it is easy to show that 

~ oo 

/ dx Ix^Wt = e- 2in *\n)\n) = |F*» , (1.58) 
e.g. |1» = |F » = J R dx \x)o\x) G3 and |J» = |F f )) = J R dx \x)~\x) f . 



1.4.2 Two-mode mixing 

The simplest example of two-mode interaction is the linear mixing described by Hamiltonian terms of the form 
H oc a)b + b^a. For two modes of the radiation field it corresponds to a beam splitter, i.e. to the interaction taking 
place in a linear optical medium such as a dielectric plate. The evolution operator can be recast in the form 



U(C) =exp{Ca t ò-C*a& t } , 



(1-59) 



where the coupling £ = <fi e 10 E C is proportional to the interaction length (time) and to the linear susceptibility 
of the medium. Using the Schwinger two-mode boson representation of SU(2) àlgebra |18|, i.e. J+ = a)b, 
J_ = (J+y = ab^ and J3 = |[J+, J-] = \{a) a — b'b), it is possible to disentangle the evolution operator 
1191 1201 |2"T1 . thus achieving the normal ordering either in the mode a or in the mode b 

U(Ç)= exp{C J+-CJ-} 

= exp|ï|| tan | C| J+} exp {log(l + tan |C| 2 ) J 3 } exp {-£ tan |C| J-} 



= exp \é B tan <j> a^b} (cos 2 
= exp {— e~ %e tan^aò^} (cos 2 



cxp{— e ie tanc/>aò'} 



exp 



{é 6 tant/xz 1 ^} 



(1.60) 



Eq. Jl .60i are often written introducing the quantity t = cos 2 <f>, which is referred to as the transmissivity of the 
beam splitter. Mode evolution under a unitary action can be obtained using the Hausdorff recursion formula 



, 1 ü e~ aA = B + a{A,B} + ^ [A, [A, B}] + ^ [A, [A, [A, B}]] 

a' 



(1.61) 
(1.62) 



where {{A, B}} = [A, B] and {{A k , B}} = [A, {{A k ~ 1 ,B}}]. Eq. Q2l generalizes to e aA B n e 



„olA r>n~—aA 



B" and 



e aA e B e~ aA 



1 22 1. The Heisenberg evolution of modes a and b under the action of U(Ç) is thus given 



tft(c) 



U(0 = B i 



(1.63) 



where the unitary matrix Bç is given by 



B c = 



cos< 



e sini 

COBÓ 



(1.64) 



Correspondingly, we have W{() S U(Ç) = N ( S and W{Ç )RU(() = N' ( R, where N' c = P 23 N c P 23 S. 
The 4x4 orthogonal symplectic matrix, obtained from dl.64> as described in Eq. 11.271 . is given by 



3te[.B c ] -9im[B c ] 
3m[B ç ] ïïe[B c ] 



(1.65) 



and describes the symplectic transformation of two-mode mixing, whereas P 2 3 is the permutation matrix 

/ 1 \ 

10 

23 ~ 10 

\ 1 ) 



(1.66) 
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The two-mode covariance matrices evolve accordingly, i.e. as a — ► TV'ç cr N'ç and V — > Nç V Nç T . If p is the 
two-mode density matrix before the mixer and g' — U (Ç) g W (() that of the evolved state it is straightforward to 
show, using Eq. (II .631 . that 

(a,ia) e , = ((Ja) e cos 2 <f> + (b^b) g sin 2 4> + (§(^6 é d + tia e~ ie )) e sin(20) , (1.67) 

(tfb) e , = ( a ïa) e sm 2 <P+ (tfb) g cos 2 <P- (±( a ïbe l9 + tfae- l9 )) e sm(2(j)), (1.68) 

and therefore 

(aU + tfb)g, =(ala + tfb) e , (1.69) 

(a?a-tfb) e > = (a^be 10 + tfa e~ ie ) g sm(2^) . (1.70) 

Eq. dl .69t says that the total number of quanta in the two modes is a constant of motion: this is usually summarized 
by saying that a two-mode mixer is a passive device. It also implies that the vacuum is invariant under the action 
of Í7(C), '•<?. U(Ç) |0) = |0), where |0) = |0) ® |0). The two-mode displacement operator evolve as follows 

U(CY D(X) £7(0 = D(B\X) , (1.71) 

and thus the evolution of coherent states is given byC7(C)|a)) = \B\a)). Analogously, E/t(C) D(A) U{() = 
L»(A7^ 1 A) and D(K) U(Ç) = D(N^ 1 A), 

1.4.3 Single-mode squeezing 

We observe the phenomenon of squeezing when an observable or a set of observables shows a second moment 
which is below the corresponding vacuum level. Historically, squeezing has been firstly introduced for quadrature 
operators \ 23\, which led to consider the squeezing operator analyzed in this section. Squeezing transformations 
correspond to Hamiltonians of the form H oc (a^) 2 + h.c. The evolution operator is usually written as 

S(0 = exp{i£(at) 2 -ifa 2 } , (1.72) 

corresponding to mode evolution given by 

S , (í)tt5(0 = Ma + ^ , I 5 , K)« t ^) = ^ t + "*a I (1-73) 

where fi £ R, v G C, /i — coshr, v = e 1 ^ sinhr, £ = re*^. Using the two-boson representation of the SU(1, 1) 
àlgebra K + = i a t2 , K_ = {K + )\ K 3 = -^[K + ,K_] = \{a^a + it is possible to disentangle 
achieving normal orderings of mode operators 

S(0=exp{(K + -CK-} 

= cxp X+} cxp {log(l - tanh 2 \£\)K 3 } cxp {-jg K_} 

= ex P ( a t) 2 } ^o+h) cxp aJ , (1.74) 

from which one also obtain the action of the squeezing operator on the vacuum state |£) = 5(£)|0). The state |£) 
is the known as squeezed vacuum state. Expansion over the number basis contains only even components i.e. 




(1.75) 



Despite its name, the squeezed vacuum is not empty and the mean photon number is given by (Ç|aTa|£) = \v\ 2 , 
whereas the expectation value of quadrature operator vanishes (Ç\xg |£) = 0, V6*. Quadrature variance Ax 2 , is thus 
given by 

Ax 2 = {t\x 2 \t) = ^[e 2r cos 2 (e-\^) + e- 2r ún 2 {8-^/2)} . (1.76) 

Squeezed vacuum is thus a minimum uncertainty state for the pair of observables x^/ 2 and x^ /2+tt/2> f° r which 
wehave Aa;^, 2 = (4K 2 ) _1 e 2r and Aa; 2 / 2 +ix/2 ~ (4K 2 ) _1 e~ 2r , respectively. Applying the displacement operator 
to the squeezed vacuum one obtain the class of squeezed states \a, £) = D(a)S(Ç)\0). Squeezed states are still 
minimum uncertainty states for the pair of observables x^/ 2 and 2^/2+^/2- However, the photon distribution is 
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no longer characterized by the odd-number suppression of the squeezed vacuum. Notice that the evolution of the 
displacement operator is given by S^{^)D(X)S(0 = D(/j,\ - v\*), and that S{Ç)D{a) = D^a + va*)S(Ç). 
Therefore, application of the squeezing operator to coherent states leads to a squeezed state of the form 5(£)|a) = 
\fia + va*,Ç). 

Properties of quantum states obtained by squeezing number [24 1 and thermal state 1251 have been extensively 
studied. In general, if g' — S(£)gS(£) is the state after the squeezer, the total number of photon is given by 

(a + a)<?' = sinh 2 r + (2 sinh 2 r + l)((Ja) e + sinh(2r)(a 2 e"^ + a t2 e^) e . (1.77) 

Mode evolutions in Cartesian representation are given by R — > H^R and <x — > cr (S = R and <x = V 
since we have only one mode) where the symplectic squeezing matrix is given by 

/ Steíi/] Qm[v] \ 

s É = /ii a + fle fle= . (1.78) 

y ksm[i>\ - s Jte[^J j 

1.4.4 Two-mode squeezing 

Two-mode squeezing transformations correspond to Hamiltonians of the form H oc a^tf + h.c. The evolution 
operator is written as 

&(0 = expiat 6t-ra&} , (1.79) 
where the complex coupling £ is again written as £ = re 1 ^". The corresponding two mode evolution is given by 

4(0 ( £ ) S 2 (0 = S 2 J £ \ , (1.80) 

where S 2 ç denotes the matrix 



** = U- ; ) ■ ("« 

As for single mode squeezing we have /i = cosh r and 1/ = e 1 ^ sinh r. A different two-boson realization of the 
SU(1,1) àlgebra, namely K+ = aW, K_ = {K + )\ K 3 = -\[K + ,KJ\ = \{a)a + ò f ò + 1), allows to put 
S 2 (£ ) in the normal ordering for both the modes 

S 2 {Ç) =exp{^at6t| ^(a'a+btb+i) exp |_ii: a6 | (182) 

A two-mode squeezer is an active devices, i.e. it adds energy to the incoming state. According to Eqs. dl.80t and 
OH . with g' = S 2 (Ç)gSl(Ç) we have 

(a'a) g / = cosh 2 r{a)a) e + sinh 2 r(\ + (b'b) g ) 

+ ±sinh(2r)(aòe-^ +a t 6 t e^) e , (1.83) 
(6*6) s » = sinh 2 r(l + (a f a) e ) + cosh 2 r(ò t ò) e 

+ i sinh(2r)(a 6 e"^ + e^) e , (1.84) 

and therefore 

(a t a + b^b) g/ = 2 sinh 2 r (1 + (a t a + 6 t 6) e ) 

+ sinh(2r) (a b + aW e^) B , (1.85) 

(0+0-6+6)^» = (o t o-6t6) e . (1.86) 

The difference in the mean photon number is thus a constant of motion. The action of S 2 (£) on the vacuum can be 
evaluated starting from Eq. dl.82> . The resulting state is given by 

-. 00 , s k 

5 2 (£)|0) = — - \ k )®\ k ) d-87) 

fc^o 

and is known as two-mode squeezed vacuum or twin-beam state (TWB). The second denomination refers to the 
fact that TWB shows perfect correlation in the photon number, i.e is an eigenstate of the photon number difference 
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a)a—b^ b, which is a constant of motion. TWB will be also denoted as | A) ) where, adopting the notation introduced 
in Section fï~2l A is the infinite matrix A = i/l — |A| 2 À a a , with A = v/fi = tanh r. Often, by a proper choice 
of the reference phase, it will be enough to consider A as real. On the other hand, the first name is connected to a 
duality under the action of two-mode mixing. Consider a balanced mixer with evolution operator U(Ç = ? e 10 ), 
then we have 

UHj e 19 ) Sa(0 V(l e ie ) = S(Çe ie ) ® S(-Çe~ l9 ) , (1.88) 

where are single-mode squeezing operators acting on the evolved mode out of the mixer. In other words, a 
TWB entering a balanced beam-splitter is transformed into a factorized states composed of two squeezed vacuum 
with opposite squeezing phases 11261 . Viceversa, a TWB may be generated using single-mode squeezers and a 
linear mixer II27I . Using Eq. ( 11.581 1 we may also write 

|A)) = v^^2|A|° ta |F. )) = V^W^f" |F V ,)) • 

Finally, the symplectic transformation associated to the two-mode squeezer is represented by the block matrix X^. 
We have S^-RS^ = S 2Ç fí and S^SS^ = P23^2^P23S With 

where is defined as in dl .781 . and the inverse is evaluated using Eq. ( 11.261 1. 



1.4.5 Multimode interactions: SU(p, q) Hamiltonians 

Let us consider the set of Hamiltonians expressed by 



p 1 



H Pq = E 7« a * a l + E + EE a ^ + h - c - > ( L9 °) 

í<fc=l Kk=l fc=l 1=1 

where we have partitioned the modes in two groups a^, k = 1, . . . ,p and bi, l = 1, . . . , q, where p + q = n, with 
the properties that the interactions among modes of the two groups takes places only through terms of the form 
akbi + h.c. Hamiltonians dl .90t form a subset of Hamiltonians of the form ( 11.291 1. The conserved quantity is the 
difference D between the total mean photon number of the a modes and the b modes, in formula 

P Q 

fc=i 1=1 

The transformations induced by Hamiltonians ( ll.29l > correspond to the unitary representation of the SU(p, q) àlge- 
bra 1281 . Therefore, the set of states obtained from the vacuum coincides with the set of SU(p, q) coherent states 
i.e. 

\C pq ) = exp {-iH pq t} |0) = exp I E E ^ a k b \ ~ h - a } l°> . ( L92 > 
where (3ki are complex numbers parametrizing the state. Upon defining 



kl- 



tanh (ELil^l 2 ) 



E*=iI/3h| 2 
\Cpq) in Eq. dl .92t can be explicitly written as 



{m} {t} k=l 1=1 



r— 1 r— 1 r— 1 / 

where tk q = rrik — YlhL·i Pkh, { m } = { m i, m 2, m-p} and the sums over m and t are extended over natural 
numbers. In the special case q = 1, Eq. ( I1.93l > reduces to a simpler form, we have that \C p i) = \C P ) is given by 

10 aT ' a? '-^: vi r; m ] + - + "" )! iw; í» ,1.94, 

{ ™ } ^m 1 \m 2 \...m p \ ^ 

where N p = 1 — \ a k\ 2 is a normalization factor. 
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1.5 Characteristic function and Wigner function 

The characteristic function of a genèric operator O has been introduced in Eq. (II .381 . For a quantum state g we 
have x[í?](A) = Tr[gZ?(À)]. In the following, for the sake of simplicity, we will sometime omit the explicit 
dependence on g. The characteristic function x(A) is a ls° known as the moment-generating function of the signal 
g, since its derivatives in the origin of the complex plane generates symmetrically ordered moments of mode 
operators. In formula 



QP+Q 
^ d\ld\* qX ( Xj 



= Tr 


9 






A=0 






S- 



(1.95) 



For the first non trivial moments we have [a^a]s = ^{a)a + aa^), [aet^ 2 ]s = -|(a^ 2 a + aa) 2 + a^aa*), [a^a 2 ]s — 
i(a 2 a^ + aïa 2 + a^a) 1291 . In order to evaluate the symmetrically ordered form of genèric moments, one should 
expand the exponential in the displacement operator 

oo oo k / , \ 

m = E^r w - = E h E d ) Afe ^ 

OO OO \ fe A 

k=0 2=0 

Using Eqs. 11.371 . dl .381 and Jl .42l > it can be shown (see Section fï.5.1l for details) that for any pair of genèric 
operators acting on the Hilbert space of n modes we have 

Tr[Ox 2 ] = / d 2 "A x [Oi](A) X [0 2 ](-A) , (1.97) 

7T" J C n 

which allows to evaluate a quantum trace as a phase-space integral in terms of the characteristic function. Other 
properties of the characteristic function follow from the definition, for example we have xP\ (0) = Tr[0] and 

L Wr x[om = Tr [ ° n] L ~ |x[0](A)| = Tr[ ° 2] ' (L98) 

where II = ®JJ =1 (— ) a k ak = (— )££=i a ï ak is the tensor product of the parity operator for each mode. 

The so-called Wigner function of the operator O, and in particular the Wigner function associated to the quan- 
tum state q, is defined as the Fourier transform of the characteristic function as follows 

W[0}{a) = J exp |A f a + aUj X [0](A) ■ (1.99) 

The Wigner function of density matrix g, namely W[g] (a), is a quasiprobability for the quantum state. Using the 
formula on the right of Eqs. dl .98i we have that x[o\ (A) is a square integrable function for any quantum state g. 
Therefore, the Wigner function is a well behaved function for any quantum state. In other words, although it may 
assume negative vàlues, it is bounded and regular and can be used to evaluate expectation vàlues of symmetrically 
ordered moments. Starting from Eq. dl ,95ft and using properties of the Fourier transform it is straightforward to 
prové that 

í d 2n aW[g](a)a k (a*) 1 = Ti[q [{a^) l a k ]] . (1.100) 

More generally (see Section [T.5.H we have that 

Tr [O x 2 ] = Tr™ / d 2n aW[0 1 ](ot)W[0 2 ](a) . (1.101) 

Notice that the identity operator for n modes has a Wigner function given by W[I](a) = ir~ n . Indeed we have 
Tr[0] = J c „ d 2n a W[0](a). The analogue of Eq. dl ,37b reads as follows 

= 2"/ d 2n aW[0](a) D(a)UD^{a) , (1.102) 

JC" 

from which follows a trace form for the Wigner function 

W[0](a) = (-} Tt[0 D(a)UD\a)] . (1.103) 
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Other forms of the Eqs. (11.1 02b and Jl . 103b can be obtained by means oftheidentity D(a)UD^(a) = D(2a)U = 
n£»t(2a). 

The Wigner function in Cartesian coordinates is also obtained from the corresponding characteristic function 
by Fourier transform. Let us define the vectors 

X = (xi,vi, ■ ■ ■ ,x n ,y n ) T , Y = (xi, . . . ,x n ,yi, ■ ■ ■ ,y n ) T , (1.104) 

where ctk — K 2 (xk + iyk)- Notice that the scaling coefficients k 2 and k 3 are not independent one each other, but 
should satisfy 2k 2 k 3 = 1. To show this, consider the n-mode extension of Eq. Jl .44l > 

P rj 2n \ 

6^ 2n) (a) = — ^ exp{i(A*a + a*A)} (1.105) 



~2n 



from which follows that 



The identity 



P j2n \ 

5 (2n) {X)= -—^(2 K2 K 3 ) 2n exp{2t K2K3 A T X} . (1.106) 

J R 2„ (271-)^™ 



*(*) = f p 



(1.107) 



implies then that 2K2K3 = 1, as we claimed above. The corresponding definition of the Fourier transform allows 
to obtain the Wigner function in Cartesian coordinates as 

í d 2n A 

W[0](X)= -— ^ exp{iA T X} X [0](A), (1.108a) 

í d 2n K 

W[0](Y)= -— ^ exp{iK T Y} X [0](K). (1.108b) 

JR 2 ™ { Z7r > 

Notice that in the literature different definitions equivalent to Eq. (11.1 05i of the n-mode complex <5-function 
are widely used, which correspond to a change of coordinates in Eq. (II .1081 . As an example, if one consider 1 16 1 

í d 2n \ 

5^ 2n \a.)= —5— exp{À*a-a*A} (1.109) 



it follows that 

w[0}(x)= í -^r exp{iA T nx} x [0](n T A), (í.no) 

the same observation being vàlid for W[0] (Y). 

Let us now analyze the evolution of the characteristic and the Wigner functions under the action of unitary 
operations coming from linear Hamiltonians of the form l ll.29l l. If g is the state of the modes befaré a device 
described by the unitary U, the characteristic and the Wigner function of the state after the device g' = Ug W can 
be computed using the Heisenberg evolution of the displacement operator. The action of the displacement operator 
itself corresponds to a simple translation in the phase space. Using Eq. dl .45di we have 

X {D(z)gD\z)}(\) = xlQ}W expjzU-A^} , (1.111a) 
W[D(z) gD^(z)](X) = W[g](a - z) . (1.111b) 

In the notation of Eq. ( 11.301 1. Q — I and d = z, thus we have no change in the covariance matrices. In general, for 
the interactions described by Hamiltonians of the form dl .291 and, excluding displacements, we have 

X [UgW](A) = X [g](P2 3 F- l P 23 A) , 

(1.112a) 

x[UgW](K) = X [g](F- 1 K), 
W[UgW](X) = W[g](P 23 F- 1 P 23 X) , 

(1.112b) 

W[Ug W] (Y) = W[g] (F~ Y) , 

where F is the symplectic transformation associated the unitary U. Eqs. ( 11.112b say that the characteristic and the 
Wigner functions transform as a scalars under the action of U. For two-mode mixing, single-mode squeezing and 
two-mode squeezing the symplectic matrices are given by in Eqs. (II .651 . (II .781 and ( 11.891 respectively. 
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In summary, the introduction of the Wigner function allows to describe quantum dynamics of physical systems 
in terms of phase-space quasi-distribution, without referring to the wave-function or the density matrix of the 
system. Quantum dynamics may be viewed as the evolution of a phase-space distribution, the main difference being 
the fact that the Wigner function is only a quasi-distribution, i.e. it is bounded and normalized but it may assume 
negative vàlues. Unitary evolutions induced by bilinear Hamiltonians correspond to symplectic transformations 
of mode operators and, in turn, of the phase-space coordinates. Evolution of the characteristic and the Wigner 
functions then corresponds to transformation ( 11.1 12i . whereas non-unitary evolution induced by interaction with 
the environment will be analyzed in details in Chapter|4] 

1.5.1 Trace rule in the phase space 

The introduction of the characteristic and the Wigner functions allows to evaluate operators' traces as integrals in 
the phase space. This is useful in order to evaluate correlation functions and the statistics of a measurement since 
we are mostly dealing with Gaussian states and, as we will see in Chapter [5] also many detectors are described 
by Gaussian operators. In this Section we explicitly derive Eqs. (II . lOli and dl .97i . for the trace of two genèric 
operators in terms of their characteristics or Wigner function. The starting points are the Glauber expansions 
of an operator in terms of the characteristic or the Wigner functions, i.e. formulas (II . 371 and ( 11.1021 . For the 
characteristic function we have 

j2n \ f j2n 1 



f A \ r d X 

Tr[0 1 2 ]=/ -—L x [O x ]{\í) — ^ X[02](Ai) Tr[Z)(Ai) D(X 2 )] , 

JC" w J C" n 

^^ xtOiKAOx^KAr) 

xTr[^(A 1 + A 2 )]exp{A t 1 A 2 - A^Ai} , 

d 2n \ 

— - X [Oi](A) x [0 2 ](-A), (1.113) 



C" 



where we have used the trace rule for the displacement Tr[D(7)] = 7r"i5' 2n · ) (7). For the Wigner function we have 
Tr[Oi 2 ] - 2 2 " / d 2n oc 1 W[Oi]{eti) í d 2n a 2 W[0 2 ](a 2 ) 

JC" JC" 

x Tr[D(2ai)nn£)(-2a 2 )] , 
= 2 2ïl / d 2n a 1 d 2n a 2 W[0 1 }(a 1 )W[02}(a 2 ) 

JC 2 " 

x Tr[D(2ai — 2c*2)] exp |2a{c*2 — 2c4 a i j 7 
= n n í d 2n a.W[0 1 \{a.)W[0 2 ]{a) , (1.114) 

JC" 

where we have used the relations n 2 = I and 8 { · 2n \a^) — \a\~ 2n 5 (2n ï (j) with séí 



1.5.2 A remark about parameters k 

In order to encompass the different notations used in the literature to pass from complex to Cartesian notation, we 
have introduced the three parameters Kh, h = 1, 2, 3, in the decomposition of the mode operator, the phase-space 
coordinates and the reciprocal phase-space coordinates respectively. We report here again their meaning 

a-k = Ki(qk + ipk) , ak = K 2 (x k + iy k ) , X k = n 3 (a, k + ib k ) . (1.115) 

The three parameters are not independent on each other and should satisfy the relations 2/Í1K3 = 2k 2 ks = 1, i.e. 
K\ = «2 = (2k3)~ 1 . The so-called canonical representation corresponds to the choice k\ = k 2 = «3 = 2" 1 / 2 , 
while the quantum optical convention corresponds to m = n 2 = 1, K3 — 1/2. We have already seen that 
2/Í2K3 = 1; in order to prové that 2K1K3 = 1, it is enough to consider the vacuum state of a single mode 
and evaluate the second moment of the "position" operator (q 2 ) = Tr \g q 2 ~\ , which coincides with the variance 
(Aq 2 ), since the first moment (q) = vanishes. Starting from the commutation relation [q,p] = (2k 2 )^ 1 it is 
straightforward to show that the vacuum is a minimum uncertainty state with 



(Aq 2 ) = (4K 2 )" 1 



(1.116) 
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On the other hand, the Wigner and the characteristic functions of a single-mode vacuum state are given by 

W (x,y) = ^exp{-2 K 2 {x 2 +y 2 )} , Xo (a, b) = exp {-^(a 2 + b 2 )} . (1.117) 

Therefore, using the properties of W as quasiprobability, and of \ as moment generating function, respectively, 
we have 



(Aq 2 ) = f dxdy 4 x 2 W (x, y) = (4* 2 ,)- 1 , 
Jw 2 

d 2 

(Aq 2 ) = - — Xo ( a ,b) 



(1.118a) 

0- 

Xo(a,b) 

a=b=0 



4, (1.118b) 



from which the thesis follows, upon using Eqs. dl.ll6> . J 1 . 1 1 8 ai and dl . 1 1 8bt and assuming positivity of the 
parameters. Now, thanks to these results and denoting by er = Vq = (4/í 2 )~ 1 l2 n = (4/í|) _1 l2n = ^3 lz n the 
covariance matrix of the ?i-mode vacuum, we have that the characteristic and the Wigner functions can be rewritten 

as 

Xo (A)=exp{-ÍA T <T A} , Xo (K) = exp {-\K T V K} , (1.119) 

and 

, exp{-iX T cr n - 1 X} /2\" f 1 T , 1 
W (X)= 1 2 J = - expl--X T (Tn 1 X} , (1.120a) 

W- /0 „ t /t, =Tf7 T = - exp{-lY^ -V} , (1.120b) 



(27r)"K 2n ^/Det [V ] V 7 *" 

respectively, independently on the choice of parameters k^. This form of the characteristic and Wigner function 
individuates the so-called class of Gaussian states. The simplest example of Gaussian state is indeed the vacuum 
state. Thermal, coherent as well as squeezed states are other examples. The whole class of Gaussian states will be 
analyzed in detail in Chapter|2] 
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Chapter 2 

Gaussian states 



Gaussian states are at the heart of quantum information processing with continuous variables. The bàsic reason is 
that the vacuum state of quantum electrodynamics is itself a Gaussian state. This observation, in combination with 
the fact that the quantum evolutions achievable with current technology are described by Hamiltonian operators 
at most bilinear in the quantum fields, accounts for the fact that the states commonly produced in laboratories 
are Gaussian. In fact, as we have already pointed out, bilinear evolutions preserve the Gaussian character of the 
vacuum state. Furthermore, recali that the operation of tracing out a mode from a multipartite Gaussian state 
preserves the Gaussian character too, and the same observation, as we will see in the Chapter|4] is vàlid when the 
evolution of a state in a Standard noisy channel is considered. 



2.1 Definition and general properties 

A state g of a continuous variable system with n degrees of freedom is called Gaussian if its Wigner function, or 
equivalently its characteristic function, is Gaussian, i.e. in the notation introduced in Chapter^ 

r .. , exp{— ifa — a) T cr' 1 (at - a)} 

W[g](a) = ^ 2V a l (2.1) 

(27r)«yDrtKÏ 

with a. = K2 X, a — K2 X, where X is the vector of the quadratures' average vàlues. The matrix cr^ 1 is related 
to the covariance matrix er defined in Eqs. dl.7> by <j a = k 2 t. In Cartesian coordinates we have: 



W[g](X) = 1 2V ' , 1 ^ , (2.2) 

(2 7 r)« K l»V^R 

or equivalently: 

W[q] (Y) = 1 2V ; ; V ^- . (2.3) 

(27r)« 4 n v / Det[vT 

Correspondingly, the characteristic function is given by 1 

Xo(A) - exp {-ÍA T crA + iA T X} , X o(K) = exp {-\K T VK + iK T Y) . (2.4) 

In the following, since we are mostly interested in the entanglement properties of the state, the vector X (or 
Y) will be put to zero. Indeed, entanglement is not changed by local operations and the vectors X (or Y) can be 
changed arbitrarily by phase-space translations, which are in turn local operations. Gaussian states are then entirely 
characterized by the covariance matrix cr (or V). This is a relevant property of Gaussian states since it means that 
typical issues of continuous variables quantum information theory, which are generally difficult to handle in an 
infinite Hilbert space, can be faced up with the help of finite matrix theory. 

Pure Gaussian states are easily characterized. Indeed, recalling that for any operator Ok, which admits a well 
defined Wigner function Wk(a), we can write Tr[Oi O2] in terms of the overlap between Wigner function [see 
Eq. dl.lOlH it follows that the purity /i ~ Ti'[g 2 ] of a Gaussian state is given by: 



M(<t) = tt" KÍ n / d ín XW\X) = . (2.5) 



'Recali that for every symmetric positive-definite matrix Q 6 M(n, R) the following identity holds 
( d n Xc x p{-^X T Q- 1 X + iA T X} = v / (27r)"Det[Q] exp { - \ A T QA} . 
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Hence a Gaussian state is pure if and only if 

Det[cr] = (2k 2 )" 4 ". (2.6) 

Another remarkable feature of pure Gaussian states is that they are the only pure states endowed with a positive 
Wigner function I30II31I . In order to prové the statement we consider system with only one degree of freedom. The 
extension to n degrees of freedom is straightforward. Let us consider the Husimi function Q(a) = ir^ 1 1 (ip\a) | 2 
where \a) is a coherent state, which is related to the Wigner function as 

Q(a) = - í d 2 (3W(P)exp\-2\a- (3\ 2 } . (2.7) 

Eq. i2.7\ implies that if Q(ao) = for at least one ao then VF(a) must have negative regions, because the 
convolution involves a Gaussian strictly positive integrand. But the only pure states characterized by a strictly 
positive Husimi function turns out to be Gaussian ones. Indeed, consider a genèric pure state expanded in Fock 
basis as \ip) = ^c„|?i) and define the function f(a) = esH (ip\a) — X) c *i"^=r- Clearly, f{a) is an analytic 
function of growth order less than or equal to 2 which will have zeros if and only if Q(a) has zeros. Hence it is 
possible to apply Hadamard's theorem 1321 . which states that any function that is analytic on the complex plane, 
has no zeros, and is restricted in growth to be of order 2 or less must be a Gaussian function. It follows that the 
Q(a) and W^a) functions are Gaussian. 

Gaussian states are particularly important from an applicative point of view because they can be generated 
using only the linear and bilinear interactions introduced in Section fTT^I Indeed, the following theorem, due to 
Williamson, ensures that every covariance matrix (every real symmetric matrix positive definite) can be diagonal- 
ized through a symplectic transformation 1 33 1, 

Theorem 1 (Williamson) Given V € M(2n, E), V T = V, V > there exist S e Sp(2n, E) and D g M(n, E) 
diagonal and positive defined such that: 

V = S^° l) S . (2.8) 

Matrices S and D are unique, up to a permutation of the elements of D. 
Proof. 

By inspection it is straightforward to see that Eq. J2.8I > implies that 

S=(D®D)-^ 2 OV- 1/2 , 
with O orthogonal. Requiring symplecticity to matrix S means that 

OV-^JV-^O? = ( _ ° _! D Q 1 ) , (2.9) 

J being defined in Eq. dl .6i . Since V and J are symmetric and antisymmetric, respectively, it follows that 
V -1 / 2 JV^ 1 / 2 is antisymmetric, hence there exist a unique O such that Eq. j2.9t holds. □ 



The elements dk of D = Diag(c?i, . . . , d n ) are called symplectic eigenvalues and can be calculated from the 
spectrum of iJV, while matrix S is said to perform a symplectic diagonalization. Changing to íí-ordering, i.e. in 
terms of the covariance matrix er defined in Eq. dl.7> . the decomposition J2.8t reads as follows 



S T WS 



(2.10) 



where W — dk I2, I2 being the 2x2 identity matrix. 

The physical statement implied by decompositions d2.8l > and (12. lQi is that every Gaussian state q can be ob- 
tained from a thermal state v, described by a diagonal covariance matrix, by performing the unitary transformation 
Í7s associated to the symplectic matrix S, which in turn can be generated by linear and bilinear interactions. In 
formula, 

q=U s pU í s , (2.11) 
• • Cg> L> n is a product of thermal states of the form dl.9> for each mode, with parameters (3^ 



where v = v\ 
given by 



Pk = In 



~ d k + 1 + (2k 2 ) 
dk - (2k 2 )" : 



(2.12) 
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in terms of the symplectic eigenvalues dk- Correspondingly, the mean number of photons is given by Nk — 
dk — (2«2)~ 2 . The decomposition d2.8l > allows to recast the uncertainty principle (II .81 . which is invariant under 
Sp(2n, R), into the following form 

d k > (2k 2 )- 2 . (2.13) 

Pure Gaussian states are obtained only if v is pure, which means that Vf. = |0) (0|, Vfc, i.e. dk = (2/í2)~ 2 - Hence a 
condition equivalent to Eq. \2.6\ for the purity of a Gaussian state is that its covariance matrix may be written as 

V = (2«; 2 )~ 4n SS T . (2.14) 

Furthermore it is clear from Eq. (12. 1 3I > that pure Gaussian states, for which one has that dk = (2«2) ~ 2 Vfc 2 , are 
minimum uncertainty states with respect to suitable quadratures. 



2.2 Single-mode Gaussian states 

The simplest class of Gaussian states involves a single mode. In this case decomposition ( 12. 1U reads as follows 
1551 : 

g = D(a)S(OvS\Ç)DHa), (2.15) 

where a = Ks(x + iy) (for the rest of the section we put K2 = 2" 1 / 2 ), v is a thermal state with average 
photon number N, Dia) denotes the displacement operator and S(£) with £ = r e lv the squeezing operator. A 
convenient parametrization of Gaussian states can be achieved expressing the covariance matrix u as a function 
of N, r, ip, which have a direct phenomenological interpretation. In fact, following Chapter^ Le. applying the 
phase-space representation of squeezing 1351 1361 . we have that for the state fl2. 15l > the covariance matrix is given 
by cr = S^er^Xç where <r„ is the covariance matrix dl.lObt of a thermal state and the symplectic squeezing 
matrix. The explicit expression of the covariance matrix elements is given by 



2N + 1 
2N + 1 



cosh(2r) + sinh(2r) co$((f) 
cosh(2r) — sinh(2r) cos(<p) 



(2.16a) 
(2.16b) 



2N + 1 

fia = 02i = x sinh(2r) sin(^) , (2.16c) 



and, fromEq. (12. 5> . it follows that [25 37l a = (2ÍV+1) _1 , which means that the purity of a genèric Gaussian state 
depends only on the average number of thermal photons, as one should expect since displacement and squeezing 
are unitary operations hence they do not affect the trace involved in the definition of purity. The same observation 
is vàlid when one considers the von Neumann entropy Sy of a genèric single mode Gaussian state, defined in 
general as 

S v (g) = -TT[glng]. (2.17) 

Indeed, one has 

*w-M^) +ta < w+1) -^ ta (£;)- ta (iT;) ■ <118) 

Eq. ( I2.18> . firstly achieved in Ref. 1 38 1, shows that the von Neumann entropy is a monotonically increasing function 
of the linear entropy (defined as Sl = 1 — p), so that both of them lead to the same characterization of mixedness, 
a fact peculiar of Gaussian states involving only one single mode. 

Examples of the most important families of single mode Gaussian states are immediately derived considering 
Eq. ( I2.15> . Thermal states v are re-gained for a = r = p — 0, coherent states for r = <p = N = 0, while 
squeezed vacuum states are recovered for a = N = 0. For N = we have the vacuum and coherent states 
covariance matrix. The covariance matrix associated with the real squeezed vacuum state is recovered for ip = 
and is given by cr = -|Diag(e~ 2r , e 2r ). 



2.3 Bipartite systems 

Bipartite systems are the simplest scenario where to investigate the fundamental issue of the entanglement in 
quantum information. In order to study the entanglement properties of bipartite Gaussian systems it is very useful 

2 This is an immediate consequence of Eq. 12.1 II together with purity condition Tr[g 2 ] = 1. 
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to introduce normal forms to represent them. This section is for the most part devoted to this purpose. The main 
concept to be introduced in order to derive useful normal forms is that of local equivalence. Two states Q\ and 
Q2 of a bipartite system Ha ® ÏÏb are locally equivalent if there exist two unitary transformations Ua and Ub 
acting on TÍa and Hb respectively, such that Q2 = Ua <8> Ub QiU a (g> U B . The extension to multipartite systems 
is straightforward. 

Let us start introducing the following 

Theorem 2 (Singular vàlues decomposition) Given C € M(n, C) then there exist two unitary matrices U and 
V, such that 

C = UHV , 

£ = Diag(y / pï, . . . , y/p»), where p\ t {k — 1, . . . , n) are the eigenvalues ofthe positive operator C^C. 
Proof. 

Let V be the unitary matrix that diagonalizes C^C; we have 

VC^CV ] = s 2 

C ] C = £ 2 V 

&C = V t SC/ t í7SV, 

and, from the last equality, one has C ] = V f £ Í7 t and C = U £ V, provided that U = CÇE Vy 1 (for a 
detailed proof see Ref. 1 39 1). □ 

Let us now consider a genèric bipartite state 

\C))=J2 c ^h)\^k). (2.19) 

h,k 

Thanks to the singular vàlues decomposition Theorem |2j the coefficients' matrix C can be rewritten as C = 
U £ V, so that 

Chk = 2J U hr &rs V s k , (2.20) 

r,s 

where a rs — y/p^5 r . s - In this way the bipartite state |C)) reads 

\c))=j2vpk~\<t>k)m, (2.2D 

with 

|fo> =J2 Vsk ^' IV-*) =53^1*.)- (2-22) 

s s 

Note that (tph\ipk) = àh,k and (<fih\4>k) — ^h,k- Eq. (12.2 1 i is known as "Schmidt decomposition", while the 
coefficients ^/pk are called "Schmidt coefficients". By construction the latter are unique. 

Let us consider now Gaussian pure states for m+n canonical systems partitioned into two sets A = {Ai , . . . , A m } 
and B = {Bi, . . . , B n } in their Schmidt form 

\Í>)ab = ^2 ^/Pk\<Pk) AWk) b ■ (2.23) 

k 

In general the Schmidt decomposition has an "irreducible" structure: generally speaking, Eq. J2.23Í cannot be 
brought into a simpler form just by means of local transformations on set A and B. In the case of Gaussian 
bipartite systems however a remarkably simpler form can be found 1401 14 1 1 f42l . As a matter of fact, a Gaussian 
pure state \iP)ab rnay always be written as 

W)AB = \Í>i) Àl 5, ■ ■ ■ \&)a.b. \°)a v |0>S. . í 2 - 24 ) 

where A = {A\ . . . , A m } and B = {Bi, . . . , B n } are new sets of modes obtained from A and B respec- 
tively through local linear canonical transformations, the states |i/>fc) are two-mode squeezed states for modes 
fe = 1, . . . , s, for some s < min[m, n] and |0)^ and |0) g are products of vacuum states of the remaining modes. 
In order to prové Eq. (12.24b we consider the partial density matrices obtained from the Schmidt decomposition 

QA = ^2pk\(t)k)(<Pk\, Q B =^2pk\<Pk)(<fk\- (2.25) 

fe k 
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Since qa and qb are Gaussian, they can be brought to Williamson normal form ( 12. lli through local linear canonical 
transformations. Suppose that there are s modes in A and t modes in B with symplectic eigenvalue d ^ (2k2)~ 2 . 
Since the remaining modes factor out from the respective density matrices as projection operators onto the vacuum 
state, we may factor \i/j)ab as \i(>)ab \0) à \®)b wnere (V')ab is a genèric entangled state for modes A\, . . . , A s 
and Bi, . . . , B t . The partial density matrices of the state \iP)ab rnay be written as 



QA 



E: 



n A Tr [. 

where we have used the notation 

CA = 



e ~/3 A -N A 



\nA)(nA\, QB=^2 



Tr [e 



CB 



'Si' ' 



-Pb-Nb 



> C B t 



|ns)(ns| , 



(2.26) 



hence rï^ and ne represent occupation number distributions on each side, and Nb are the number operators, 
and (3a and /3b represent the distribution of thermal parameters. In order to have the same rank and the same 
eigenvalues for the two partial density matrix, as imposed by Schmidt decomposition, there must exist a one-to- 
one pairing between the occupation number distributions Ha and Us, such that (3a ■ ^a = (3b ■ Hb. K turns out 
that this is possible only if s — t and ua — H B , provided that (3a = (3b (for a detailed proof see Ref. 11401 1. Hence, 
reconstructing the Schmidt decomposition of \i(j)ab from qa and qb we see that the form J2.24I is recovered for 
\Í>)ab. 

Let us consider now the case of a genèric bipartite mixed state. Due to the fact that the tensor product structure 
of the Hilbert space translates into a direct sum on the phase space, the genèric covariance matrix of a bipartite 
m + n modes system is a 2m + 2n square matrix which can be written as follows: 



A C 

C T B 



(2.27) 



Here A and B are 2m and 2n covariance matrices associated to the reduced state of system A and B, respectively, 
while the 2m x 2n matrix C describes the correlations between the two subsystems. Applying again the concept 
of local equivalence we can straightforwardly find a normal form for matrix \2.21\ . A genèric local transformation 
Sa © S B , with S A € Sp(2m,R) and S B e Sp(2n,R), acts oncr as 



A^S A AS T A 



B 



Sb B S B , 



C —>■ Sa C S t b 



(2.28) 



Notice that four local invariants [i.e. invariant with respect to transformation belonging to the subgroup Sp(2m, R)® 
Sp(2n,R) C Sp(2m + 2n,R)] can immediately be identified: h = Det [A], I 2 = Det [B], I 3 = Det [C], 
li = Det [tr] . Now, the TheoremHallows to choose Sa and S B such to perform a symplectic diagonalization of 
matrices A and B [see Eq. J2.10H . namely 



S a AS t a = W a = @ d A , k la , S B BS T B = W B = ® d BM l a , 



(2.29) 



k=l 



k=l 



where Wa(b) i s a diagonal matrix. Thus any covariance matrix er of a bipartite m x n system can be brought into 
the form 

( W A E \ 

where E = Sa C S t b . A further simplification concerns the case m = n if we focus our attention on the 



2x2 diagonal blocks of E, which we call E n , with h = 1, . 



Matrices E \ being real 2x2 matrices, 



admit a singular value decomposition by suitable orthogonal (and symplectic) matrices 0\ and O b : E = 
A E h (O b ) t . Such Oa( B \'s transformations do not affect matrices Wa and W B , being their diagonal blocks 
proportional to the identity matrix. Collecting this observations we can write the following normal form for a 
genèric n x n covariance matrix 

f W A È \ 



W B 



(2.31) 



where 



ei,i 






e2,2 



ei,2n-l 
e2,2n-l 



ei,2n \ 



e2n-l,l 62n-l,2 
\ e2n,l e2n,2 



e2n-l,2n-l 

e2n,2n / 



(2.32) 
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Due to the relevance in what will follow, we write explicitly the normal form J2.3H for the case of a 1 x 1 
system. It reads as follows: 

fa ci \ 
a c 2 
ci 6 
c 2 b 



\ 



(2.33) 



/ 



where the correlations a, b, c±, and c 2 are determined by the four local symplectic invariants 1\ = a 2 , 7 2 = b 2 , 
I3 = C1C2 and I4 = (ab — cf )(ab — cl). When a = b the state is called symmetric. The normal form J2.33t allows 
to recast the uncertainty principle dl.8t in a manifestly Sp(2, K.) ® Sp(2, K) invariant form |43 1: 



/1 + h + 2/3 < 



l/(8« 2 ). 



(2.34) 



In order to prové this result it is sufficient to note that it is true for the normal form 12.331 . then the invariance of 
Ineq. 12.341 ensures its validity for every covariance matrix. 

Finally we observe that Ineq. 12.341 can be recast in an even simpler form, in the lxl modes case. Indeed the 
symplectic eigenvalues of a genèric covariance matrix can be computed in terms of the symplectic invariants (we 
putKi = 2- 1 / 2 ) |44|: 



V2d ± = h+I 2 + 2/3 ± V(h+h + 2I 3 ) 2 -4I 4 



1/2 



[defining d\ — d + and d% = cL in the relation ( 12. 8H . The uncertainty relation ( 12.131 1 then reads 



(2.35) 



(2.36) 



Therefore, we may see that purity of two mode states corresponds to I4 = 1/16 and I\ + 1% + 21% = 1/2, where 
the first relation follows from Eq. ( 12. 6> and the second one from the fact that a pure Gaussian state has minimum 
uncertainty. Notice also that bipartite pure states necessarily have a symmetric normal form (i.e., a = b), as can be 
seen by equating the partial entropies Sv calculated from Eq. J2.33Í . 

Eq. ( I2.35> allows also to express the von Neumann entropy Sv of Eq. ( 12.171 in a very simple form. Indeed the 
entropy of a genèric two-mode state is equal to the entropy of the two-mode thermal state obtained from it by a 
symplectic diagonalization, which in turn corresponds to a unitary operation on the level of the density operator g 
hence not affecting the trace appearing in the definition ( I2.17> . Exploiting Eq. ( 12. 1 8i and the additivity of the von 
Neumann entropy for tensor product states, one immediately obtains: 



S v (cr)=f(d + ) + f(d-) , 

where f(d) = {d+ |) In(d +\)-{d-\) ln(d - §). 

A relevant subclass of Gaussian states is constituted by the two-mode squeezed thermal states 3 

Q = S 2 (C) Qv S\ (C) , Qv = v k ®v B , 



(2.37) 



(2.38) 



where Vk, k — A, B, are thermal states with mean photon number Ni and N2 respectively, whose covariance 
matrix rr v is given in Eq. dl . 10bi . Following Chapter[2 i.e., applying the phase-space representation of squeezing, 
we have that for the state d2.38l the covariance matrix is given by cr = X^c v^2£, where is the symplectic 
two-mode squeezing matrix given in Eq. dl .89i . In formula, 



1 

i^ 2 



Al 2 Cfíç 

CRç Bl 2 



being defined in Eq. (II . 891 and 

A = A(r, N U N 2 ) = cosh(2r) + 2 TYj cosh 2 r + 2N 2 sinh 2 r , 
B = B(r, Ni,N 2 ) = cosh(2r) + 2 Ni sinh 2 r + 2 N 2 cosh 2 r , 
C = C(r, Ni,N 2 ) = (1 + Ni + N 2 ) sinh(2r) . 



The TWB state |A)) described in Section fï^^l is recovered when q u is the vacuum state (namely, Ni 
and <p = 0, leading to 

A = B = cosh(2r) , C = sinh(2r) . 



(2.39) 



(2.40a) 
(2.40b) 
(2.40c) 

N 2 = 0) 
(2.41) 



3 For a general parametrization of an arbitrary bipartite Gaussian state, by means of a proper symplectic diagonalization, see Ref. 1441 . 
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2.4 Tripartite systems 

In the last Section of this Chapter we deal with the case of three-mode tripartite systems, i.e. lxlxl systems. 
The genèric covariance matrix of a three-mode system can be written as follows 



0" 11 & 12 Cl3 
O" 21 T22 T23 
""31 ""32 ""33 



(2.42) 



where each <Thk is a real 2x2 matrix. Exploiting the local invariance introduced above and following the strategy 
that led to the normal form J2.3 1 i for the bipartite case, it is possible to find a local invariant form also for matrix 
dZ42l 1451. 

In the following we will consider a local symplectic transformation belonging to the subgroup Sp(2,R) £g> 
Sp(2, R) <g> Sp(2, E) C Sp(6, R), referred to as S = Si ® S 2 © S 3 . The action of S on the covariance matrix er 
is given by 

S h er hh Sl for h = k 
S h cr hk S T k íoxh^k 

Performing, now, a symplectic diagonalization, we can reduce the diagonal blocks as 



(SaS T ) hk = 



(2.43) 



ShChhSÏ — a/j 1 2 



(2.44) 
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Concerning the remaining three blocks, one may follow the procedure that led to Eq. J2.3H . In fact it is still 
possible to find three orthogonal symplectic transformation Oi, O2 and O3 able to put two of the three blocks in 
a diagonal and in a triangular form, respectively, leaving unchanged the diagonal blocks. The covariance matrix 
\2A2\ can then be recast into the following normal form 



(2.45) 



where we can identify 12 independent parameters. 

Three-mode tripartite systems have been studied in different contests, from quantum òptics 1461 1471 . to con- 
densate physics 1 11 1. A study was also performed in which the mode of a vibrational degree of freedom of a 
macroscopic object such as a mirror has been considered |48 1. As examples we consider here the two classes of 
states generated by means of the all optical systems proposed in |46| and |47|. The first generation scheme is a 
very natural and scalable way to produce multimode entanglement using only passive optical elements and single 
squeezers, while the second one is the simplest way to produce three mode entanglement using a single nonlinear 
optical device. They both can be achieved experimentally 1491 1501 . As concern the first class of states, it is gen- 
erated with the aid of three single mode squeezed states combined in a "tritter" (a three mode generalization of a 
beam-splitter). The evolution is then ruled by a sequence of single and two mode quadratic Hamiltonians. As a 
consequence, being generated from vacuum, the three-mode entangled state is Gaussian, and its covariance matrix 
is given by (for the rest of this section we set k 2 — 2~ 1 / 2 ): 
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(2.46) 



where 



1 2 

7^ ± = cosh(2r) ± - sinh(2r) , S = -- sinh(2r) . 

O ó 



(2.47) 



and r is the squeezing parameter (with equal squeezing in all initial modes). 

The second class of tripartite entangled states is generated in a single non linear crystal through a special case 
of Hamiltonian H pq in Eq. ( I1.90> . namely 



71 44 



724°3 



h.c. 



(2.48) 
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which describes two interlinked bilinear interactions taking place among three modes of the radiation field coupled 
with the support of two parametric pumps. It can be realized in media by a suitable configuration exposed 
in Ref. 1 50 1 . The effective coupling constants 7&, k = 1, 2, of the two parametric processes are proportional to 
the nonlinear susceptibilities and the pump intensities. As already seen in Section fï.4.51 if we take the vacuum 



|0) = |0)i ® |0) 2 &> |0) 3 as the initial state, the evolved state |T) 
states of SU(2, 1) and it reads [see Eq. ( I1.94H 



-iH in 



* 1 0) belongs to the class of the coherent 



|T) 



1 



OO 



N 2 



1 + iVi 



p/2 



1 + JVi 



g/2 



3 -í(p<fe+ï<fe) . 



' (p+g)l , , N 



(2.49) 



where Nk(t) = (a^(t) Ofc(í)) represent the average number of photons in the fc-th mode and <pk are phase factors. 
Notice that the latter may be eliminated by proper local unitary transformations U 2 and Ua on modes a 2 and 0,3, 
namely Uk = exp {i(f>ka k ak\, k — 2,3. The symmetry of the Hamiltonian J2.48l > implies that Ni = N 2 + N3, 
where 



Ni 



l7i| 2 



72 1 



si 4 



[cos(Oí) - iy 



A, = 



mi 
fi 2 



■sin 2 (r2í) , 



(2.50) 



with íl = -\/ 1 T2 1 2 — |7i| 2 - Also for this second class, being the initial state Gaussian and the Hamiltonian 
quadratic, the evolved states will be Gaussian. The explícit expression of its covariance matrix reads as follows 
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(2.51) 



where Tk = Nj. + i and 



A k = y/N k (l + N 1 )cos(j)k ■ 
C = \J N 2 N 3 cos(0 2 - 03) , 



B k = ^/Nkiï+m) sin 4> k , 
V = y/N 2 Na sin(0 2 - fa) ■ 



As already noticed, the covariance matrix ( 12.5 11 may be simplified by local transformations setting <f> 2 = <f>a = 0. 
Finally, if the Hamiltonian J2.48Í acts on the thermal state g v = v A ® v B (£> i/ c , with equal mean thermal photon 
number A on each mode, we obtain the following covariance matrix 



Vr.th = (2A + 1)V T 



(2.52) 



Chapter 3 

Separability of Gaussian states 



Entanglement is perhaps the most genuine "quantum" property that a physical system may possess. It occurs 
in composite systems as a consequence of the superposition principle and of the fact that the Hilbert space that 
describes a composite quantum system is the tensor product of the Hilbert spaces associated to each subsystems. 
In particular, if the entangled subsystems are spatially separated nonlocality properties may arise, showing a very 
deep departure from classical physics. 

A non-entangled state is called separable. Considering a bipartite quantum system H = Ha <8> H-b (the 
generalization to multipartite systems is immediate), a separable state g is defined as a convex combination of 
product states, namely I5Ï1: 

fc 

where pk > 0, J2k "Pk = 1' an ^ Pk^' Sk^ b e l° n g to 'tÍAi T^-b respectively. The physical meaning of such a deíï- 
nition is that a separable state can be prepared by means of operations acting on the two subsystems separately (i.e. 
local operations), possibly coordinated by classical communication between the two subsystems 1 The correlations 
present, if any, in a separable state should be attributed to this communication and hence are of purely classical 
origin. As a consequence no Bell inequality can be violated and no enhancement of computational power can be 
expected. 

The separability problem, that is recognizing whether a given state is separable or not, is a challenging question 
still open in quantum information theory. In this chapter a review of the separability criteria developed to date will 
be presented, in particular for what concern Gaussian states. We will profusely use the results obtained in Chapter 
|2|regarding the normal forms in which Gaussian states can be transformed. 



3.1 Bipartite pure states 

Let us start by considering the simplest class of states, for which the separability problem can be straightforwardly 
solved, that is pure bipartite states belonging to a Hilbert space of arbitrary dimension. First of all, recali that such 
states can be transformed by local operations into the normal form given by the Schmidt decomposition J2.23l >. 
namely 

d 

\iP)ab = VPfc \<j>k) A\<Pk) B ■ (3.2) 
k=l 

Therefore, since the Schmidt coefficients are unique, it follows that the Schmidt rank (i.e., the number of Schmidt 
coefficients different from zero) is sufficient to discriminate between separable and entangled states. Indeed, a 
pure bipartite state is separable if and only if its Schmidt rank is equal to 1. On the opposite, a pure state is said 
to be maximally entangled if its Schmidt coefficients are all equal to dr 1 / 2 (up to a phase factor). In order to 
understand this definition, consider the partial traces qa = Trg [g] and qb = Tr^ [g] of the state in Eq. (13.21 . where 
g = \ip){ip\. From Eq. Í3.2Ï it follows that 

QA = y^Pfc \4>k) AA{<t>k\ , QB = *^2pk\(Pk)BB(<fk\ , (3.3) 

fc fc 

hence it is clear that the partial traces of a maximally entangled state are the maximally chaotic states in their 
respective Hilbert space. From Eq. J3.3I > it also follows that the von Neumann entropies (12. 17l > of the partial traces 

'Quantum operations obtained by local actions plus classical communication is usually referred to as LOCC operations. 
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are equal one each other, in formula: 

S$ = S$ = -J2 Pk log dPk . (3.4) 

k 

It is possible to demonstrate that J3.4-Í is the unique measure of entanglement for pure bipartite states 1521 . It 
ranges from 0, for separable states, to 1, for maximally entangled states. 

Let us now address the case we are more interested in, that is infinite dimensional systems. Consider for the 
moment a two-mode bipartite system. Following the definition of maximally entangled states given above, it is 
clear that the twin-beam state (TWB) 

oo 

|A» = yfl - A 2 ]T X n \n,n), (3.5) 

n=0 

where À = tanh r, r being the TWB squeezing parameter, is a maximally entangled state. In fact, its partial traces 
are thermal states, te., the maximally chaotic state of a single-mode continuous variable system, with mean photon 
number equal to the mean photon number in each mode of the TWB, namely (a k a,k) = | A | 2 / (1 — |À| 2 ) = sinh 2 r, 
k = A, B, in the notation of Section fT.4.41 The unique measure of entanglement is then given by Eq. (12. 1 81 . that is 
the von Neumann entropy of a genèric Gaussian single mode state. Remarkably, these observations are sufficient 
to fully characterize the entanglement of any bipartite m x n pure Gaussian state. Indeed in Section l2~3l we have 
demonstrated that such a system can be reduced to the product of TWB states and single mode local state at each 
party. As a consequence the bipartite entanglement of a Gaussian pure state is essentially a 1 x 1 entanglement. 

We mention here that besides the separability criterion given by the Schmidt rank, for pure bipartite system 
another necessary and sufficient condition for the entanglement is provided by the violation of local realism, for 
some suitably chosen Bell inequality |53|. 

3.2 Bipartite mixed states 

The problem of separability shows its complexity as soon as we deal with mixed states. For example, there exist 
states that do not violate any inequality imposed by local realism, but yet cannot be constructed by means of LOCC. 
The first example of such a state was given by Werner 1511 . Despite the efforts, a general solution to the problem of 
separability in the case of an arbitrary mixed state has not been found yet. Most of the criteria proposed so far are 
generally only necessary for separability, even if for some particular classes of states they provide also necessary 
conditions for entanglement. Fortunately, these particular cases are of great relevance in view of the application to 
quantum information, in fact they include 2x2 and 2x3 finite dimensional systems and m x n and lxlxl 
infinite dimensional systems in case of Gaussian states. 

Most of the separability criteria relies on the key observation that separability can be revealed with the aid of 
positive but not completely positive maps. Let us explain these point in more details. Every linear map q ^ C[q], 
in order to be an admissible physical transformation, has to be trace preserving and positive in the sense that it 
maps positive semidefinite operators (statistical operators) again onto positive semidefinite operators. However, a 
physical transformation has not only to be positive: if we apply the transformation only to one part of a composite 
system, and leave the other parts unchanged, then the overall state after the operation has to be described by a pos- 
itive semi-definite operator as well. In other words, all the possible extensions of the map should be positive. Such 
a map is called completely positive (CP). A map which is positive but not CP doesn't correspond to any physical 
operation, nevertheless these maps have become an important tool in the theory of entanglement. The reason for 
this will be clear considering the most prominent example of such a map, transposition (T). Transposition applied 
only to a part of a composite system is called partial transposition (in the following we will use the symbol T with 
a subscript that indicates the subsystem with respect to the transposition is performed). Positivity under partial 
transposition has been introduced in entanglement theory by Peres 1 54 1 as a necessary condition for separability. 
In fact, consider a separable state as defined in Eq. J3 - 1 1 and apply a transposition only to elements of the first 
subsystem A. Then we have: 

k 

Since the transposed matrix (g^) T — (q^)* is non-negative and with unit trace it is a legitimate density matrix 
itself. It follows that none of the eigenvalues of q Ta is negative if g is separable. This criterion is often referred 
to as ppt criterion (positivity under partial transposition). Of course, partial transposition with respect to the 
second subsystem B yields the same result. 

If we consider systems of arbitrary dimensions ppt criterion is not sufficient for separability, but it turns out to 
be necessary and sufficient for systems consisting of two qubits 1 55 1, that is a system described in the Hilbert space 



3.2 Bipartite mixed states 
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H = C 2 C 2 . This is due to the fact that there exist a general necessary and sufficient criterion for separability 
which saying that a state g is separable if and only if for all positive maps £, defined on subsystem TL A , (£<S>I) [g] is 
a semi-positive defined operator 1 55 1 . Due to the limited knowledge about positive maps in arbitrary dimension this 
criterion turn out to be inapplicable in general. Nevertheless, in the case of C 2 it is known that all the positive maps 
can be decomposed as Vi + V%T, where V\, V2 are CP maps |56|. Hence the sufficiency of partial transposition 
criterion for two qubits follows. 

The ppt criterion turns out to hold also for 7i = C 2 ® C 3 systems, but for higher dimensions no criterion 
vàlid for every density operator g is known. Indeed, in general there exist entangled states with positive partial 
transpose, the so called bound entangled states. The first example of such a state was given in Ref. 1571 . 

At first sight it is not clear how the ppt criterion, developed for discrete variable systems, can be translated 
to continuous variables. Furthermore, considering that ppt criterion ceases to be sufficient for separability as the 
dimensions of the system increases, one might expect that it will provide only a necessary condition for separa- 
bility in case of continuous variables. In fact, Simón |43| showed that for arbitrary continuous variable case this 
conjecture is true. However, Simón also demonstrated that for lxl Gaussian states the ppt criterion represents 
also a sufficient condition for separability. 

Simon's approach relies on the observation that transposition translates to mirror reflection in a continuous vari- 
ables scenario. In fact, since density operators are Hermitian, transposition corresponds to complex conjugation. 
Then, by taking into account that complex conjugation corresponds to time reversal of the Schroedinger equation, 
it is clear that, in terms of continuous variables, transposition corresponds to a sign change of the momentum 
variables, i.e. mirror reflection. In formula, 

R — > AR, S^AS, (3.7) 

where 

A = Diag(l, — 1, . . . , 1, —1) , A = l n 8(-1„). (3.8) 

The action of transposition on the covariance matrices V and er of a genèric state reads as follows: V — > A VA 
and er — > AerA, respectively. For a bipartite system H — Ha ®1~Íb partial transposition with respect to system 
A will be performed on the phase space through the action of the matrices A^ = Affil and A A = A © 1, where 
the first factor of the tensor product refers to subsystem A and the second one to B. Following now the strategy 
pursued above in case of discrete variables, a necessary condition for separability is that the partial transposed 
operator is semi-positive definite, which in terms of covariance matrix is now reflected to the following uncertainty 
relation 

A a <tA a + íl > , A A VA A --^J>0. (3.9) 

We may write these conditions also in the equivalent form 

«t>-^Oa, V>^J A) (3.10) 

where íl A = A A ÍIA A and J A = A A JA A . 

Let us consider, in particular, the case of 1 x 1 Gaussian states. We have already seen that, by virtue of the 
normal form ( 12.331 1. relation Eq. A3.9I has the simple local symplectic invariant form given by Eq. ( 12.341 . Recalling 
the definition of the four invariants given in Section l2~3l we have 

h=h, Ï2 = h, h = -h, h=h, (3.11) 

where 1^ are referred to matrix A A er A A , while 1^ to er. Notice that of course these relations would not have 
changed if we had chosen to transpose er with respect to the second subsystem B. Hence, a separable Gaussian 
state must obey not only to Ineq. J2.34l > but also to the same inequality with a minus sign in front of I3. This leads 
to a more restrictive uncertainty relation. Together with J2.34l > they summarize as follows 

Ji+7 2 + 2|7 3 | <&4h + ^j ■ (3.12) 

Moreover, notice that for states with I3 > 0, this relation is subsumed by the physical constrain given by the 
uncertainty relation J2.34Í . Relation Eq. ( I3.12t . being invariant under local symplectic transformations, does not 
depend on the normal form A2.33Í . nevertheless it is worthwhile to rewrite it in case of a correlation matrix given 
in the normal form Eq. d2.33i . In fact, Eq. ( 13. 121 then simplifies to: 

8K?(aò - c\)(ab - cf.) > a 2 + b 2 + 2\c l c 2 \ - i . (3.13) 



28 



Chapter 3: Separability of Gaussian states 



As pointed out in Chapter|2j the uncertainty relation for a covariance matrix can be summarized by a condition 
imposed on its minimum symplectic eigenvalue. Hence, in terms of the symplectic eigenvalues d± of the partially 
transposed covariance matrix the ppt criterion becomes 

d- > (2 Kl )- 2 . (3.14) 

Viewed somewhat differently, the ppt criterion can be translated also in term of expectation vàlues of variances 
of properly chosen operators. In fact, it is equivalent to the statement that for every four-vectors e and e' the 
following Inequality is true: 

,/M2\ ■ 1 



([Aw(e)y) + ([A«(e')] a > > ^(\^e A ,e' A 



\Z(e B ,e' B )\) , (3.15) 
(ei,e 2 ,e 3 ,e 4 ), e A = (ei,e 2 ), e B = (e 3 , e 4 ), 3(e k , e' k ) =e T k Je' k , 



where w(e) — e T R, and we defined e 
with k — A, B. 

Here we have shown that the ppt criterion is necessary for separability, as concern its sufficiency we remand 
to the original paper by Simón |43 1. 

Another necessary and sufficient criterion for the case of two-mode bipartite Gaussian states has been devel- 
oped in Ref. |58 1, following a strategy independent of partial transposition. It relies upon a normal form slightly 
different from [Eq. J2.33H 
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ai - 1/4 a 2 — 1/4 
bi - 1/4 ~~ b 2 - 1/4 ' 

Mil - |da| = V(ai - 1/4) (òi - 1/4) - ^{a 2 - 1/4) (ò 2 - 1/4) . 



(3.17a) 
(3.17b) 



Every two mode covariance matrix can be put in this normal form by combining first a transformation into the 
normal form (12.331 . then two appropriate local squeezing operations. In terms of the elements of (13 . 1 61 the criterion 
reads as follows: 

{{Au,f) + {(Av,f)>^(al + \) , (3.18) 
where Auq indicate the variance of the operator uq, and 



u = a qi 
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lai - 1/4 
bi - 1/4 



(3.19) 



Without the assumption of Gaussian states, an approach based only on the Heisenberg uncertainty relation of 
position and momentum and on the Cauchy-Schwarz inequality, leads to an inequality similar to Eq. (13. 181 . It 
only represents a necessary condition for the separability of arbitrary states, and it states that for any two pairs of 
operators and Bk, with k = x,p, such that [A x , A p ] = [B x , B p ] = i/ (2ftf), if g is separable then 

((Au) 2 ) + ((A,) 2 ) - ' ' ' ' 



> 



(Va > 0) 



where 



aA x T -B x 
a 



aA n ±-B, 



(3.20) 



(3.21) 



In order to demonstrate the equivalence between the necessary condition given by Simon's and Duan et al. 
criteria let us compare Ineqs. ( I3.15> and ( 13.201 . We follow the argument given in Ref. |59|. It is immediate to 
see that when Ineq. ( 13 . 1 5I > is respected then so is Ineq. ( 13.20b . In fact, for any given a it is sufficient to consider 
Ineq. ( I3.15t for e = (a, 0, ±a _1 , 0) and e! = (0, a, 0, =Fa _1 ), and to identify q A = A x , p A = A p and qs = B x , 
Pb = B p . The reverse, can be seen as follows: denoting with e and e' the two vectors for which Ineq. (I3.15> is 
violated, then there exist a, X and a pair of symplectic transformations S A , Sb € Sp(2, R) such that S A (a, 0) = 



Xe A , S A (0,a) = Xe' A S^a^O) = Xe B S B (0,a- 
Q _ X í e i e'i 



Xe' B , namely 

Sb = Xa 



e3 
e 4 



(3.22a) 
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e 3 e' 4 - e 4 e 3 



(3.22b) 
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The existence of À and a is ensuredby the fact that a violation of Ineq. (I3.15> implies that 3(eA, e' A )Z{e.B, e' B ) < 
[otherwise Ineq. ( I3.15t would correspond to the uncertainty principle, consequently it should be respected by any 
state]. Notice now that if Ineq. P.15> is violated for e and e', so is for Ae and Ae', implying that 

([Aui(Ae)] 2 ) + ([Au;(Ae')] 2 ) > -L (|J(Ae A , Xe' A )\ + \Z(Xe B , Xe' B )\) . (3.23) 

By inspection of the left hand side of the Ineq. J3.23b one can identify A x = w(S A (í, 0) t ), A p — w(S A (0, 1) t ), 
B x = ±w(S B (l,0) T ), B p = tw(S*s(0,1) t ), where for simplicity we indicated w(S A )(l,0) T ee w((S a 8 
I)(l, 0, 0, 0) T ) and so on. Being Sa and S B symplectic, the operators introduced satisfy the commutation relation 
[A x ,A p ] = [B X ,B P ] = í/{2k\) and 3(\e A , Xe' A ) — a 2 , 3(\e B , \e' B ) — a~ 2 . Consequently Ineq. ( 13.201 is 
violated. 

Other criteria based on variances of suitable operators can be found in Refs. 1601 1611 l62l . These criteria, 
though only necessary for separability, are worthwhile in view of an experimental implementation. In fact, in order 
to apply criteria ( 13.91 or (13 . 1 8i it is necessary to measure all the entries of the covariance matrix. Although this is 
achievable, e.g. by quantum tomography, it may experimentally demanding. The criteria in Refs. |60 I61ll62ll63l 
allow instead to witness entanglement measuring only the variances of appropriate linear combinations of all the 
modes involved. An experimental implementation of such a criterion can be found in Ref. |49|. 

As an example consider the TWB, whose covariance matrix is given in Eq. ( I2.4H . It is immediate to see that 
the criterion given by Eq. ( 13 . 1 3i implies that sinh 2 (2r) < 0, which is violated for every squeezing parameter r. 
The application of criterion Eq. ( 13 . 1 8i is also straightforward. 

Concerning more than one mode for each party, it is possible to demonstrate that the criterion given by 
Ineq. ( 13.91 gives a necessary and sufficient condition for separability only for the case of 1 x n modes |64|. 
The simplest example where the criterion ceases to be sufficient for separability involves a 2 x 2 system, where 
bound entangled states can be found. For a general n x m Gaussian state there is also a necessary and sufficient 
condition, which states that a covariance matrix er correspond to a separable state if and only if there exist a pair of 
correlation matrices cr A and cr B , relative to subsystems A and B respectively, such that the following inequality 
holds l64l : 

(T>(T A ®<T B . (3.24) 

Unfortunately this criterion is difficult to handle in practice, due to the problem of finding such a pair of correlation 
matrices. A more practical solution has been given in Ref. 1651 . It gives an operational criterion based on a 
nonlinear map, rather on the usual linear partial transposition map, hence independent of ppt criterion. Consider 
a genèric covariance matrix tro, decomposed as usual in the following blocks: 

" " (cl S) ■ <" 5 > 

Define now a sequence of matrices {cr k }, k = 0, . . . , oo, of the form (13.251 . according to the following rule: if cr k 
is not a covariance matrix [i.e., if <r k ^ — í(4k 2 ) -1 í~2] then <Jk+i = 0, otherwise 

A k+1 = B k+1 =A k - $te[D k ] (3.26a) 
C k+1 = -Zm[D k ] (3.26b) 

where D k = C k [B k + i(4K 2 ) -1 f2] -1 Cfc (the inverse should be meant as pseudo-inverse). The importance of 
this sequence is that er is separable if and only if a k is a vàlid separable covariance matrix. Then the necessary 
and sufficient separability criterion states that if, for some k > 1 

1. A k ^ — í(4k 2 ) _1 íí, then er is not separable; 

2. A k — ||Cfc|| op 1 > — í(4ki ) _1 íí, then tro is separable.; 

|| O || p stands for the operator norm of O, i.e. the maximum eigenvalue of yOH). Thus, one just has to iterate 
the map ( 13.261 until he finds that either A k is no longer a covariance matrix or A k — ||Cfc||op 1 is a covariance 
matrix. Moreover it is possible to demonstrate that these conditions occur after a finite number of steps, and that 
in case of a separable er decomposition J3.24Í can be explicitly constructed. We finally mention that recently it 
has been shown 1 66 1 that ppt criterion is necessary and sufficient for a subclass ofmxn Gaussian states, namely 
the bisymmetric ones. The latter are defined asmxn Gaussian states invariant under local mode permutations on 
subsystems A and B. This result is based on the observation that bisymmetric states are locally equivalent to the 
tensor product of a two-mode entangled state and of m + n — 2 uncorrelated single-mode states. 

As for the quantification of entanglement, no fully satisfactory measure is known at present for arbitrary mixed 
two-mode Gaussian states. There are various measures available such as the entanglement of distillation and of 
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formation II67I . They quantify the entanglement of a state in terms of the pure state entanglement that can be dis- 
tilled out of it and the one that is needed to prepare it, respectively. Another computable measure of entanglement 
is the "logarithmic negativity" based on the negativity of the partial transpose |68|. Physically it is related to the 
robustness of the entanglement when the state under consideration evolves in a noisy environment. The negativity 
of a quantum state g is defined as 



where ||0||tr = Tr [v Ot O] stands for the trace norm of an operator O. The quantity Af(g) is equal to | ^ fe Afc|, 
the modulus of the sum of the negative eigenvalues of q Ta , and it quantifies the extent to which q Ta fails to be 
positive. Strictly related to Aí is the logarithmic negativity E^j-, defined as Ej^ = ln || q Ta ||t r - The negativity has 
been proved to be convex and monotone under LOCC |68 1. For two-mode Gaussian states it can be easily shown 
that the negativity is a simple function of cL, which is thus itself an (increasing) entanglement monotone; one has 
infact [68| 

Ej^{ct) = max jo, - ln [(2k 2 ) 2 íL] } . (3.28) 

This is a decreasing function of the smallest partially transposed symplectic eigenvalue d_. Thus, recalling 
Eq. (13. 14-1 . the eigenvalue cL completely qualifies and quantifies the entanglement of a two-mode Gaussian state 
tr. 



3.3 Tïipartite states 

When systems composed by n > 2 parties are considered, the separability issue becomes more involved. An 
immediate observation concerns the fact that situations can occur in which some parties of the total system may 
be entangled one each other but separable from the rest of the system. Thus, a classification of all the possible 
situations must be firstly considered. We adopt the classification introduced in Ref. |69| which exploits all the 
possible ways to group the n parties into m < n subsets, which are then themselves considered each as a single 
party. Now, it has to be determined whether the resulting m-party state can be written as a mixture of m-party 
product states. The complete record of the m-separability of all these states then characterizes the entanglement 
of the n-party state. Let us investigate in particular the case we are more interested in, that is tripartite systems. 
For these systems, we need to consider four cases, namely the three bipartite cases in which AB, AC, or BC are 
grouped together, and the tripartite case in which all A, B, and C are separate. In total, we have the following five 
different entanglement classes: 

Class 1 (Fully inseparable states or genuinely entangled states) States which are not separable for any grouping 
of the parties. 

Class 2 (1 -party biseparable states) States which are separable if two of the parties are grouped together, but 

(r) (st) 

inseparable with respect to the other groupings. In general, such a state can be written as Y2hPh Sh ® Bh 
for one party r. 

Class 3 (2-party biseparable states) States which are separable with respect to two of the three bipartite splits but 
inseparable with respect to the third, i.e. they can be written as ^ h Ph É?l <8> for two parties r. 

Class 4 (3-party biseparable states) States which are separable with respect to all three bipartite splits but cannot 
be written as a mixture of tripartite product states. 

Class 5 (fully separable) States that can be written as a mixture of tripartite product states, J2hP h &h ® 8h ® 

(C) 

Qh ■ 

Needless to say, the most interesting class is the first one. In fact fully inseparable states are necessary to implement 
genuinely multipartite quantum information protocols able to increase the performances with respect to classical 
ones |46|. 

In general, it is hard to identify the class to which a given state belong. The problem arises even in case of 
pure states, because a Schmidt decomposition doesn't exist in general. The state vector then cannot be written as 
a single sum over orthonormal basis state. Concerning discrete variable systems, it is known that are only two 
inequivalent classes of pure fully inseparable three-qubits states, namely the GHZ |70| and the W states [71 1 



GHZ) = (|000) + |Hl))/\/2 |W) = (|100> + |010) + |0G1»/V3 



(3.29) 
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In other words, any pure fully inseparable three-qubit state can be transformed via stochastic LOCC (where 
stochastic means that the transformation occurs with non-zero probability) to either the GHZ or the W state. 
Hence a satisfactory knowledge for this case has been achieved. 

When arbitrary mixed states are considered there is no general necessary and sufficient criterion to ensure 
genuine entanglement. The diffïcult of the subject is well exemplified if one exploits the issue of nonlocality 
considering multi-party Bell inequalities. Indeed, violations of such inequalities ensures only that the state under 
investigation is partially entangled. Reversely, fully inseparable states do not necessarily violate multi-party Bell 
inequalities. As an example, the pure genuinely n-party entangled state 

|V>) = cosa|0...0) +sina|l...l) (3.30) 

for sina < 2~(™~ 1 )/ 2 does not violate any n-party Bell inequality \72\, if n is odd, and does not violate Mermin- 
Klyshko inequalities 1731 1741 1? 5 1 for any n. 

Nevertheless, regarding the case of tripartite three-mode Gaussian states, the separability has been completely 
solved. Extending Simon's ppt approach Giedke et al. |76 1 gave a simple criterion that allows to determine which 
class a given state belong to. Hence genuine entanglement, if present, can be unambiguously identified. Observing 
that for these systems the only partially separable forms are those with a bipartite splitting of 1 x 2 modes, it 
follows that the ppt criterion is necessary and sufficient. We have the following equivalences: 

<t > — — -t;ÇIa , c It —-r^^B i o > — — -k£Ic Class 1 (3.31a) 

Anf 4rt| 4/tj 

a> % -xtí A ,<r> Kà B ,(T> ^íí c <^Class2 (3.31b) 

4/q 4/q 4/q 

er > %Sl A , cr > , er > ^fi c ^ Class 3 (3.31c) 

4/í| \k\ 4k| 

rr > --^Üa , cr > — t^Ob , er > --^ft c & Class 4 or 5 (3.31d) 

4kj \k\ 4kj 

Analogue formulas may be written for the covariance matrix V. Notice that in classes 2 and 3 all the permutations 
of the indices A, B, and C must be considered. Classes 4 and 5 cannot be distinguished via the ppt criterion. 
An additional criterion has been given in Ref. 1761 to distinguish between these two classes. It is based on the 
consideration that necessary and sufficient for full separability is the existence of three single mode covariance 
matrices cr\, cr B , <Tç such that 

er > a\ © a x B ® <r l c . (3.32) 

Obviously, for the identification of fully inseparable states, only class 1 has to be distinguished from the rest, thus 
the ppt criterion alone suffices. 

As examples, consider the states given in Section l2~4l The separability issue of state J2.46Í has been addressed 
in Refs. 1.431 1761 . In particular, in Ref. |76| the authors analyzed a generalization of state J2.46Í . in which some 
noise has been added. The state considered is described by the covariance matrix er 3 i(L1 = 0-3 + ^1. Depending 
on the value of the squeezing parameter r and of the noise coefficient /i, <x 3 M belongs either to class 1, 4 or 5. 
If fj, = the state is fully inseparable for any value of r. In fact, applying the ppt criterion we find that matrix 
er 3 + i(4:Ki)~ 1 riA always has a negative minimum eigenvalue A m in given by (k\ = 2 -1 / 2 ) 

A min = cosh(2 r) v/3 + 3 cosh(4 r) + 8 y/2 sinh(2 r) . (3.33) 

V6 

From the symmetry of the state full inseparability follows. On the contrary, if u > 1 then Ineq. ( I3.32l > is satisfied 
with <j\ = <j l B = <Tç = 5!) hence the state is separable. A detailed inspection considering a fixed squeezing 
parameter r shows that two threshold value of the noise u , /ii can be identified, such that <T3 iM is fully inseparable 
for fi < /ió and separable for fi > fi\. When fio < fi < fií it belongs to class 4, hence it is an example of a bound 
entangled state, having every partial transpose positive, nevertheless being inseparable. 

Let us focus now on state J2.52I . The symmetry of this state under the exchange of modes 02 and a 3 allows 
to study the separability problem only for modes ci\ and 02- Furthermore, as already pointed out, we can set 
<p2 = 03 = without affecting the entanglement properties of the state under investigation. Concerning the 
first mode, from an explícit calculation of the minimum eigenvalue of matrix Vt,üi — | J ' A (we consider again 
m = 2~ 1/2 ) it follows that 



Ai nin = N + (1 + 2A0 [iVi - ^1(^1+!) 



(3.34) 
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As a consequence mode ai is separable from the others when 

N > Ni + s/ïh(Ph+T) . (3.35) 

Calculating the characteristic polynomial of matrix Vr,th — \Jb one deals with the following pair of cúbic 
polynomials 

çi(A, Ni, N 2 , N) = A 3 - 2 [2(1 + JVi) + N(3 + 4JVi)] A 2 

+ 4 [1 + N 2 + 2N 3 + N(4 + AN 2 + 6N 3 + N(3 + 4iVi))] A 

-8N[1 + N 2 + N(2 + N + 2N 2 )} , (3.36) 

q 2 (X, N 1} N 2 , N) = A 3 - 2 [1 + 27Vj + iV(3 + 4^)] A 2 

+ A[N 2 + 2N{1 + Nt) + N 2 {3 + 4iVi)] A 

- 8(1 + N) (N 2 - 2N 2 - 2NN 2 ) . (3.37) 

While the first polynomials admits only positive roots, the second one shows a negative root under a certain 
threshold. It is possible to summarize the three separability thresholds of the three modes involved in the following 
inequalities 

N > N k + ^N k {N k + l) . (3.38) 

If the inequality d3 .38i is satisfied for a given k, then mode a k is separable. Clearly, it follows that the state |T) 
evolved from vacuum (i.e., N — 0) is fully inseparable. 

When one deals with more than three parties and modes the separability issue becomes more involved, even 
remaining in the frameworkof Gaussian states. As an example, considerthe case of four parties and modes, labeled 
by A, B, C and D. The one-mode bipartite splittings can still be tested via the ppt criterion, involving 1x3 
modes forms. In the Gaussian language it is necessary and sufficient to consider whether <x ^ — í(An 2 )~ 1 rLs (for 
S = {A, B, C, D}). However, also bipartite splittings of the 2x2 mode type must be taken into account. We 
have already mentioned above that for this case the ppt criterion ceases to be sufficient for separability. Hence 
to rule out the possibility of bound entanglement one have to rely on the operational criteria given in Ref. |65|. 
In general, in order to confirm genuine n-party entanglement, one has to rule out any possible partially separable 
form. In principle, this can be accomplished by considering all possible bipartite splittings and applying either 
the ppt criterion or the criterion from |65|. Although a full theoretical characterization including criteria for 
entanglement classification has not been considered yet for more than three parties and modes, the presence of 
genuine multipartite entanglement can be confirmed, once the complete correlation matrix of the state is given. 



Chapter 4 

Gaussian states in noisy channels 



In this Chapter we address the evolution of a n-mode Gaussian state in noisy channel where both dissipation and 
noise, thermal noise as well as phase-sensitive ("squeezed") noise, are present. At first we focus our attention on 
the evolution of a single mode of radiation. Then we extend our analysis to the evolution of a n-mode state, which 
will be treated as the evolution in a global channel made of n non interacting different channels. For the single 
mode case a thorough analysis may be found in |77 1. 

4.1 Master equation and Fokker-Planck equation 

The propagation of a mode of radiation (the system) in a noisy channel may be described as the interaction of the 
mode of interest with a reservoir (bath) made of large number of external modes, which may be the modes of the 
free field or the phonon modes of a sòlid. We denote by bj such mode operators and assume a weak coupling gj 
between the system and the bath modes. Interaction Hamiltonian is written as Hj = aB^ (t)e~ luJat + a) B(t)e ÍUJat , 
B(t) = J2t 9kbke~ ÍUJht being the collective mode of the bath and u a the frequency of the system. The global 
density matrix R t = R, describing both the system and the bath at time t, evolves, in the interaction picture, 
according to the equation R = i[R, Hj], while the reduced density operator g for the system only is obtained by 
partial trace over the bath degrees of freedom. Upon a perturbative expansion to second order and assuming a 
Markovian bath, i.e. (b(u>h) b(oJk))R — M5(2ui a — ojh — ut) and {b\u>h) b(uik))n = NS(u>h — ^fc) the dynamics 
of the reduced density matrix is described by the following Master equation 

g = -{(N + l)£[a\ + N£[ai} - M*V[a] - M2?[a t ]} g, (4.1) 

where T is the overall damping rate, while N S M and M £ C represent the effective photons number and the 
squeezing parameter of the bath respectively. £[0]p = 20gO^ — O^Og — gO^O and 2?[0]p = 20gO — OOg — 
gOO are Lindblad superoperators. The terms proportional to C[a] and to £[a^] describe losses and linear, phase- 
insensitive, amplification processes, respectively, while the terms proportional to V[a] and ï?[a^] describe phase 
dependent fluctuations. The positivity of the density matrix imposes the constraint |M| 2 < ^(A^ + 1). At thermal 
equilibrium, i.e. for M ~ 0, N coincides with the average number of thermal photons in the bath at frequency oj a . 

4.1.1 Single-mode Gaussian states in noisy channels 

Let us now focus on Gaussian states and start with single mode states. The first step is to transform the Master 
equation ( 14. 1> into a Fokker-Planck equation for the Wigner function. Thanks to Eq. dl.l03> it is straightforward 
to verify the correspondence 

ag ^ {a + \d a ,)W[g]{a) , a^g (a* - ±d a )W[g](a) , (4.2a) 
ga ^ (a - ld a *)W[g](a) , ga^ -> (q* + \d a )W[g]{a) . (4.2b) 

Eqs. ( I4.2> . together with the composition rules L[Oi0 2 ] = L[0{\ L[0 2 ] and R[Oi0 2 ] = R[0 2 ] R[0{\, where 
L and R denote action on the density matrix from the left and from the right respectively, allows to evaluate the 
differential representation of superoperators in Eq. (14. II . We have 

C[a]g -> [3 a a + 8 a *a* + d 2 aa ,] W[g\{a) , (4.3a) 
C[o!]q - [d a a + d a *a* - d 2 aa ,] W[g](a) , (4.3b) 
V[a]g - -dl. a ,W[g](a) , V[ct]g -> -d 2 aa W[g](a) . (4.3c) 
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From now on we put W(a) = W[g](a). In this way, the Master equation ( 14. 1 i transforms into the following 
Fokker-Planck equation for the Wigner function 

d t W(a) = ^{d a a + d a *a* + (2N + l)d 2 aa , + M*d 2 a , a , + Md 2 aa }w{a) . (4.4) 

Passing to Cartesian coordinates a — «2 (x + íy), d a = {2K2)~ x {d x — id y ), we have 



d a a + d a *a* = d x x + d y y , d 2 aa , = (d% x + d 2 y ) , 

d 2 -—(d 2 -2id 2 -d 2 ) d 2 = — (d 2 + 2id 2 -d 2 ) 

w aa \ u xx """ay w yy I ) w a*a* | K 2 ^ u T yy ' 

and Eq. J4.4I rewrites as 

dtW(x, y) = Ud x x + d v y + ^N + l)(d 2 xx + d 2 yy ) 



v 2 

or, in a more compact form, as 



^(m S [M]{dl x -dl y ) + 2^m[M]dl y ) \]VÍ,r.u). (4.5) 



d t W(X) = E ( d T x X + d x a x d x ) W(X) , (4.6) 



2 

where X = (x, y) T , d x = (d x , d y ) T , and we introduced the diffusion matrix 

(l + N ) + 3?m[M] 
2k z 2 y 3m[M] (i + N) - 3te[M] 



(4.7) 



The diffusion matrix is determined only by the bath parameters and, as we will see, represents the asymptotic 
covariance matrix when the initial state is Gaussian. 

The Wigner function at time t, Wt(X), i.e. the general solution of Eq. J4.6I > can be expressed as the following 
convolution 

W t (X)= í d 2 ZGt(X\Z)W (Z) (4.8) 

JR 2 

where W (X) is the initial Wigner function and the propagator G t (X\Z) is given by 

exp \-\{X - e-^ Tt Zy Xt\X - er^Z)} 
GAX\Z) = ^ , L , (4.9) 

with Sí = (1 — e~ rí ) (Joo- The solution (14. 8> holds for any initial Wq{X). For an initial Gaussian state, since 
the propagator is Gaussian, Eq. (14. 8> says that an initial Gaussian state mantains its character at any time. This 
fact is usually summarized saying that the Master equation ( 14. 1> induces a Gaussian map on the density matrix of 
a single-mode. 

From now on, we put «2 = 1 and consider an initial Gaussian state. Using Eq. ( 14. 6> . the evolution of X is 
given by 

X~= { d 2 XXd t W(X) 

JR 2 

= £/ d 2 XXd x XW(X) + J - [ d 2 XXd T x a 00 d x W{X), (4.10) 

* JR 2 * JR 2 

The first integral is easily evaluated by parts, leading to — |r, while the second gives no contribution. Eq. ( 14.1 Ol i 
thus becomes 

X~ = ~X, (4.11) 

i.e. X is damped to zero. 



4.1 Master equation and Fokker-Planck equation 
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Now we address the evolution of the covariance matrix er. Since 

(T ii =ï 2 -2m, à yy =y 2 -2yy, & xy = xy - x y - x y , (4.12) 

we should only evaluate x 2 , y 2 and xy. These evolve as follows 

/ x 2 \ „ Ix 2 

2 = / d 2 X\ y 2 \d t W(X) = -2\ lp |+2l <7&> I, (4.13) 



y 

\ xy 



xy 



X 2 \ 




(° 




+ 2 




xy ) 







where, in solving Eq. (14. 1 31 . we have substituted (14. 6> and integrated by parts. Therefore, the evolution equation 
for er simply reads 

à = -T (cr - (Toc) , (4.14) 

which yields 

a(t) = e~ rt <x(0) + (1 - e- rt ) cr^ , (4.15) 



in agreement with Eqs. (I4.8Í and ( I4.9Í . Eq. J4. 15i says that the evolution imposed by the Master equation is a 
Gaussian map with er^ as asymptotic covariance matrix. <r(t) satisfies the uncertainty relation (II .81 iff these are 
satisfied by both er^ and er(0). 



4.1.2 n-mode Gaussian states in noisy channels 

In this Section we extend the above results to the evolution of an arbitrary n-mode Gaussian state in noisy channels. 
We assume no correlations among noise in the different channels. Therefore, the dynamics is governed by the 
Master equation 

g = J2 -f{( N h + l)C[a h ] + N h C[a{] - M* h V[a h ] - M h V[a{]} g , (4.16) 

h=l 

where Nh and Mh have the same meaning as in Eq. ( 14. li and each channel has a damping rate IV The positivity 
of the density matrix imposes the constraint |i\l/i| 2 < Nh(Nh + 1) V/i. At thermal equilibrium, i.e. for Mh = 0, 
the parameter Nh coincides with the mean number of thermal photons in the channel h. 

As for the single mode case, we can convert the Master equation ( I4.16> into a Fokker-Planck equation for the 
Wigner function. In compact notation we have 

d t W(X) = ^{d^TX + d^T cr 00 dx)w(X) , (4.17) 

with T = 0' h l =1 T h l 2 - Eq. fyj} is formally identical to Eq. but now X = [xx, j/i, ... , x n , y n ) T , dx = 
(d Xl , d yi , . . . , d Xn , d Vn ) T and the diffusion matrix is given by the direct sum (7^ = ©^ =1 <Th ,oo where 

J_ í (l + N h ) + ^e[M h ] Sm[M h ] 
<T?1 ' DO ~ 2«i ^ Sm[M h ] (| + N h ) - MM h . 

is the asymptotic covariance matrix of the /i-th channel. The general solution of ( 14. 17t is an immediate generaliza- 
tion of d4.8i and therefore, also for the n-mode case, we have that Gaussian states remains Gaussian at any time. 
For an initial n-mode Gaussian state of the form \2.2\ the Fokker-Planck equation ( I4.17t corresponds to a set of 
decoupled equations for the second moments that can be solved as for the single mode case. Notice that the drift 
term always damps to the first statistical moments, i.e. 

n 

X(t) = <Gl /2 X{0) with <G t = e~ Tht 1 2 • (4.19) 

h=l 

The evolution imposed by the Master equation preserves the Gaussian character of the states. The covariance 
matrix at time t is given by 

o-(í) = G t 1/2 o-(O) G t 1/2 + (1 - Gt) <Too . (4.20) 

Eq. ( 14.201 describes the evolution of an initial Gaussian state er(0) into the Gaussian environment cíoq. Since 
Eq. ( 14.201 is formally similar to Eq. ( 14. 151 . the considerations we made about the evolved covariance matrix for 
the single mode also hold for the n-mode state. 



(4.18) 
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Sometimes it is useful to describe a system by means of its characteristic function x(A). Since the Wigner 
function is defined as the FT of x(A) we have that x(A) obeys 

n 

d t x(A) = - 2 ( AT]r ^A + A T r ct^A) X (A) , (4.21) 

h=l 

with 9a = (d ai ,3b!,..., d a „ , db n ) T and ò-qo = (/Í2/K3) 2 cr^ Finally, Eq. J4.2 li can be integrated leading to the 
solution 

X (A) = A T G t 1/2 exp {-i A T (1 - <E t ) A} Xo (A) , (4.22) 

Xo(A) being the initial characteristic function. Eq. 14.221 confirms that the evolution imposed by the Master 
equation maintains the Gaussian character of states. 



4.2 Gaussian noise 

In this Section we address the noise described by the map 

Qa{q) = j d 2 "7 ex P{-J rA '^} Dh) ^t (7) , (4.23) 

A being the covariance matrix characterizing the noise and D(~f) is the displacement operator (ll.32> . The map 
Ga is usually referred to as Gaussian noise. Using Eq. dl . 1 1 1 al . the characteristic function of the state Ga{q) is 
given by 

x IGa(q)] (A) = / d 2 " 7 eM ' 1 / l^ ï} eV ' À ~ ÀtT x[ff](A) 

= X[0\W exp{-A T AA} , (4.24) 
whereas, thanks to Eq. dl . 1 1 lbl . its Wigner function reads 

W [0 A(í) ] („) = / ^ 7 exp{-( 7 -arA-y 7 -a)} ^ ^ 

i.e. a Gaussian convolution of the original Wigner function. The average number of photons of a state passing 
through a Gaussian noise channel is obtained using Eq. (11.1 21 

00 00 

J2( a l ak ^Me) = 5]<4 a fc} e + ^Dct[A] , (4.26) 

k=l k=l 

Sfe ( a Í a k) g being the average number of photons in the absence of noise. 

When W[g](ot) itself describes a Gaussian state, i.e. has the form Eq. i l2.lt . then W [Ga{q)] ( a ) is Gaussian 
too, with covariance matrix 

o-gn = cr Q + ÍA. (4.27) 

The Gaussian noise map can be seen as the solution of the Master equation J4. 16i in the limit of large thermal 
noise and short interaction time. In order to derive this result, let us consider Th = T, Nh = N and M = Wi in 
the Eq. J4. 17b . Then, in the limit Tt -C 1 Eq. (14.201 reads cr(t) = cr a + Tt (cr^ — cr a ), which, assuming N ^> 1, 
namely considering cr a negligible with respect to cr^, reduces to 

o-(í) =a a +Tt(T 00 . (4.28) 

By comparing d4.27l > and ( 14.281 . one has that, for N ^> 1 and Tt <C 1, the evolution imposed by the Master 
equation d4. 16i is equivalent to an overall Gaussian noise with covariance matrix given by 

A = 2Tt(T oc . (4.29) 



4.3 Single-mode Gaussian states 

In this Section we address the evolution of single-mode Gaussian states in a noisy channel described by the Master 
equation ( 14. li . In particular, in the following two Sections, we analyze the evolution of purity and nonclassicality 
as a function of time and noise parameters. 
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4.3.1 Evolution of purity 

As we have seen in Chapter[2] the purity /i of a quantum state g is defined as /i = Tr[g 2 ]. For continuous variable 
systems one has < fi < 1. Since /i is a nonlinear function of the density matrix it cannot be the expectation value 
of an observable quantity. On the other hand, if collective measurements on two copies of the state are possible, 
then the purity may be directly measured [78 1. For instance, collective measurement of overlap and fidelity have 
been experimentally realized for qubits encoded into polarization states of photons |79 80 1. 

Purity fi can be easily computed for Gaussian states. In fact, using Eqs. ( 11.1011 and Eq. i2.5\ . for an n-mode 
Gaussian state we have fi(cr) = [(2ft2) 2ll \/Det M] -1 - 

Here we focus our attention only on the evolution of the purity in the case of a single mode Gaussian state of 
1811 : the purity of a two-mode Gaussian state is studied in Ref. |82|. We assume an initial state with zero first 
moments X = 0, i.e. a state of the form S(ro, (po) v S^(ro, (fo). Using ( 14. 1H we conclude that X t = Ví 
and that three parameters are enough to describe the state at any time. These may be either the three independent 
elements of the covariance matrix, the three parameters r(t), ip(t) and N(t), or as it will be the following the three 
parameters r(í), ip(t) and fi{t). 

Let us first consider the case M — 0, for which the initial state is damped toward a thermal state with mean 
photon number N 1351 1251 . In this case ip is constant in time and does not enter in the expression of fi. The 
quantities fi(t) and r(t) in Eqs. ( 14.3 lï solve the following system of coupled equations 



fi = r [n 



fi 2 cosh(2r) 



LJL 

2 floc 



sinh(2r) , 



(4.30) 



which, in turn, can be directly found working out the bàsic evolution equation fi = 2Tï[q q] as a phase-space 
integral; /ioo is defined as fioo = (2N + It is easy to see that, as t — > oo, fi(t) — > fioo and r(t) — > 0, as 

one expects, since the channel damps (pumps) the initial state to a thermal state with mean photon number N. 
Therefore, the only constant solution of Eq. J4.30I is /i = fi^, r = 0, i.e. only initial non-squeezed states are left 
unchanged by the evolution in the noisy channel. The general solution of (I4.30Í is given by 



(1 



,- rt ) 2 + e- 2Tt + ^ e- rt (1 - e- rí ) cosh(2r ) 



-1/2 



(4.31) 



with 



cosh[2r(í)] = fi(t) 



-Vi 



_ rt cosh(2r ) ^ 



/'o 



(4.32) 



Eq. d4 . 3 1 i shows that /i(í) is a decreasing function of tq: in a non-squeezed channel (M = 0), a squeezed state 
decoheres more rapidly than a non-squeezed one. The optimal evolution for the purity, obtained letting r = in 
Eq. (14.3 II . reads 

fi{t) = _ r r . (4.33) 

Obviously, fi(t) is not necessarily a decreasing function of time: if fiQ < fi^ then the initial state will undergo a 
certain amount of purification, asymptotically reaching the value fioo which characterizes the channel. In addition, 
fi(t) is not a monotonic function for any choice of the initial condi tions. Letting /i = in Eq. J4.30I . and exploiting 
Eqs. <E3p and P3H . one finds the following condition for the appearance of a zero of fi at finite positive times: 
cosh(2ro) > max[/io//ioo, Moo/mo]- If this condition is satisfied, then fi(t) shows a local minimum. 

Let us now consider the case M ^ 0, corresponding to a squeezed thermal bath. The general solution for 
purity can be written as 



= Mol — (1 - e 



-rt 



2 + 



2^ (i _ e -rí) e -rt 



smh(2roo) sinh(2r ) 003(2^00 — 2^0) + cosh(2roo) cosh(2r ) 



-1/2 



(4.34) 



where we have already inserted the asymptotic vàlues of the parameters fi, r and tp, i.e. 

Moo = [(2iV+l) 2 -4|M| 2 ]- 1/2 



cosh(2r co ) = ^l + AfiUM\ 2 . 
Q?m[M] 



tan(2(p 00 ) 



(4.35a) 
(4.35b) 

(4.35c) 
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These vàlues characterize the squeezed channel. We see from Eq. 14.341 that /i(í) is a monotonically decreasing 
function of the factor cos(2ip oci — 2(fo), which gives the only dependence on the initial phase tp Q of the squeezing. 
Thus, for any given tp^ characterizing the squeezing of the bath, (po = tp^ + ir/2 is the most favorable value of 
the initial angle of squeezing, i.e. the one which allows the maximum purity at a given time. For such a choice, 
[i(t) reduces to 

M (í) = Mo (iÍ(l-e- rí ) 2 + e- 2rt + 2-^-cosh(2roo - 2r ) (l - e'") e - rt l . (4.36) 

This is a decreasing function of the factor cosh(2roo — 2ro), so that the maximum value of the purity at a given 
time is achieved for the choice ro = roo , and the evolution of the purity of a squeezed state in a squeezed channel 
is identical to the evolution of the purity of a non-squeezed state in a non-squeezed channel. 

In conclusion, for a general channel characterized by arbitrary /loo, r m , (foo and T, the initial Gaussian state 
for which purity is best preserved in time must have a squeezing parameter r = r M and a squeezing angle 
<Po = foo + t/2, i.e. it must be anti-squeezed (orthogonally squeezed) with respect to the bath. The net effect for 
the evolution of the purity is that the two orthogonal squeezings of the initial state and of the bath cancel each other 
exactly, thus reproducing the optimal purity evolution of an initial non-squeezed coherent state in a non-squeezed 
thermal bath. 



4.3.2 Evolution of nonclassicality 



As a measure of nonclassicality of the quantum state g, the quantity r, referred to as nonclassical depth, has been 
proposed in Ref. 1 83 1 

t = , (4.37) 



where s is the maximum s for which the generalized quasiprobability function 

j2n , 



f d ln A 

(X)= -^ rX (A)exp{ l A-X + SK3 |A| 2 } , 

JR 2 " T 



(4.38) 



is a probability distribution, i.e. positive semidefinite and non singular. As one should expect, r = 1 for number 
states and r = for coherent states. The nonclassical depth can be interpreted as the minimum number of thermal 
photons which has to be added to a quantum state in order to erase all the 'quantum features' of the state. 1 While 
quite effective in detecting nonclassicality of states, the nonclassical depth is not easily evaluated for relevant 
quantum states, with the major exception of Gaussian states. In fact, for a Gaussian state characterized by a 
covariance matrix er, the explícit expression for the nonclassical depth reads 



l-2u 



.0 



(4.39) 



u being the minimum of the eigenvalues of cr. In the case of a single mode Gaussian state, this smallest eigenvalue 
turns out to be simply u = e~ 2r / \i 1771 . In this way, thanks to Eq. ( I4.39> . we obtain the following expression for 
the nonclassical depth: 



-2r 



1 - 



I' 



.0 



Therefore, we define the quantity re(í) as 

«(*) 

the time evolution of the nonclassical depth is given by 



cosh(2r ) £ _ rt ^ cosh(2r 00 ) . _ ^_ rt ^ 

A*0 Moo 



r(t) 



1 - K(t) + y/ K{t) 2 - fl(t)- 



(4.40) 



(4.41) 



(4.42) 



which increases with both /i(í) and n(t). The phase maximizing r(í) at any time is again cpo = ip^ + ir/2, 
as for the purity. The maximization of r(í) in terms of the other parameters of the initial state is the result of 
the competition of two different effects: on the one hand a squeezing parameter ro matching the squeezing 
maximizes the purity thus delaying the decrease of r(t); on the other hand, a bigger value of ro obviously yields 
a greater initial r(0). Numerical analysis unambiguously shows |77| that, in non-squeezed baths, the nonclassical 
depth increases with increasing squeezing ro and purity fiQ, as one should expect. 

'This statement can be made more rigorous by assuming that a given state owns 'quantum features' if and only if its P-representation is 
more singular than a delta function (which is the case for coherent states) 1831 . 
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4.4 Two-mode Gaussian states 

In this Section we address the separability of two-mode Gaussian states propagating in a noisy channel. In particu- 
lar, we consider the effect of noise on the twin-beam state of two modes of radiation | A) ) = Vl — A 2 ^2 h ^ h \h)\h), 
|A| < 1 <, A € R, whose Gaussian Wigner function has the form d2.2> with n = 2, X = 0, and the covariance 
matrix given by 

1 / cosh(2r)l 2 sinh(2r)<T 3 \ 
o-twb = -t^t , (4.43) 

^ sinh(2r) er 3 cosh(2r) 1 2 j 

where er 3 = Diag(l, — 1) is a Pauli matrix, and r = tanh -1 A the squeezing parameter of the TWB. 



4.4.1 Separability thresholds 



The Wigner function of a TWB is Gaussian and the evolution in a noisy channel preserves such character, as we 
have seen in Section l4~Tl Therefore, we are able to characterize the entanglement at any time and find conditions 
to preserve it after a given propagation time or length. As we have seen in Chapter|3]a Gaussian state is separable 
iff its covariance matrix satisfies the relation S = <x + ï(4/í 2 ) _1 íí j 4 > 0. Let us now focus on the separability of 
the TWB evolving in generalized Gaussian noisy channels described by the Master equation J4. 16i . The evolved 
covariance matrix is simply given by Eq. (14.201 with <r(0) = <r TWB , and, assuming M as real, its explicit expression 
is 



^ = 2^ 
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(4.44) 



where 
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(4.45) 



Dl(t) 



1 + 2N± 2M 



-rt\ 



(4.46) 



In deriving Eq. 14.441 we have put M% 
when 



M-2 = M and Ni = N2 = N. The conditions ( 13.101 1 are then satisfied 



S? S 2 > — . 



1 

16 



v 2 v 2 > 

^2 ^3 > 16 » 



(4.47) 



which do not depend on the sign of M. 

From now on we put k-2 = 1. If we assume the environment as composed by a set of harmònic oscillators 
excited in a squeezed-thermal state of the form g — v S^(í;), we can rewrite the parameters N and M in 
terms of the squeezing and thermal number of photons N s = sinh 2 £ and Nth., respectively. In this way we get 

M = (1 + 2 N th ) ^N s (l + N s ) , N = N th + N s (l + 2 N th ) . (4.48) 



Now, by solving inequalities (14.471 with respect to time t, we find that the two-mode state becomes separable for 
t > t B , where the threshold time t s = í s (r, T, iVth> N s ) is given by |85 1 
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1 + 2N th 




N s (í + N s ) 
iVt h (l + iVth) 



(4.49) 



and we defined 



f = f(r,N th ,N s ) 



(l + 2N th ) [l + 2 7V th -e- 2r (l + 2jV s )] 
4AT th (l + ^ th ) 



(4.50) 



As one may expect, í s decreases as Nth and N s increase. Moreover, in the limit N s 
reduces to the case of a non squeezed bath, in formula 1 86 , 87 1 



0, the threshold time Í4~4"9"ï 



íd = í s (r,r,JV t h,0) = -ln 



1 



1-e 



2N th 



(4.51) 
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which is always longer than í B . We conclude that coupling a TWB with a squeezed-thermal bath destroys the 
correlations between the two channels faster than the coupling with a non squeezed environment. 

One can also evaluate the threshold time for separability in the case of an out-of-phase squeezed bath, i.e. for 
complex M = \M\ e 10 . The analytical expression is qui te cumbersome and will not be reported here. However, in 
order to investigate the positivity of S as a function of 8, it suffices to consider the characteristic polynomial qs (x) 
associated to S, and study the sign of its roots. This polynomial has four real roots and three of them are always 
positive. The fourth becomes positive adding noise, and the threshold decreases with varying 0. In other words, 
the survival time becomes shorter 1 87 1. 



4.5 Three-mode Gaussian states 

As an example of propagation in a three-mode noisy channel, we investigate now the evolution of state \T) in- 
troduced in Section l2~4l We will refer to the model described in Section 14.1.21 and we will consider three noisy 
channels with the same damping constant Y, same thermal noise A^h and no phase-dependent fluctuation [M/. = 
in Eq. J4.18H . From now on we denote with Nh, h = 1,2,3, the mean photon numbers that characterize the 
state \T) via the Eqs. J2.50I . Accordingly to Eq. (I4.20Í . rearranging the order of the entries, the evolution of the 
covariance matrix Vt is given by the following convex combination of Vt itself and of the stationary covariance 
matrix Voc = (Nth + 5)15 (for the rest of the Section we put k 2 = 2 -1 / 2 ): 

V(t) = e~ n V T + (1 - e- rt ) V m . (4.52) 

Consider for the moment a pure dissipative environment, namely Nth — 0. Applying the separability criterion 
( 13 .3 lt . one can show that the covariance matrix V(t) describes a fully inseparable state for every time t. In fact, 
defining Vk(Í] = V(t) — iJic> with K = A, B,C corresponding to channel (mode) 1, 2 or 3, respectively, we 
have that the minimum eigenvalue of Va (t) is given by 



A^ in = 2e~ Ft 



N! - y/NM + 1) 



(4.53) 



Clearly, \™ m is negative at every time í, implying that mode A is always inseparable from the others. Concerning 
mode B, the characteristic polynomial of Vb(í) factorizes into two cúbic polynomials: 

q 1 (X,T,N 1 ,N 2 ,N 3 ) = -X 3 + 4:[l + e- Tt N 1 \ A 2 

+ 4 [-1 - e- rt (2N 2 + 3N 3 - e~ rt Ni)] X + 8e- rt N 3 (l - e - rt ) , (4.54a) 
q 2 (A, r, iVi, N 2 , N 3 ) = -A 3 + 2 [l + 26-^JVi] A 2 

+ 4 [-e~ rt (2iV2 + N 3 ) + e~ 2rt ^i] A - 8e~ 2r *iV2 . (4.54b) 

While the first polynomial has only positive roots, the second one admits a negative root at every time. Due to the 
symmetry of state \T) the same observation apply to mode C, hence full inseparability follows. Notice that this 
result resembles the case of the TWB state in a two-mode channel studied in the previous Section [see Eq. J4.5 li 
foriV th ^0]. 

When thermal noise is considered (Nth 7^ 0) separability thresholds arise, again resembling the two-mode 
channel case. Concerning mode in channel A, the minimum eigenvalue of matrix Va{í) is negative when 




^(JVl + l)-^! 



(4.55) 



Remarkably, this threshold is the same as the two-mode one given in Eq. ( 14.5 li . if one consider both of them as a 
function of the total mean photon number of the TWB and of state \T) respectively. This consideration confirms 
the robustness of the entanglement of the tripartite state \T) . Concerning mode B, the characteristic polynomial of 
Vs(t) factorizes again into two cúbic polynomials. As above, one of the two have always positive roots, while the 
other one admits a negative root for time t below a certain threshold, in formula: 

- 8e- 2rt 7V 2 + 8(e~ rt - l) e - rí (e- rí ATi - 2N 2 - N 3 )N th 

+ 8(e~ rt - 1) 2 (1 + 2e- rí A r 1 )A7 t 2 h - 8(e~ rt - l) 3 7V t 3 h < . (4.56) 



Mode C is thus subjected to an identical separability threshold, upon the replacement N 2 «-> A^. 



Chapter 5 



Quantum measurements on continuous 
variable systems 



In this Chapter we describe some relevant measurements that can be performed on continuous variable (CV) sys- 
tems. These include both single-mode and two-mode (entangled) measurements. As single-mode measurements 
we will consider direct detection of quanta through counters or on/off detectors and homodyne detection for the 
measurement of the field quadratures. As concerns two-mode entangled measurements, we will analyze the joint 
measurement of the real and the imaginary part of the normal operators Z± — aztb^, a and b being two modes of 
the field, through double homodyne, six-port homodyne or heterodyne-like detectors. Throughout the Chapter we 
will mostly refer to implementation obtained for the radiation field. This is in order since it is in this context that 
they have been firstly developed, and are available with current technology. However, it should be mentioned that 
the schemes analyzed in this Chapter can also be realized, or approximated, also for other fields, as for example 
in atòmic or condensate systems. The measurement schemes will be described in some details in order to evaluate 
their positive operator-valued measure (POVM), as well as the corresponding characteristic and Wigner functions, 
both in ideal conditions and in the presence of noise, i.e. of non-unit quantum efficiency of the detectors. In Sec- 
tion l5.1l we review the concept of POVM and its relations with customary measurement of observables, whereas 
in Section l5~2l we briefly review the concept of moment generating function. Direct detection of the field, either by 
counting or by on/off detectors, is the subject of Section l5~3l while homodyne detection of the field quadratures is 
analyzed in Section l5~5l Finally, the joint measurement of 5Re[Z] and 9m[Z] is analyzed in Section l5"6l 



5.1 Observables and POVM 

In order to gain information about a quantum state one has to measure some observable. The measurement process 
unavoidably involves some kind of interaction, which couples the mode under examination (the signal) to one or 
more other modes of the field (the probe). Therefore, one has to admit that, in general, the measured observable 
is not defined on the sole Hilbert space of the signal mode. Rather, it reflects properties of the global state which 
results from the interaction among the signal mode and the set of the probe modes. In some cases, it is possible 
to get rid of the probe modes, such that the statistics of the outcomes can be described in terms of an observable 
defined only on the Hilbert space of the signal mode. As we will see, this is the case of homodyne detection of 
a field quadrature. More generally, eliminating the probe modes by partial trace, we are left with a more general 
object, that is a spectral measure of an observable to describe the statistics in terms of the signal's density matrix. 
Let us denote by Tí the Hilbert space of the signal, by K. the Hilbert space of the probe modes, and by X the 
measured observables on Tí (g) JC. The spectral measure of X is given by x — > dE(x) — \x) (x\dx with x E X C K 
(the spectrum of X) and (x\x') — S(x — x') 1 . The probability density of the outcomes is thus given by 

p(x) = Tr[g® aE(x)} , (5.1) 

where g and a are the initial preparations of the signal and the probe respectively and the trace is taken over all the 
Hilbert spaces. Eq. ( 15. U can be written as 

p(x) = Tr H [g Trjç [a E(x)] ] = Tr„ [g U(x)} , (5.2) 

where H(x) = Tr^[cr E(x)] is usually referred to as the positive operator-valued measure (POVM) of the mea- 
surement scheme 2 . From the definition we have that a POVM is a set {n(x)} a , gAr of positive, H(x) > (hence 

'For observables with a discrete spectrum the spectral measure reads k —* ïl k = \k) (fc| with (k\l) = <5 fei 
2 POVM are also sometimes referred to as probability operator measure (POM) 
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selfadjoint), and normalized J x dx Tl(x) — I operators while, in general, they do not form a set of orthogonal 
projectors. 

In summary, a detection process generally corresponds to the measurement of an observable defined in the global 
Hilbert space of the signal and the probe. If we restrict our attention to the system only, the statistics of the out- 
comes is described by a POVM. The converse is also true, i.e. whenever a set of operators satisfying the axioms 
for a POVM is found, then the following theorem assures that it can be seen as the measurement of an observable 
on a larger Hilbert space H88II89I . 

Theorem (Naimark) If{Tl(x)} x( - x is a POVM on TL then there exist a Hilbert space IC, a spectral measure E(x) 
on Tí ® JC and a density operator a on IC such that IL(x) = Trjç [a E(x)]. 

In addition, the number of these Naimark extensions is infinite, corresponding to the fact that a given POVM may 
results from different physical implementations. We will see an example of this property in Section l5"6l 



5.2 Moment generating function 



For a detection scheme measuring the observable X, with eigenvalues x £ X C K, the so-called moment generat- 
ing function (MGF) is defined as 



M x {y)=Tx[Ré vX ] , 



(5.3) 



where R is the overall quantum state (signal plus probe) at the detector. MGF generates the moments of the 
measured quantity X according to the formula 



<*"> = H)"^ M x (y)\ y=0 . 
The MGF Mx(y) also provide the distribution of the outcomes p(x) through its Fourier transform. In fact, 



(5.4) 



Z7T 



"•H j^x-x) 



= Tr [R S(X - x)} = Tr [R \x) (x\] = p(x) . 



(5.5) 



The density matrix R in Eq. i5.5\ should be meant as the global quantum state, system plus probe, at the input of 
the detector. In turn, the trace should be performed over the global Hilbert space Ti®JC describing all the degrees 
of freedom of the detector. Comparing Eq. (15. 5> with Eq. d5.2t we have that for R = g ® a the POVM of the 
detector, obtained by tracing out the probes, can be expressed as follows 



Tl{x) = Tifc 



2tt 



(5.6) 



Eq. d5.5l > can be generalized to the multidimensional case. In particular, any detector measuring a couple of 
commuting operators [X, Y] = can be seen as measuring the complex normal operator Z = X + ïY . The MGF 
is defined as 



M Z (X) = Tr 



Re 



XZ^-X'Z 



(5.7) 



with À 6 
transform 



whereas the probability distribution of the outcomes a 6 C is obtained as the complex Fourier 

d 2 \ 



p(a)= I ^e x ' a - a * x M z (X) 



Again, in Eq. i5.7\ . R denotes the overall quantum state at the detector; the POVM can be evaluated as 

n(a) = Tr/c 



d X c a(ZÏ-X')-a*(Z-X) 



(5.8) 



(5.9) 



5.3 Direct detection 

By direct detection we mean the measurement of the quanta of the field, either by effective counting (i.e. discrim- 
inating among the number of incoming quanta) or just by revealing their presence or absence (on/off detection). 
In the following we analyze in some details the detection of photons. Analogue schemes have been developed for 
atòmic systems. 



5.3 Direct detection 
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5.3.1 Photocounting 

Light is revealed by exploiting its interaction with atoms/molecules or electrons in a sòlid: each photon ionizes a 
single atom or promotes an electron to a conduction band, and the resulting charge is then amplified to produce 
a measurable pulse. In practice, however, available photodetectors are not ideally counting all photons, and their 
performances are limited by a non-unit quantum efficiency Ç, namely only a fraction £ of the incoming photons 
lead to an elèctric signal, and ultimately to a count: some photons are either reflected from the surface of the 
detector, or are absorbed without being transformed into elèctric pulses. 

Let us consider a light beam entering a photodetector of quantum efficiency £, i.e. a detector that transforms 
just a fraction £ of the incoming light pulse into elèctric signal. If the detector is small with respect to the coherence 
length of radiation and its window is open for a time interval T, then the probability p(m; T) of observing m counts 
is a Poissonian distribution of the form 1 90 1 



p m (T) =Tr 



(5.10) 



where g is the quantum state of light, : ■ • ■ : denotes the normal ordering of field operators, and I{T) is the beam 
intensity 

I(T) = ^ £ dtE^\t) ■ E^ + \t) , (5.11) 

given in terms of the positive, E^ + \ and negative, E^~\ frequency part of the elèctric field operator. The quantity 
p(T) = ÇTt[qI(T)] equals the probability of a single count during the time interval (T,T + dt). Let us now 
focus our attention to the case of the radiation field excited in a stationary state of a single mode at frequency u>. 
Eq. d5 . 1 Oi can be rewritten as 



Prn(v) = Tr 



(770*0)*" , f , ' 
g : : — exp |— rja' o) : 



(5.12) 



where the parameter r) denotes the overall quantum efficiency of the photodetector. By means of the identities 

: {a)a) m : = (aï) m a m = 0*0(0*0 - 1) . . . (a*a - m + 1) and :e -3!ato : = (1 - x) aia |91 1, one obtains 

Pm(v) = £ Skk QV TO (1 " V) k - m , (5-13) 



where g^k = (k\g\k) — Pk(i] = 1)- Hence, for unit quantum efficiency, a photodetector measures the photon 
number distribution of the state, whereas for non-unit quantum efficiency the output distribution of counts is given 
by a Bernoulli convolution of the ideal distribution. Eq. ( 15.131 1 can be written as p m {rj) = Tx\g Il m (r7)] where the 
POVM of the photocounter is given by 



k= 



Notice that Ii m (rj) > and H m (rj) = ï, but [H m (rj), üfc C 7 ?)] 7^ 0, i.e. they do not form a set of orthogonal 
projectors. The corresponding characteristic and Wigner functions can be easily obtained from that of a number 
state |fc)(A;|, namely 

X [\k)(k\](X) = (k\D(X)\k) = e- 1 ^ 2 L k (\X\ 2 ) , (5.15a) 
W[\k)(k\](a) = -(k\(-f a D(2a)\k) = -{-) k e' 2 ^ L fc (4|a| 2 ) , (5.15b) 

where Lk{x) is a Laguerre polynomials. We have 

x[n,„(,)](A) = ÍL m ('^')expj-í 5 -2|A| ï j , (5.16a) 

w-m») - Ij^L· u («!) =xp{-^l m»} . (5..*) 

The effects of non-unit quantum efficiency on the statistics of a photodetector, namely Eqs. J5.131 and (15. 14-1 . 
can be also described by means of a simple model in which the realistic (not fully efficient) photodetector is 
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replaced with an ideal photodetector preceded by a beam splitter of transmissivity cos 2 </>, with the second mode 
left in the vacuum state. The reflected mode is absorbed, whereas the transmitted mode is photodetected with unit 
quantum efficiency. The probability of measuring m clicks in such a configuration is given by 



p m {4>) = Tr a 



UéQ< 



U\ \m){m\ (g) I 



(5.17) 



where we denoted by a and b the two involved modes, U<f, is the unitary evolution of the beam-splitter (see Section 
ll.4.2> and g the initial preparation of the signal. Using the cyclic properties of the (full) trace and, then, performing 
the partial trace over the vacuum mode, we have 



Pm{4>) 



Tr ab 
Tr„ 



g® |0)(0| t/J|m)<m|®II7^ 
g (0| U\ \m) (m\ ® I |0> 1 = Tr Q [g TL„ 



(5.18) 



Eq. J5.17> reproduces the probability distribution of Eq. ( I5.13> with -q = cos 2 cf>. We conclude that a photodetector 
of quantum efficiency rj is equivalent to an ideal photodetector preceded by a beam splitter of transmissivity 77 
which accounts for the overall losses of the detection process. 

If we have more than one mode impinging on a photocounter, we should take into account that each click may 
be due to a photon coming from each of the mode. The resulting POVM assumes the form 



n m(T7) = ■- ® -■ ® n k n (Vn) 6 í ^Tk s - ml , 

fc 1= fc„=u \s=l / 

where we have also supposed that each mode may be detected with a different quantum efficiency. 



(5.19) 



5.3.2 On/off photodetectors 



As mentioned above, in a photodetector each photon ionizes 
a single atom and, at least in principle, the resulting charge 
is amplified to produce a measurable pulse. Taking into ac- 
count the quantum efficiency, we conclude that the resulting 
current is proportional to the incoming photon flux and thus 
we have a linear detector. On the other hand, detectors op- 
erating at very low intensities resort to avalanche process in 
order to transform a single ionization event into a recordable 
pulse. This implies that one cannot discriminate between a 
single photon or many photons as the outcomes from such 
detectors are either a click, corresponding to any number of 



Q 



V 



Figure 5.1: Model of a realistic on/off photodetector 
with non-unit quantum efficiency r\, and non-zero dark 
counts. 

photons, or nothing which means that no photons have been revealed. These Geiger-like detectors are often re- 
ferred to as on/off detectors. For unit quantum efficiency, the action of an on/off detector is described by the 
two-value POVM {Ilo = |0) (0| , lli = I — üo}, which represents a partition of the Hilbert space of the signal. In 
the realistic case, when an incoming photon is not detected with unit probability, the POVM is given by 



Hofa) = X> " \k){k\ , IMrç) = I - 110(77) , 



(5.20) 



k=0 



with 77 denoting quantum efficiency. The corresponding characteristic and the Wigner functions can be easily 
obtained from that of a number state [see Eqs. ( 15. 15al > and <5.15bH . We have 



Xpo(»7)](A) - 
W[n (r7)](a) = 



exp 
2 



2,, 



tt2- 



V 



■ exp 



|A| 2 
2 '/ 



" X [n 1 (r/)](A)=7r^ 2 )(A)-x[no(r ? )](A), 
W[n 1 (r,)](a) = --W[U (r 1 )}(a) 



2-7? 



(5.21a) 
(5.21b) 



Besides quantum efficiency, i.e. lost photons, the performance of a realistic photodetector are also degraded by 
the presence of dark-count, namely by "clicks" that do not correspond to any incoming photon. In order to take into 
account both these effects we use the simple scheme introduced in the previous Section and depicted in Fig. 15. II 
A real photodetector is modeled as an ideal photodetector (unit quantum efficiency, no dark-count) preceded by 
a beam splitter of transmissivity equal to the quantum efficiency 77, whose second port is in an auxiliary excited 
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state v, which can be a thermal state, or a phase-averaged coherent state, depending on the kind of background 
noise (thermal or Poissonian) we would like to describe. When the second port of the beam splitter is the vacuum 
v = |0) (0|, we have no dark-counts and the POVM of the photodetector reduces to that of Eq. J5.20b . On the other 
hand, when the second port of the BS excited in a genèric mixture 

v = ^v ss \s)(s\ , 

S 

then the overall POVM describing the on/off photodetection is expressed as the following generalized convolution 

OC OO / i \ 

nS(rO=TtJc^ e ®i/tf}l®|0)(0|l = J> -»?)"$>» ( ) \n)(n\ , (5.22) 



rtir^ ln\/nll _ Wi _J»V„ ^í n + s 

s 

whereas the characteristic and the Wigner functions read as follows 

|2 



xlKivm = -exp(-^|A| 2 )5> s L s f (1 * )|A| ) , (5.23a) 



The density matrices of a thermal state and a phase-averaged coherent state (with mean photons) are given by 

1°°/ \ s OO / w 

„ 1 ^1 Ub \ l„\/„l ,. _„-n^ W 



4.TEUrXT l-X-l . v P = e-"»J2^\s){s\. (5.24) 
n b + í^\n b + lj n s! 

s—O v 7 s—O 

In order to reproduce a background noise with mean photon number TV we consider the state v with average 
photon number rij, = N / (1 — rj). In this case we have 

1 x / \ n 

n—O ^ 



n p ( v ,N) = e- N Y,(l-v) n L n {-^) \n){n\, (5.25b) 



where i and p denotes thermal and Poissonian respectively. The corresponding characteristic and Wigner function 
are given by 

ït/„ s mm_ 1 / 2(l + iV) -77,^,2 



x[IÍ (r?,iV)](A) = -exp|-^ ^ í |A| 2 }> , (5.26a) 



x [ng(r ? ,Ar)](A) = Íexp|-^|A| 2 }j Í2J-^|Apj (5.26b) 



1 2 2ï ? ,_.,a 



and 

W ^(V,N)](a) = - 2{1 + - N) _ v ^{- 2{i ;^_ v \«q ■ (5.27a) 
I^[n^iV)](a) = i^-cxp|-^-(7V+H 2 )|/ (^^) , (5.27b) 
respectively, where lo (x) is the 0-th modified Bessel function of the first kind. 

5.4 Application: de-Gaussification by vacuum removal 

As we have already pointed out, Gaussian states are very important for continuous variable quantum information. 
However, there are situations, as for example in testing nonlocality with feasible measurements (see Chapter|6jl, 
where one needs to go beyond Gaussian states. Indeed, when the Gaussian character is lost, then immediately 
the Wigner function of the state becomes negative, for pure states, hence stronger nonclassical properties should 
emerge. An effective method to "de-Gaussify" a state is through a conditional measurement, and, in particular, 
by elimination of its vacuum component leading to a state which is necessarily described by a negative Wigner 
function. In the next two Sections this strategy will be applied both to the TWB and the tripartite state |T) given 
in Eq. ( I2.49> through the on/off detection scheme introduced in the previous Section. 
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5.4.1 De-Gaussification of TWB: the IPS map 

In this Section we address a de-Gaussification process onto the twin-beam state (TWB) of two modes of radiation 
|A))o6 = \/l — \ 2 J2^=o A n |n, n )ab, where we assume the TWB parameter A = tanhr as real, r being referred 
to as the squeezing parameter. The corresponding Wigner function is given by 



W r (a, 0) = exp{-2A(|a| 2 + |/3| 2 ) + 2B(af3 + a*(3*)} . 



(5.28) 



with A = A(r) = cosh(2r) and B = B(r) = sinh(2r). 

The de-Gaussification of a TWB can be achieved by subtracting photons from both modes through on/off 
detection 1 92 93 , 94 1 . Since the scheme does not discriminate the number of subtracted photons, we will refer to 
this process as to inconclusive photon subtraction (IPS). 

The IPS scheme is sketched in Fig. 15.21 The modes a 
and b of the TWB are mixed with vacuum modes at two 
unbalanced beam splitters (BS) with equal transmissivity 
t = cos 2 <j>\ the refiected modes c and d are then revealed 
by avalanche photodetectors (APD) with equal efficiency 
-q. APD's can only discriminate the presence of radia- 
tion from the vacuum. The positive operator-valued mea- 
sure (POVM) {n (?7), n 1 (?7)} of each detector is given 
in Eq. (15.201 . Overall, the conditional measurement on 
modes c and d, is described by the POVM (assuming 
equal quantum efficiency for the photodetectors) 




Figure 5.2: Scheme of the IPS process. 



nbi(»í) = Do,c(»?)<8ni i(í (»/), 

n n (í/) =n hc (r 1 )<g)U hd {r 1 ) . 



n 00 (?7) = n 0!C (?7) <g> n 0)(í (r?) , 

n 10 (r?) = n ljC (r?) ® Tlo,d(v) , 
When the two photodetectors jointly click, the conditioned output state of modes a and b is given by |92 9 

Tr cd [£/ gc (0) g U bd {4>) R g |Q) CC (Q| g |0) drf (0| U^W g U}j$) l a ®l b ® llufa)] 



(5.29a) 
(5.29b) 



£{R) = 
where U ac (<j>) = exp{- 



Pu(r, <t>,ri) 



(5.30) 



iJr 3 — nr^ 



)} and Ubd{<t>) are the evolution operators of the beam splitters and R the 
density operator of the two-mode state entering the beam splitters (in our case R = p TWB = |A)) a b(, a ((A|). The 
partial trace on modes c and d can be explicitly evaluated, thus arriving at the following decomposition of the IPS 
map 3 . We have 



£{R) 



OC 



m p (<t>, 77) M pq {4>) R Aft (0) m g ((f>, 77) 



where 



= t^cf > {l (1 T,)"] = aPbq{cQS(j) fa +b H_ 



(5.31) 



(5.32) 



Now we explicitly calculate the Wigner function of the state g IPS = £(g vm ), which, as one may expect, is no longer 
Gaussian and positive-definite. The state entering the two beam splitters is described by the Wigner function 



Wr(a,(3,C,0 =W r (a,f3)— ex P {-2|(| 2 - 2|£| 2 } , 



(5.33) 



where the second factor at the right hand side represents the two vacuum states of modes c and d. The action of 
the beam splitters on W^"' can be summarized by the following change of variables (see Section [T.4.2> 



a — > ot cos <fi + Ç sin 4> , Ç - 
(3 — ► /3 cos cj) + £ sin <f> , £ - 

and the output state, after the beam splitters, is then given by 



• £ cos <j) — a sin <f> , 
£ cos cj) — /3 sin cj) , 



(5.34a) 
(5.34b) 



W%(a, (3, C, = ^ W r ,+ (a, (3) exp {-a|£| 2 +w£ + w*C } 

x exp { - a|C| 2 + (v + 2BÍ, sin 2 0)C + (v* + 2BÇ* sin 2 0)C* } , (5.35) 

should be intended as a polyindex so that 15.311 reads 



3 Eq. 15.311 is indeed an operator-sum representation of the IPS map: {p, q] 

S{R) ='£ e AeRAlmÜiAe = [pii{r,<f>,ri)]~ 1/2 rn p {4>,v) M pq {4>). 



5.4 Application: de-Gaussification by vacuum removal 



47 



where 

W r ^(a,(3) = 4 exp{-b(\a\ 2 + \(3\ 2 ) + 2Bcos 2 (j)(a(3 + a*l3*)} (5.36) 



and 



a = a(r, <j>) = 2(Asin 2 <j> + cos 2 0), (5.37a) 

b = b(r, 4>)=2(A cos 2 <j) + sin 2 <jy) , (5.37b) 

v = v(r,4>) = 2cos^ sin0[(l - (5.37c) 

to = ic(r, 0) = 2 cos sin 4> [(1 - A)/3* + Ba] . (5.37d) 



At this stage on/off detection is performed on modes c and d (see Fig. 15. 2> . We are interested in the situa- 
tion when both the detectors click. The Wigner function of the double click element 1111(77) of the POVM [see 
Eq. d5"29l l is given bv I92ll96l 

w v (C,0 = wpiiC^íCO = 4{i - Qv(0 - Qv(0 + Qr,(C)Qv(0] , (5-38) 



with 



) t ,(z) = -^-cxpJ --^-|z| 2 . (5.39) 
2-?y 2-?/ 



Using Eq. ( I5.30> and the phase-space expression of trace for each mode [see Eq. dl . 1011 1. the Wigner function of 
the output state, conditioned to the double click event, reads 

W r ^ n (a, (3) = , (5.40) 

where f(a,(3) = f r ,<j,,r,(a, (3) with 

f(a, (3)=tt 2 f d 2 C d 2 Ç ^ WrA*, P) E ^T 1 G rl „(«> & C, . (5-41) 

and pn(r,(j),rf) is the double-click probability reported above, which can be written as function of f(a,(3) as 
follows 

p n {r,(t>,r]) =n 2 í d 2 a d 2 (3 f(a, (3) . (5.42) 
Jc 2 



The quantities G* l (a, (3, Ç, £) in Eq. ( I5.41> are given by 



G l lja : (3, C, = exp { - x k \(\ 2 + (v + 2BÇ sin 2 <j>)C + (v* + 2BÇ sin 2 <j>)Ç } 

xexp{-7/ fe |e| 2 +< + W *r} , (5.43) 



where the expressions of x k = Xk(r, 4>, 77), y k = yk{ r i cj>, 77), and Ck(rj) are reported in Table l5.ll 

The mixing with the vacuum in a beam splitter with transmissivity r followed by on/off detection with quantum 
efficiency r\ is equivalent to mixing with an effective transmissivity [92 1 

r cff ee r cff (0, 77) = 1 - ïy(l - r) (5.44) 

followed by an ideal (z'.e. efficiency equal to 1) on/off detection. Therefore, the state ( 15.401 can be studied for 
77 = 1 and replacing r = cos 2 0=1 — sin 2 with r e ff. Thanks to this substitution, after the integrations we have 



16C fc 



7T 2 ïfcT/fc - 452(1 - T Gff ) 2 

x exp{(/ fc - Ò)H 2 + ( 5fc - 6)|/3| 2 + {2Bt cS + h k )(a(3 + a*f3*)} (5-45) 

and 

/ v ^ lQ[x k y k -AB 2 {l~T cíi ) 2 ]- l C k 

Pii{r, T eS ) = > — — r — , (5.46) 

l~{( b ~ fk)( b ~ 9k) - (2Bt cS + h k y 
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Table 5.1: Expressions of Ck, Xk, and yk appearing in Eq. 15.431 



where we defined C k = C k (í) and 

f k = f k (r, T cff ) = N k [x k B 2 + AB 2 (l - A)(l - r off ) +y k {l- Af] 
g k = g k (r, r cff ) = N k [x k (l - Af + AB 2 {\ - A)(l - r off ) + y k B 2 } 
h k = h k {r,T cS ) = N k {{x k + y k )B(l - A) + 2B(B 2 + (1 - A) 2 )(í 

Nk = Nk{r ^ s) = XkVk -^ { i-r cS y 

In this way, the Wigner function of the IPS state can be rewritten as 

4 



1 



7T 2 Pn(r, r cff ) 



J2C k W k (a,(3), 



k=l 



where we introduced 



C fe = C k (r, t oS ) 



4C k 



eff)] 



x k y k - 4B 2 (1 - r cff ) 2 ' 
and defined 

W k (a,0) = exp{(/ fc - 6)|a| 2 + (g k ~ ò)|/3| 2 + {2Bt cS + h k )(a/3 + a*(3*)} . 
Finally, the density matrix corresponding to W ws (a 1 j3) reads as follows |92 1 

1 - A 2 



Pu(r,T eS ) 

in[ri 

E 



n,m— 
Minin, rn 



h.k=0 



l-r eff \ fe+fc 

Teff / 



(5.47a) 
(5.47b) 
(5.47c) 

(5.47d) 



(5.48) 



(5.49) 



(5.50) 



x \n - k) a \n - h) bb (m - h\ a (m - k\ , (5.51) 



with A = tanh r. 

The state given in Eq. J5.48Í is no longer a Gaussian state. Its use in the enhancement of the nonlocality 
1971 198 95 1 and in the improvement of CV teleportation |92 1 will be investigated in Chapter|6]and0 respectively. 



5.4.2 De-Gaussification of tripartite state: the TWBA state 

In this Section we consider the tripartite state |T) given in Eq. ( I2.49Í as a source of two-mode states. In particular, 
we analyze two-mode non-Gaussian state obtained by a conditional measurement performed on it. Due to the 
structure of the state |T), its vacuum component can be subtracted by a conditional measurement on mode a^, the 
same observation being vàlid for mode a^- Let us consider on/off detection performed on mode a^. The three-mode 
two-valued POVM is {üg 3 ^ (77), IT^ (??)}, with the element associated to the "no photons" result given by 



(l-^)»33<n| 



(5.52) 
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The probability of a "click" is 

P 1 = Pi(N 3 ,ri) = Tr 123 [\T)(T\ l·lf^r?) 
while the conditional output state reads as follows 

)t){t\u[ 3 \ v ) 

1 + vN 3 



V N 3 



(1 + r)N 3 ) ' 



(5.53) 



1 T 

— Tr 3 

Pi 



{1 + N± + N 2 ) V N 3 ^ 



Eí— 



p 1 - íl - n) p + 

— ^ r ^(4riA)} 1221 «A|< 



(5.54) 



where we denote by | A)) 12 the TWB state of the modes a\ and a 2 with parameter A = J N2 / ( 1 + ÏVx ) (see 
Section lT.4.4> . We indicated this state with a subscript TWBA (i.e. TWB added) since it corresponds to a mixture 
of TWBs with additional photons in one of the modes. In order to evaluate its Wigner function we use J5 . 2 1 al > for 
the characteristic function of H\ (77), hence the characteristic function of g™ BA is given by 



X[íw] (Ai, A 2 ) = — i X [\T) (T\}(Xx, A 2 , 0)- 



V JC T 



X [|T)(T|](A 1 ,A 2 , M )exp 



2 V 



- N 



(5.55) 



After some àlgebra the Wigner function associated with state g TV/BA can now be calculated. It reads as follows 



^(^ = 1+^ (ï) 2 í -=À= exp Í-Y* (V> T y X y] 
4riN 3 \irj [ ^Det [V T ] l V ; - 1 

-- ; 2 expí-y^p- 1 )^}! , (5.56) 

where Y = (xi, x 2 , yi, ï/ 2 ) t , and D = V T + Diag(0,0, ^,0,0, 2=2), F T being defined in Eq. (ESB- In order 
to simplify the notation we have indicated with O' the 4x4 matrix obtained from the 6x6 matrix O deleting the 
elements corresponding to the third mode (3-rd row/column and 6-th row/column), due to the trace over the 3-rd 
mode. Nonlocality properties of the TWBA state will investigated in Chapter[6] 



5.5 Homodyne detection 

Homodyne detection schemes are devised to provide the measurement of a single-mode quadrature x$ through the 
mixing of the signal under investigation with a highly excited classical field at the same frequency, referred to as 
the local oscillator (LO). Homodyne detection was proposed for the radiation field in Ref. 1991 . and subsequently 
demonstrated in Ref. 1 1 001 . For the radiation field quadrature measurements can be achieved by balanced and 
unbalanced homodyne schemes, whereas realizations for atòmic systems have also been proposed 1 101 1. 

5.5.1 Balanced homodyne detection 

The schematic diagram of a balanced homodyne detector is reported in Fig. 15.31 The signal mode a interferes with 
a second mode b excited in a coherent semiclassical state (e.g. a làser beam) in a balanced (50/50) beam splitter 
(BS). The mode b is the LO mode of the detector. It operates at the same frequency of a, and is excited in a coherent 
state \z) with large amplitude z. The BS is tuned to have real coupling, hence no additional phase-shift is imposed 
on the reflected and transmitted beams. Moreover, since in all experiments that use homodyne detectors the signal 
and the LO beams are generated by a common source, we assume that they have a fixed phase relation. In this case 
the LO phase provides a reference for the quadrature measurement, namely we identify the phase of the LO with 
the phase difference between the two modes. As we will see, by tuning <fi = arg[z] we can measure the quadrature 
x^ at different phases <\>. After the BS the two modes are detected by two identical photodetectors (usually linear 
photodiodes), and finally the difference of photocurrents at zero frequency is electronically processed and rescaled 
by 2\z\. According to Eqs. dl.63> and dl.70> . denoting by c and d the output mode from the beam splitter, the 
resulting homodyne photocurrent H is given by 
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Figure 5.3: Schematic diagram of the balanced homodyne detector. 



Notice that the spectrum of the operators a% + tfa is discrete and coincides with the set Z of relative integers. 
Therefore the spectrum of the homodyne photocurrent H is discrete too, approaching the real axis in the limit of 
highly excited LO (\z\ 3> 1). We now exploit the assumption of a LO excited in a strong semiclassical state, i.e. 
we neglect fiuctuations of the LO and make the substitutions ò — > z, b> — > z* . The moments of the homodyne 
photocurrent are then given by 



H 



H 2 =xl 



4ÜP 



xl n ~ 2 f x 2 , 



a) a 
4ÜP 



(5.58) 



which coincide with the quadrature moments for signals satisfying (a^a) <C 4|z| 2 . In this limit the distribution of 
the outcomes h of the homodyne photocurrent is equal to that of the corresponding field quadratures. The POVM 
{fi/,} of the detector coincides with the spectral measure of the quadratures 



n, 



* II(x) = = SÇxj, - x) , 



(5.59) 



i.e. the projector on the eigenstate of the quadrature x$ with eigenvalue x. In conclusion, the balanced homodyne 
detector achieves the ideal measurement of the quadrature x^ in the strong LO regime. In this limit, which sum- 
marizes the two condi tions i) \z\ ^> 1 to have a continuous spectrum and ii) \z\ 2 ^> (ala) to neglect extra terms 
in the photocurrent moments, the probability distribution of the output photocurrent H approaches the probability 
distribution p(x, <j>) = <t > {x\g\x) l p of the quadrature x^ for of the signal mode a, The same result 1 102 1 can be 
obtained by evaluating the moment generating function Aí#(/x) = Tr [g <g) \z)(z\ e lfJ,H ] . Using the disentangling 
formula for SU(2) dl .60l > we have 



M H (p) 



i tanf -^r-T )b* a 

■d \2\z \ J 



JJ_ 

2\z\ 



o ítan( 5 ^ T )atb 



(5.60) 



Since mode b is in a coherent state | z) the partial trace over b can be evaluated as follows 



M h (jí) 



cos 



JL 

2\z\ 







-6 f 6 




\ ( » V 






cos — — - 






. \2\z\J. 





Now, rewriting ( 15. 611 in normal order with respect to mode a we have 

Mk( M ) = ^e lzsin (*) at exp|-2sin 2 (afa + \z\ 2 ) \ < ' 

In the strong LO limit J5.62Í becomes 



lim M H (p) = { e^ 8 '*"' exp 



= (cxp{i l· ix c f,}) a . 



The generating function in (I5.63Í then corresponds to the POVM 



dfi. 



= / 7T~ exp{i^j,(xé - x)} = 8{xé - x) = \x)éé{x\ , 

Z7T 



(5.61) 



(5.62) 



(5.63) 



(5.64) 
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which confirm the conclusions drawn in Eq. ( I5.59> . 

In order to take into account non-unit quantum efficiency at detectors we employ the model introduced in 
the previous Sections, i.e. each inefficient detector is viewed as an ideal detector preceded by a beam splitter 
of transmissivity 77 with the second port left in the vacuum. The homodyne photocurrent is again formed as the 
difference photocurrent, now rescaled by 2\z\q. We have 



1 

2\z~\ 



2r] 



(u + v) 



ò f + h. 



(5.65) 



where only terms containing the strong LO mode 6 are retained, and u and v denote the additional vacuum modes 
introduced to describe loss of photons. The POVM is obtained by replacing 



1-77 
2-q 



(U0 + V,),) 



(5.66) 



in Eqs. ( 15. 641 . with — ^(w^e 1 ^ + we w — u,v, and tracing the vacuum modes u and v. One then obtains 



n v (x) 



2nS 2 



exp 



(X0 - xf 



2^ 



dy exp 



251 



\v)4>4>(v\ 



where 



3 = 



1 — 7] 
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(5.67) 



(5.68) 



Thus the POVMs, and in turn the probability distribution of the output photocurrent, are just the Gaussian convo- 
lutions of the ideal ones. 

The Wigner functions of the homodyne POVM is given by 



W[U v (x)](a) 



1 



2^1 



: exp 



2^ 



(5.69) 



which leads to l¥[II(x)](c 



7J— 1 



8(x — |(ae ^ + a*e 1 ^)) in the limit of an ideal homodyne detector. 



5.5.2 Unbalanced homodyne detection 



The scheme of Fig. l5.4l is known as unbalanced homodyne detec- 
tor and represents an alternative method to measure the statistics 
of a field quadrature. The signal under investigation is mixed with 
the LO at a beam splitter with transmissivity r = cos 2 cf). The re- 
flected beam is then absorbed, whereas the transmitted beam is 
revealed through a linear photocounter. If a is the signal mode 
and b the LO mode the transmitted mode c can be written as 
c = a cos <j) + b sin <j). The detected photocurrent is given by 





1 T 


9 v - 


LO 



Figure 5.4: Schematic diagram of unbalanced ho- 
modyne detector. 



era cos 2 



sin 2 (j) + (a^ b + b' a) sin <p cos <j) 



(5.70) 



and the unbalanced homodyne photocurrent is obtained by a simple rescaling 



2skr 



where z is again the LO amplitude. Upon tracing over the local oscillator, in the limit of \z\ (aJa), we have for 
the first two moments 



I H ) = l 



1 



\z \ tant 
1 

\z \ tan< 



-(x e )+0(\z\- 2 ), 



1 



(x e 



UI 2 tan 2 , 



■0(\z\ 



(5.71) 
(5.72) 



and therefore (A/^) = (\z\ 2 tan 2 (f)^ 1 (/S.X 2 ,), where 9 is the shift between signal and LO. The procedure can be 
generalized to higher moments, thus concluding that through unbalanced homodyne one can recover the statistics 
of the field quadratures. In order to minimize the effect of LO, the regime </> <C 1 with \ z\(j> finite should be adopted. 



52 



Chapter 5: Quantum measurements on CV systems 



5.5.3 Quantum homodyne tomography 

The measurement of the field quadrature x<p for all vàlues of the phase </> provides the complete knowledge of the 
state under investigation, i.e. the expectation vàlues of any quantity of interest (including quantities not directly 
observable). This kind of measurement is usually referred to as quantum homodyne tomography [102 103 1 for 
reasons that will be explained at the end of this Section. 

In order to see how the knowledge of p(x, 4>) — ^(x\g\x)^ allows the reconstruction of any expectation value 
let us rewrite the Glauber formula dl.37> changing to polar variables A = (~i/2)ke 1 ^ 

0= T— f — Tr [0 e* kx *} e~ ikx * , (5.73) 

JO T JR 4 

which shows explicitly the dependence on the quadratures x$. Taking the ensemble average of both members and 
evaluating the trace over the set of eigenvectors of x$, one obtains 

(0}= í ^ / dxp(x,<p)TZ[0}(x,<p), (5.74) 
Jo k JR 

The function lZ[0](x,cf>) is known as kernel or pattern function for the operator O, its trace form is given by 
TZ[0](x, <j>) = Tï[OK{x ( j } — x)] where K(x) writes as 

f dk 1 

K{x) = / — \k\e lkx = -3te / dk ke lkx . (5.75) 
Jr 4 2 Jo 

Therefore, upon calculating the corresponding pattern function, any expectation value can be evaluated as an 
average over homodyne data. Remarkably, tomographic reconstruction is possible also taking into account nonunit 
quantum efficiency of homodyne detectors, i. e. upon replacing p(x,(f>) with p v (x,cf>). Indeed, we have 

(Q)= [ ^ í dx Pn (x, 0) K n [0]{x, 0) , (5.76) 



where the pattern function is now lZ n [O] (x,(f) = Tr[(9 K v (x^ — x)), with 

K v (x) = iaíe J dk k exp ji^fc 2 + ikx^ . (5.77) 

Notice that the anti-Gaussian in J5.77Í causes a slower convergence of the integral J5 . 76b and thus, in order to 
achieve good reconstructions with non-ideal detectors, one has to collect a larger number of homodyne data. As 
an example, the kernel functions for the normally ordered products of mode operators are given by 1 104| |105I 

■RJa^a m ]( X , <j>) = e Krn-n)^Hn+ J^x) 

V(27?)"+-("+ m ) 

where H n (x) is the n-th Hermite polynomials, whereas the reconstruction of the elements of the density matrix in 
the number representation g nm = Ti[g \n)(m\] corresponds to averaging the kernel 



n v [\n)(n + d\]{x, ( j))=e^ + V^j—^ l J^dk\k\e^ ík ~ l2kx k d L d n {k 2 ) , (5.79) 

where L„ (x) denotes the generalized Laguerre polynomials. Notice that the estimator is bounded only for rj > 1/2, 
and below the method would give unbounded statistical errors. 

The name quantum tomography comes from the first proposal of using homodyne data for state reconstruction. 
For a single mode a relevant property of the Wigner function W[g] (a) is expressed by the following formula 

p(x, <j>) = tpixlQlx)^ = í ^W[g] ((x + iy)e**) , (5.80) 
Jr 7r 

which says that the marginal probability obtained from the Wigner function integrating over a genèric direction in 
the complex plane coincides with the distribution of a field quadrature. In conventional medical tomography, one 
collects data in the form of marginal distributions of the mass function m(x, y). In the complex plane the marginal 
r(x, (j>) is a projection of the complex function m(a) = m(x, y) on the direction indicated by the angle <f> G [0, 7r], 
namely 

r(x, <f>) = í m ((x + iyy*) . (5.81) 

.IB TT 
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The collection of marginals for different <f> is called "Radon transform". The tomographic reconstruction essen- 
tially consists in the inversion of the Radon transform ( 15. 81L in order to recover the mass function m(x, y) from 
the marginals r(x, <ft). Thus, by applying the same procedure used in medical imaging Vogel and Risken 1 106 1 
proposed a method to recover the Wigner function via an inverse Radon transform from the quadrature probability 
distributions p(x, <f), namely 

W(x,y) = [* *t [ dx' p{x',<j>) í — \k\ e **(*'-*«* . (5.82) 

Jo t J R J m 4 

In this way the Wigner function, and in turn any quantity of interest, would have been reconstructed by the to- 
mography of the Wigner obtained through homodyne detection. However, this first method is unreliable for the 
reconstruction of unknown quantum states, since there is an intrinsic unavoidable systematic error. In fact the 
integral over k in J5.82Í is unbounded. In order to use the inverse Radon transform, one would need the analytical 
form of the marginal distribution of the quadrature p(x, (j>). This can be obtained by collecting the experimental 
data into histograms and splining these histograms. This is not an unbiased procedure since the degree of splining, 
the width of the histogram bins and the number of different phases on which the experimental data should be 
collected are arbitrary parameters and introduce systematic errors whose effects cannot be easily controlled. For 
example, the effect of using high degrees of splining is the wash-out of the quantum features of the state, and, 
vice-versa, the effect of low degrees of splining is to create negative bias for the probabilities in the reconstruction 
(see Refs. 1102111031 for details). On the other hand, the procedure outlined above allows the reconstruction of 
the mean vàlues of arbitrary operators directly from the data, abolishing all the sources of systematic errors. Only 
statistical errors are present, and they can be reduced arbitrarily by collecting more experimental data. 

5.6 Two-mode entangled measurements 

In this Section we describe in some details three different schemes achieving the joint measurement of the real 
and the imaginary part of the complex normal operators Z± = a±w, a and b being two modes of the field. The 
POVMs (actually spectral measures since Z± are normal) of this class of detectors are entangled, i.e. consist of 
projectors over a set of maximally entangled states, and thus represent the generalization to CV systems of the 
so-called Bell measurement. Detection of Z has been realized in different contexts, e.g. the double-homodyne 
scheme has been employed in the experimental demonstration of C V quantum teleportation [ 571 . 

In the next three Sections we address double (eight-port) homodyne, heterodyne, and six-port homodyne re- 
spectively, whereas in Section l5.6.4l the common two-mode POVM is evaluated. In Section l5.6.4l we also derive 
the single mode POVMs corresponding to situations in which the quantum state of one of the mode is known and 
used as a probe for the other one. 

5.6.1 Double-homodyne detector 

Double homodyne, also called eight-port homodyne, detec- 
tor is known for a long time for the joint determination of phase 
and amplitude of the field in microwave domain, and it was 
subsequently introduced in the optical domain 1 107 1. 

A schematic diagram of the experimental setup is reported 
in Fig. 15.51 There are four input modes, which are denoted by 
a, b, 0,3, and 0,4, whereas the output modes, i.e. the modes that 
are detected, are denoted by Cfe. There are four identical pho- 
todetectors whose quantum efficiency is given by ?/. The noise 
modes used to take into account inefficiency are denoted by 
Uk- The mixing among the modes is obtained through four bal- 
anced beam splitters: three of them (denoted by BS in Fig. l5.5> 
have real coupling Ç = 7r/4, i.e. they do not impose any ad- 
ditional phase, whereas the fourth has evolution operator [see 
Eq. (I1.59H given by 

U(C±) = exp{C±o. + &- (±atf) (5.83a) 

c±= J ex pHf-^±)} > < 5 - 83b ) 

where <j>± = ±ir/2 is the phase-shift imposed by a shifter (a quarter-wave plate) inserted in one arm. We consider 
a and ò as signal modes. The mode 0,4 is unexcited, whereas 0,3 is placed in a highly excited coherent state \z) 



b *■ 



a 



BS 







BS 



c, 1 



BS 



Figure 5.5: Schematic diagram of an eight-port ho- 
modyne detector. 
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provided by an intense làser beam, and represents the local oscillator of the device. The detected photocurrents are 



c\ck, which form the eight-port homodyne observables 



Z^ = 



h-h 
2 V \z\ 



Z 2 = 



HA 



(5.84) 



The latter are derived by rescaling the difference photocurrent, each of them obtained in an homodyne scheme. In 
Eq. (15 . 84i -q denotes the quantum efficiency of the photodetectors whereas \z \ is the intensity of the local oscillator. 
In order to obtain Z\ and Z 2 in terms of the input modes we first note that the input-output mode transformation is 
necessarily linear, as only passive components are involved in the detection scheme. Thus, we can write 



Cfe 



5> 

i=i 



klO-l , 



M = - 
2 



V 



1 
1 

-iS± 
-i6± 



J6± 



1 



-1 
1 

J4>± 



1 \ 

-1 

J<t>± 

J4>± j 



(5.85) 



where ai = a, a 2 = b, 6± = \ — <\>±, and the transformation matrix M can be computed starting from the 
corresponding transformations for the beam splitters and the phase shifter. Eq. J5.85t together with the equivalent 
scheme for the inefficient detection leads to the following expression for the output modes, namely 



('k 



(5.86) 



i=i 



Upon inserting Eqs. J5.86Í in Eq. (I5.84Í . and by considering the limit of highly excited local oscillator, we obtain 
the two photocurrents in terms of the input modes. If we set the phase shifter at tj> + and tune the fourth beam 
splitter accordingly we have 



Zir,+ = q a + Çb + 



Z2r)+ = Pa- Pb + 



while if we choose 0_ we obtain 

Zi v - 



■>2ri- 



q a - qb + 

Pa+Pb + 



1-rj 



1-ri 



i-n 



V 



[q Ul - qu 2 ] + 0(\z\ V ). 



[p Ui ~ Pu 3 ] + 0(\z\ 1 ), 



[q Ul + q U2 ] + 0(\z\ x ). 



[Pu, + Pu 3 ] + Odzl- 1 ) , 



(5.87a) 
(5.87b) 

(5.88a) 
(5.88b) 



where q^ and pu in Eqs. (157: 
assume Ki = 1) 



denotes quadratures of the different modes for specific phases as following (we 
1 



1 



x o = x(a f +o) , P = x 7l/2 = — (a-a f ) . 



Using Eq. ( 15.891 1 we may write the complex photocurrent Z = Z\ 

Z- = a — b^ or Z + = a 



li 

iZ 2 as follows 



(5.89) 



(5.90) 



whereas, for non unit quantum efficiency, it becomes a Gaussian convolution of Eq. J5.90Í . as we will discuss in 
detail in Section l5.6.4l 

It is worth noticing here that the mode transformation defined by Eq. J5.85Í is distinctive for a canonical 4x4- 
port linear coupler as defined in Refs. [ 108 1. It has been rigorously shown 1 109 1 that aiVx TV -port linear coupler 
can always be realized in terms of a number of beam splitters and phase-shifters. However, this implementation is, 
in general, not unique. The interest of eight-port homodyne scheme lies in the fact it provides the minimal scheme 
for realizing a 4 x 4-multiport. 



5.6.2 Heterodyne detector 

Heterodyne detection scheme is known for a long time in radiophysics and it has been subsequently introduced in 
the domain of òptics II 101 . The term "heterodyne" comes from the fact that the involved modes are excited on 
different freqüències. 
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Figure 5.6: Schematic diagram of a heterodyne detection. Relevant modes are pointed out. 



In Fig. l5.6l we show a schematic diagram of the detector. We denote by Es the signal field, whereas E^ç, describes 
the local oscillator. The field E^ accounts for the losses due to inefficient photodetection. The input signal is 
excited in a single mode (say a) at the frequency u, as well as the local oscillator which is excited at a mode at the 
frequency uq. This local oscillator mode is placed in a strong coherent state \z) by means of an intense làser beam. 
The beam splitter has a transmissivity given by r, whereas the photodetectors shows quantum efficiency ?y. The 
heterodyne output photocurrents are given by the real Z\ and the imaginary Z2 part of the complex photocurrent 
Z. The latter is obtained after the rescaling of the output photocurrent /, which is measured at the intermediate 
frequency ujj = u — uq. By Fourier transform of Eq. ( 15. lli we have 

I(ui) = [ deu' E ( -\uj' + lo i )E ( + ) {uj') , (5.91) 

JR 

Eq 1 ^ being the positive and the negative part of the output field. In terms of the input fields Eq. ( 15.911 1 can be 
written as 



du' 



Heterodyne photocurrent is obtained by the following rescaling 

I(ui 

1111 

— 1 I»Iti 

1*1 



rfFEf^ \lü' + ui) + y/ v (l - t)E { l ^(lü' + ui) + y/l- vE { L \u'+u I ) 



lim 



(with | z | a/1 — t constant) 



(5.92) 



(5.93) 



In practice, this definition corresponds to have a very intense local oscillator, which is allowed only for a little 
mixing with the signal mode 11111 . In this limit only terms containing the local oscillator field E^q(uj ) at the 
frequency ui can survive in Eq. ( 15.921 . so that we have 



Zf)~ — Z\ v - + iZ2 V - , 



where 

Zi v + =q a + q c + 

Z2-Q+ = Pa - Pc + 

In writing Eq. ( 15.951 1 we have substituted 



1 — 7] 



[q Ul -q U2 ] + 0(\z\ 1 ). 



\p ul -p U2 ] + 0(\z\ ^ 



E ( s + \2u -u). 



Ui 



E 



(+) 



U 2 



E ( +\2iü - lu) , 



(5.94) 

(5.95a) 
(5.95b) 

(5.96) 



for the relevant modes involved. Since u\ and U2 are not excited, they play the role of noise modes accounting 
for the quantum efficiency of the photodetector. The expression d5.95i for the heterodyne photocurrents is thus 
equivalent to that of Eq. fl5.87l > for the eight-port homodyne scheme. The full equivalence of the two detection 
schemes has been thus proved. Also for heterodyne detection, a simple rearrangements of phase-shifts provides 
the measurement of the complex operators Z + instead of Z_ . 
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5.6.3 Six-port homodyne detector 

A linear, symmetric three-port optical coupler is a straightforward generalization of the customary lossless sym- 
metric beam splitter. The three input modes a^, k = 1, 2, 3, are combined to form 3 output modes Cfc, k = 1, 2, 3. 
In analogy to lossless beam splitters, which are described by unitary 2x2 matrices II 121 . any lossless triple coupler 
is characterized by a unitary 3x3 matrix lli 131 . For the symmetric case we have the form 




(5.97) 

where £ = exp {ï|tt} and each matrix element Thk represents the transmission ampli tude from the fc-th input 
port to the h-Ü\ output port, namely Ch = J2k=i ThkQ,k- Such devices have already been implemented in single- 



mode optical fiber technology and commercial triple coupler are available II 141 . Any triple coupler can be also 
implemented by discrete optical components using symmetric beam splitters and phase shifters only 1 1 13 1. As it 
has already mentioned in Section l5.6.1l this is due to the fact that that any unitary 77i-dimensional matrix can be 
factorized into a sequence of 2-dimensional transformations plus phase-shifts 1 109 1. Moreover, this decomposition 
is not, in general, unique. In Fig. I5.7l we sketch a possible implementation of a triple coupler where the input 
modes are a\ = a, ü2 = b, and a 3 . Experimental realizations of triple couplers has been reported for both cases, 
the passive elements case and the optical fiber one I113II1151 . 




Figure 5.7: Realization of a triple coupler in terms of 50/50 beam splitters (BS) and phase shifters "0". In order to obtain a 
symmetric coupler the following vàlues has to be chosen: cj>i = arccos (1/3) and (f>2 = </>i/2. 



Let us now consider the measurement scheme of Fig. 15.81 1 1 16 1. The three input modes are mixed by a triple 
coupler and the resulting output modes are subsequently detected by three identical photodetectors. The measured 
photocurrents are proportional to I n , n — 1, 2, 3, given by 

1 3 2tt 

In = 4i c n = o ^2 CX P{^n('- k )} a k a l > S n = —(n-1) , (5.98) 

k,l=l 

where a\ = a and = b. After photodetection a Fourier transform (FT) on the photocurrents is performed 

1 3 

J s = FT(7 1 ,/ 2) J 3 ) = -^=^7 n e X p{-^„( S -l)} (s = l,2,3). (5.99) 

* n— 1 

This procedure is a straightforward generalization of the customary two-mode balanced homodyning technique. 
In that case, in fact, the sum and the difference of the two output photocurrents are considered, which actually 
represent the Fourier transform in a two-dimensional space. By means of the identity 

h{s - 1) = i £exp - 1)1 , (5.100) 

n— 1 ^ ^ 

for the periòdic (modulus 3) Rronecker delta <5 3 , we obtain our final expressions for the Fourier transformed 
photocurrents: 



1 



X\ = \ a) a + b^b + a^a^ 



(^a^b + 6^a 3 + a^a^j , X 3 = —= ^a^a 3 + b^a + a^b^J 



12 ^ v - a , , ^ ^ 

I\ gives no relevant information as it is insensitive to the phase of the signal field, whereas Z2 and X 3 are hermitian 
conjugates of each other and contain the relevant information in their real and imaginary part. In the following 
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Figure 5.8: Outline of triple coupler homodyne detectors: the hexagonal box symbolizes the electronically performed Fourier 
transform (FT). 



let us assume a and b as the signal modes and 0,3 fed by a highly excited coherent state \z) representing the local 
oscillator. For large \z\ the output photocurrents are intense enough to be easily detected. They can be combined 
to give the reduced photocurrents 



Z í+ = V3 



Z- 



2+ 



12+13 =q a +q b + 0(\z\- í ), 



2\z\ 



2i\z\ 



(5.102a) 
(5.102b) 



which we refer to as the triple homodyne photocurrents. Again, the complex photocurrent Z + = Z\ + + ÍZ2+ has 
the form Z_ = a + b\ a and b being two modes of the field. 

When accounting for the non unit quantum efficiency 77 of the photodetectors the output modes are written as 
Ch = ^/rj Ylk=i Thk a k + \/ï — V Uh, h — 1, 2, 3, so that the reduced photocurrents are now given by 

Z 



2^3 - ^2 



q a + qb 



Pa-Pb 



1-ry 



n 



[<3ui +Qu 2 ] + 0(\z\~ 1 ), 



(5.103a) 



(5.103b) 



When, as it is the case, the modes are placed in the vacuum, the six-port photocurrents in Eq. d5. 103i leads 
to the same statistics of the eight-port photocurrents in Eq. ( 15.871 1. Indeed, they describe different devices leading 
to the same amount of information on the signal modes. The measurements of Z_ can also be achieved by six-port 
homodyne by a suitable choice of the phase-shifts among the modes. 



5.6.4 Output statistics from a two-photocurrent device 

Although the two pairs of single-mode quadratures [qk,Pk] = V '2, where k = a.b are two modes of the field, do not 
commute, the sum and difference quadratures do, and therefore can be measured in a single experiment. Indeed, the 
three detection schemes analyzed in this Section provide the joint measurement of the operators q a + qb and p a —pb, 
or q a — qb and p a +pb- In turn, the two cases corresponds to the measurement of the real and the imaginary part of the 
complex photocurrents Z± = a ± respectively. In both cases we have that [Z±,z\_] =0= [3?e[Z±], 9m[Z±]], 
i.e. Z± are normal operators, and therefore the spectral theorem holds 

d 2 z z \z)) ±± ((z\ , 

where \z))± with z G C are orthogonal eigenstates of Z±, respectively. 

Let us first consider Z_ = a — b>. Using the matrix notation introduced in Section [T~2l we have that 

|*))_ = 4= \ D (»))) = -7= D ( z ) ® = 4= I® ^(-**)|l)> - (5-104) 

V 7T y/lT y/TT 

where D{z) is the displacement operator and = J2 n \ n ) ® \ n )- ^ n f act 

Z-\z))- = 4= D{z)D\z){a tf)D(z)\l)) 

= 4= D(z)(a + z - &t)|i)) = 4= D{z) z\t)) = z\z))- , (5.105) 
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where we have used the fact that o®I = Orthogonality of |z))_'s follows from that of displacement 

operators 

-<{w\z))- = -Tr \dUw)D(z)] = 5 {2) {z - w) . (5.106) 

7T L J 

Notice that the eigenstates of the complex photocurrent Z + = a + ¥ may be analogously written as 

\z))+ = 4=-D(«)®I|J» , (5-107) 
V 71 " 

where [I\ pq = (— ) p S pq , i.e. |J)) = ) p li 5 ) ® Ip)- If ^ i s tne density matrix describing the quantum state of 

modes a and b the statistics of the measurement is described by the probability density 

K ± (z)=Tr ab [RE ± {z)\ , 

with E±(z) = \z))±±{{z\ denoting the overall POVM of the detector. 

Let us now consider the effects of nonunit quantum efficiency. The measured photocurrents are given in 
Eqs. KKÍ\ . KM . K95Í or Í5~T03l : using (15. 9> the POVM n í) (z) is obtained upon tracing over the vacuum 
modes used to simulate losses: for either (z) we have 

-^ U1U2 «00| exp { 7 (Zt - z*) - 7 *(Z V - z)} |00))„ 1U2 

■ 7T 



jj^ exp{ 7 (Z + -z*)- 7 *(Z_-z)} exp|-i-^| 7 | 2 | 



■ exp < Z — x 



where 



7r(l — 77) [ 1 — 77 
d 2 w ( \z — w\ 2 , 



A, 2 = ^ — 5 . (5.109) 

77 



The characteristic function of the POVM, for unit quantum efficiency, is given by 

X [E_(z)}(\ u \ 2 )=-({z\D(\ l )®D(\ 2 )\z))- 

= I e A lZ *-At, ((i| jD (A 1 ) s í>(Aa)|l)) 

7T 



I e A lZ *-A*. {{1 \ D{Xl )D T (X 2 ))) 
7T 

-e AlZ *- A i 2 Tr[L>(Ai) J D T (A 2 )] 

s^'-^'Í^ÍAi-AS). (5.110) 



Analogously, 

X [E + (z)](\ 1 ,\2) = +«*|X>(Ai) ® £>(Aa)|«)). 



I e Ai**-AJ* ((j|U(Ai)®£>(A 2 )|J» 
}_ e x lZ *~\U ((J\D(X 1 )ID T (X 2 ))) 

7T 
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— e 

7T 



,Ai — Aí 2 



Tr[n^(Ai)nD T (A 2 )] 



= e >i**-AÍ*5(2)(A 1 + ^) , (5.111) 
where II = <£>&(— ) a I- afe = (— )£* a I a fc is the multimode parity operator. Using (15. 1 101 and ( 15.1 1 1> we have 

W[E±(z)](X) = \s{ Xl ±x 2 - x,)6(yi T Vi - V.) , (5.112) 
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where z = x z + iy z and we used the Cartesian form of the Wigner function for the sake of simplicity. 

For nonunit quantum efficiency V^[n i( ±] (X) is given, according to Eq. J5. 108i . by a Gaussian convolution of 

W[E ± (z)](X),i.e. 



W[U V± ](X) 



exp < - 



(x 1 ±x 2 - x z f (y 1 TV2- Vz)' 



2A2 



2A2 



(5.113) 



Let us now consider a situation in which R is factorized, namely R = a ® t. a is the state under investigation 
and t a known reference state usually referred to as the probe of the detector (see Fig. 15.91 1. The statistics of the 
outcomes, for unit quantum efficiency, may be described as follows 



1 



K(z)=Tv ab [a®Tll(z)] = -Tr a aTr h [i ® t \D(z)))((D(z)\] 



= -Tr Q 

TT 



with 



aTr b [\D(z)rn)((D(z)\] = Tr a [a U^z)] 



Hi(z) = -D{z)t t D\z) , 



(5.114) 



(5.115) 



which is the single mode POVM of the detector viewed as a measurement of the first mode probed by the second 
mode 1 1 17]. If t = |0) (0| is the vacuum, then the POVM IIi(z) is the set of (nonorthogonal) projectors \z)(z \ over 
coherent states, and setup measures the Q-function Q{z) = Tï^ 1 (z\a\z) of the state a. Notice that, as required for 
a POVM, IT(z) is selfadjoint and normalized. The first property follows from the fact that r itself is selfadjoint. 
In fact, = r implies that r T = r* and therefore = ■n^ 1 D{z) r* (z) = n~ 1 D(z) t t D^(z) = Tli(z). 

Normalization follows from completeness of the set of displacement operators, and in particular from Eq. dl.52b> . 
The role of signal and probe may be exchanged, and the statistics can be written as follows 



K{z) = -Tr ò 

7T 



= -Tr 6 

7T 



-Tr a [a®l\D{z))){{D{z)\}_ 
-Tv a [\aD(z)))((D(z)\]} =Tr 6 [rn 2 (z)] , 



where the POVM acting on the mode b is given by 



n 2 (z) = -D(-z*)a T D\-z*) = -D T {z)a T D*{z) 

TT 7T 



(5.116) 



(5.117) 



3m[z] 




3ie[z] 



3m[z] 




Figure 5.9: (a): measurement of the two-mode POVM E(z) viewed as a single-mode measurement of the r-dependent POVM 
<5.115> on mode a; (b): the same for the POVM <5.117> on mode b. 



The action of lli (z) and II2 (z) is depicted in Fig. l5.9l (a) and (b) respectively. The Wigner functions of the POVMs 

Ilfc(z), k = 1, 2 are given by 



W[Ui{z)](a) = -W[r T ]{a - z) = -W[r](a* ~ z*) , 

TT TT 

W[n 2 (z)](a) = -W[a T }{a + z*) = -W[r]{a* + z) , 

TT TT 



(5.118) 



where we have used 11.11 Ibi and the fact that transposition corresponds to mirror reflection in the phase space 
[see Eq. ( 13. 7H . For nonunit quantum efficiency the POVMs becomes Gaussian convolutions of the ideal POVM 
(with variance equal to Af ). The Wigner functions modify accordingly. 
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Chapter 6 

Nonlocality in continuous variable systems 



In their famous paper of 1935 |118|, Einstein, Podolsky and Rosen (EPR) introduced in quantum physics two 
strictly related concepts: entanglement 1 and nonlocality, which afterward generated a longstanding debaté on the 
completeness of quantum mechanics. These two concepts have become more and more important in the subsequent 
dècades, as the recent progresses in quantum information science definitely demonstrated. 

In Chapter|3]quantum entanglement for continuous variable (CV) systems has been extensively analyzed. We 
also pointed out that the concept of entanglement coincides with nonlocality only for the simple case of bipar- 
tite pure states. As soon as we deal with bipartite mixed states, entanglement can be found which do not show 
properties of nonlocality (while the converse instead is always true) |51 1. 

This Chapter will be devoted to the issue of nonlocality for CV systems. First of all, we recali what the concept 
of nonlocality means. Usually two different notions are subsumed in it: (non-)locality and (non-)realism. A theory 
is said to be local if no action at distance, between two subsystems A and B, is contemplated in it. Hence, quoting 
Einstein 1 120 1: 

"The real factual situation of the system B is independent of what is done with the system A, which 
is spatially separated from the former." 

A realistic theory is a theory able to assign a definite counterpart to every element of reality and again following 
ÍTTÏÏ1 : 

"If, without in any way disturbing a system, we can predict with certainty (i.e., with probability equal 
one) the value of a physical quantity, then there exist an element of physical reality corresponding to 
this physical quantity." 

A theory which is not a local realistic one is simply incomplete, according to the spirit of EPR paper. Let us now 
apply this notions to a composite system of two distant partides, described by the so-called EPR wave function, i.e. 
the TWB wave function dl.87t in the limit of infinite squeezing 2 . By measuring, say, the position of one particle 
the position of the other one can be predicted with certainty, as follows from the correlations between the two. If 
there is no action at distance, this prediction can be made without in any way disturbing the second particle. Hence 
an element of physical reality must be assigned to its position. On the other hand, the same argument apply to 
the measurement of momenta. However, quantum theory precludes the simultaneous assignment of position and 
momentum without uncertainty. So EPR conclude that quantum theory is not complete. 

The debaté about whether or not quantum mechanics is a local realistic theory remained in the realm of philos- 
ophy, rather of physics, for many years. The situation drastically changed when Bell proved that EPR point of view 
leads to algebraic predictions (the celebrated Bell's inequalities) that are contradicted by quantum mechanics 1 72 1 . 
Bell formulated his inequalities in a dichotomized fashion, suitable for a discrete variable setting rather than for the 
original continuous variable one. In particular, Bell followed the simple and elegant formulation given by Bohm 
to the EPR gedanken experiment using spin-i partides. More recently, however, the increasing importance of CV 
systems leads many authors to explore the nonlocality issue in its original setting, where dichotomic observables 
to test Bell's inequalities are not uniquely determined. The attempts to translate Bell's inequalities to CV clarified 
the fact that crucial in a nonlocality test is the existence of a set of dichotomized bounded observables used to 
perform the test itself, from which the so called Bell operator is derived. The more debated question has dealt with 
the nonlocality of the normalized version of the original EPR state, i.e. the TWB state of radiation. Nonlocality 
of the TWB was not clear for along time. Using the Wigner function approach, Bell argued that the original EPR 

The word "entanglement" was introduced in this contest by Schròdinger in hisreply 119 1 to EPR paper. 

2 Using the expression 12. 411 for the TWB Wigner function it is immediate to see that the original EPR state introduced in Ref. II 181 is 
recovered in the infinite squeezing limit. 
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state, and as a consequence the TWB too, does not exhibit nonlocality because its Wigner function is positive, 
and therefore represents a local hidden variable description 1 121 1. More recently, Banaszek and Wodkiewicz 11221 
showed instead how to reveal nonlocality of the EPR state through the measurement of displaced parity opera- 
tor. Furthermore a subsequent work of Chen et al. \ 123 1 showed that TWB's violation of Bell's inequalities may 
achieve the maximum value admitted by quantum mechanics upon a suitable choice of the measured observables. 
Indeed, the amount of violation crucially depends on the kind of Bell operator adopted in the analysis, ranging 
from no violation to maximal violation for the same (entangled) quantum state. 

6.1 Nonlocality tests for continuous variables 

In this Section we recali the inequalities imposed by local realism to the situations of our interest. Let us start by 
focusing our attention on a bipartite system. Let m(a) denotes a one-parameter family of single-system observable 
quantities with dichotomic spectrum. In the following m(ai) — ±1 and m(a' 1 ) = ±1 will denote the outcomes of 
the measurements on the first subsystem and, similarly, m{a>}) = ±1 and m(a' 2 ) = ±1 for the second subsystem. 
The essential feature of this measurements is that they are local, dichotomic and bounded. The Bell's combination 

B 2 = m(ai) <g> m(a 2 ) + m(ai) <g> m{a' 2 ) + m(a' 1 ) <g> m{az) — m{a' í ) ® m(a' 2 ) , (6.1) 
under the assumption of local realism, leads to the well known Bell-CHSH inequality 1 124|: 

B 2 = \(B 2 )\ = \E( ai ,a 2 ) +E( ai ,a' 2 ) a 2 ) - £K,a 2 )| < 2 , (6.2) 

where E(ai, a 2 ) = (m(ai) ® m(a2)) is the correlation function among the measurement results. If we describe 
the system quantum mechanically, then we have that 

E(ai,a 2 ) = Tr[i2m(ai) ® m(a 2 )] , 

R being the density matrix of the system under investigation. 

Bipartite entangled pure states violate ( 16. 21 for a suitable choice of the observables m(a) and of the vàlues of 
the parameters. Bipartite entangled mixed states may or may not violate d6.2l >. Systems which involves only two 
parties are the simplest setting where to study violation of local realism in quantum mechanics. A more complex 
scenario arises if multipartite systems are considered. Studying the peculiar quantum features of these systems is 
worthwhile in view of their relevance in the development of quantum communication technology, e.g. to manipulate 
and distribute information in a quantum communication network 1461 1491 . Although the study of multipartite 
nonlocality has originated without the use of inequalities |70|, an approach to derive Bell inequalities (so called 
Bell-Klyshko inequalities) has been developed 1 73l 1741 also for these systems and applied to characterize their 
entanglement properties 1 75 1 . Being originally developed in the framework of discrete variables, these multiparty 
Bell inequalities have found application also in the characterization of continuous variable systems 11251 11261 . 
Bell-Klyshko inequalities 17311741 175 1 provides a generalization of inequality J6.2Í and are based on the following 
recursively defined Bell's combination (operator) 

B n = ~ [m[a n ) + m{a' n )] ® B n -i + ^ [m(a„) - m(a' n )] <g> B' n _ 1 , (6.3) 

where B' n denotes the same expression as B n but with all the a n and ot' n exchanged, and m(a n ) = ±1, m(a' n ) = 
±1 denote the outcomes of the measurements on the n-th party of the system. Bell-Klyshko inequalities then read: 

B n = \{B n }\ < 2. (6.4) 

In the case of a three-partite system, local realism assumption imposes the following inequality from combination 
JO : 

B 3 = \E(ai,a 2 , a' 3 ) + E(ai,a' 2 ,a 3 ) + E(a[,a 2 , a 3 ) ~ E(a' 1 , a' 2 , a' 3 )\ < 2 , (6.5) 

where again E(a\ , a 2 , a 3 ) is the correlation function between the measurement results. 

Quantum mechanical systems can violate inequalities Í6.21 and d6.5> by a maximal amount given by, respec- 
tively, B 2 < 2^2 and B 3 < 4. In general B 2 n < 2 n+1 holds (see, e.g., Ref. J75)). 

We now briefly review three different strategies to reveal quantum nonlocality in the framework of continuous 
variables systems. These nonlocality tests are the basis for the analysis the will be performed in the remaining of 
the Chapter. In order to introduce the argument, recali that in the case of a discrete bipartite system, for example 
a spin--| two particle system, the local dichotomic bounded observable usually taken into account is the spin of 
the particle in a fixed direction, say d. Hence the correlation between two measurements performed over the two 
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partides is E(di, cb) = (di • cr <g> d.2 ■ cr), where the operator cr = (cr x , cr y , cr z ) is decomposed on the Pauli 
matrices base and di, d 2 are two unit vectors. The Bell operator is then given by the expression: 

£>2o- = di ■ cr (£> d2 • cr + d'i • cr ® d2 ■ cr + di • cr ® d'2 • cr — d'i ■ cr ® d'2 ■ cr . (6.6) 

Consider now a n-partite continuous variable system. Following the original argument by EPR it is quite 
natural attempting to reveal the nonlocality of this system trying to infer quadratures of one subsystem from those 
of the others. From now on, we will refer to this procedure as a homodyne nonlocality test (H), as quadrature 
measurements of radiation field are performed through homodyne detection. Here we identify the quadrature x^, 
relative to mode k, according to the definition il . 13i . As they are local but neither bounded nor dichotomic, 
quadrature observables are not immediately suitable to perform a nonlocality test based on Bell's inequalities. The 
procedure to make them bounded and dichotomic is quite arbitrary and consist in the assignment of two domains 
D + and £>_ to each observable 1 127|. When the result of a quadrature measurement falls in the domain D± the 
value ±1 is associated to it. Usually the choice D± = K ± is considered, though a choice suitable to the system 
under investigation may be preferable. Considering a bipartite system we can introduce the following quantities 
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where P{x^,x 2 Lp ) is the joint probability distribution of the quadratures x\ and x 2 ,. We can now identify the 
homodyne correlation function E H (-d, ip) as 

£? a (i?, v ) = p + +(4,4)+p__(xi > 4)-p + _(4,4)-p_ + (4,4) , (6.8) 

which can be straightforwardly used to construct the Bell combination B2H of Eq. \6.2\ and to perform the non- 
locality test. The main problem of pursuing such a nonlocality test is that it is not suitable in case of systems 
described by a positive Wigner function, as the TWB state of radiation. Indeed, a positive Wigner function can be 
interpreted as a hidden phase-space probability distribution, preventing violation of Bell-CHSH inequality unless 
the measured observables have an unbounded Wigner representation, which is not the case of the dichotomized 
quadrature measurement described above. In fact P(x$, x 2 ^) can be determined as a marginal distribution from the 
Wigner function. From Eqs. ( 16. 1\ and dó.81 one has 

E H (iï,tp)= dx\dx\ <24+f dx 2 9+ * sgn[44] W (4>4+f ' ' (6 ' 9) 

where the integration is performed over the whole phase-space and without loss of generality we have considered 
D± = Mr 1 . Eq. ( 16. 9> itself is indeed a local hidden variable description of the correlation function, hence obeying 
inequality i6.2\ . 

In order to overcome this obstacle different strategies have been considered by many authors, based essentially 
on parity measurements. Banaszek and Wodkiewicz 1 122 1 have demonstrated the nonlocality of the TWB consid- 
ering as local observable on subsystem k the parity operator on the state displaced by au (hence we will refer to 
this procedure as a displaced parity (DP) test), which is dichotomic and bounded: 

n 

ïl(ct) = <^D k (a k )(-l) nk Dt(a k ). (6.10) 

k=l 

In the above formula, a = (a%, . . . ,a n ), while n k = a\ak and D k (ak) denote the number operator and the phase 
space displacement operator for the subsystem k, respectively. Hence the correlation function reads: 

E DP (a) = (n(a)), (6.11) 

from which Bell's combinations B20P in Eq. (16. 2> and Bj, DP in Eq. ( 16. 5> can be easily reconstructed in the cases 
n = 2,3. The reason why this procedure would be able to reveal nonlocality also in case of quantum states 
characterized by a positive Wigner function is clear using the following relation [see Eq. 11.1031 1: 



W{a) = 




(6.12) 
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Indeed, the analog of Eq. d6.9t is: 

E DP (a) = J d 2n \(J^J W(a)5 (2n) (a- A) . (6.13) 



Since the Dirac-<5 distribution is unbounded, then Ineqs. Í6.2Ï and i6.5\ are no more necessarily vàlid for B2 D p and 

B^DP- 

Another strategy, developed by Chen et al. | 123 1, shares a similar behavior as the one described above, allowing 
to reveal nonlocality for quantum states with positive Wigner function. This type of nonlocality test will be referred 
to as Pseudospin (PS) nonlocality test. It can be seen as a generalization to CV systems of the strategy introduced 
by Gisin and Peres for the case of discrete variable systems 1 53 1, hence, for the case of a pure bipartite system, it 
is equivalent to an entanglement test 1 129 \. Let us consider the following set of operators, known as pseudospins 
in view of their commutation relations, = (s k , s k ,s k ) acting on the fc-th subsystem 



-ï = £ 

n=0 



(\2n + l) kk (2n + 1| - \2n) kk (2n\) , (6.14a) 



s k x ±is k y =2s k ± , (6.14b) 
d k s k = cosi') k s k z +ÚYid k {é Vk s k _+e- lVk s%) : (6.14c) 

where s_ = J2 n \^ n )kk(2n + 1| = (sÍL)t and dfc is a unit vector associated to the angles i) k and (p k . In analogy 
to the spin-i system mentioned above and defining the vector d = (di , . . . , d n ) the correlation function is simply 
given by: 

E PS (d) = (®£ =1 d fcS fc) , (6.15) 

from which the Bell combinations B 2 ps and £>3 PS are evaluated. Also different representations of the spin-^ 
àlgebra have been discussed in the recent literature 1 1 301 1 1 3 li ÍT32 1 . In particular in Ref. 1 132| it has been pointed 
out that different representations lead to different expectation vàlues of the Bell operators. Hence, the violation 
of Bell inequality for CV systems turns out to depend, besides to orientational parameters, also to configurational 
ones. In the following Sections we will also consider the pseudospin operators Ilfe = (11^,11^,11^) taken into 
account in Ref. 1 132 1, which have the following Wigner representation (for k\ = 2~ 1 / 2 ): 

W[U k x ](a k ) = sgn[3?e[a fc ]] , W[n k ](a k ) = -S^ (3te[a fc ]) V ' 



y "~ K ' ' V "~ L K ' J ' »m[a fc ] ' (6.16) 
W[ïl k z ](a k ) = -n5W(a k ), 

where V stands for the "principal value". The correlation function obtained using operators Y\. k will be indicated 

as£ PS (d) = (®£ =1 d fc n fc ). 



6.2 Two-mode nonlocality 

In this Section we will analyze the nonlocality properties of two-mode states. First we concentrate on the TWB 
state of radiation, then we will consider non-Gaussian states and apply to them all the strategies introduced in the 
preceding section. 

6.2.1 Twin-beam state 

As already mentioned, the more debated question concerning nonlocality in continuous variable systems involved 
the TWB state, due to its importance both from an applicative point of view and from a fundamental perspective, as 
it is a normalized version of the original EPR state. Since it is not suitable for homodyne test, the TWB nonlocality 
will be investigated exploiting the DP and PS tests. 

Displaced parity test 

Let us follow Ref. 1 122|. Using Eq. J2.4H for the Wigner function of the TWB and Eq. J6.12> . it is immediate to 
evaluate the correlation function E DP (pt\ , a 2 ) of Eq. ( 16. lli . In Ref. 1 122 1 the following parameterization has been 
considered to construct the Bell combination B 2dp 

ai=a 2 =0, a' 1 = -a' 2 = Vj. (6.17) 
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It follows that 

B 2DP = 1 + 2cxp{-2J r cosh(2r)} - cxp{-4 l 7e 2r } . (6.18) 




Figure 6.1: (a) Plot of the combination Bïdp defined in Eq. 16. 1 8i and (b) according to the parametrization given by Eq. <6.2QL 
Only vàlues exceeding the bound imposed by local theories are shown. 

As depicted in Fig. l6.lf a). B 2dp in J6.18> violates the upper bound imposed by local theories. For increasing r, 
the violation of the Bell's inequality is observed for smaller J . Therefore an asymptotic analysis for large r and 
J < 1 may be performed. Then a straightforward calculation shows that the maximum value of B2 (for this 
particular selection of coherent displacements) is obtained for 

Je 2r = ~\n2, (6.19) 

corresponding to B2 D p = 1 + 3 • 2~ 4 / 3 w 2.19. Thus, in the limit r — > 00, when the original EPR state is 
recovered, a significant violation of Bell's inequality takes place. Notice that in order to observe the nonlocality of 
the EPR state, very small displacements have to be applied, decreasing as J oc e~ 2r . As pointed out in Ref. II 1221 
the results above have been obtained without any serious attempt to find the maximum violation. For this purpose 
one should consider a general quadruplet of displacements. An analysis to obtain the maximum violation of Bell 
inequalities within this formalism is performed in [ 133 1. Choosing 

ai = -a 2 = iVjct[ = -a 2 = -3iy/j (6.20) 

an asymptotic violation of B 2 — 2.32 can be obtained (see Fig. l6.li b)). This shows that the EPR state does not 
maximally violate Bell's inequalities in a DP test. The reason for this has been addressed in Ref. 11291 . and it is 
attributed to the fact that the displaced parity operator does not completely flip the parity of the entangled quantities 
characterizing the TWB (Le., the number states). 

Pseudospin test 

Let us now focus on the "pseudospin nonlocality test". Considering a TWB state, it is known that the correlation 
function (16 . 1 5 i has the following expression (setting to zero the azimuthal angles) 11231 : 

Epsfòij&z) = costfi cosí9 2 + /twb siní?i sin0 2 , (6.21) 

where / TOB = tanh(2r). Choosing i?j = 0, = ir/2 and $ 2 = — §' 2 , we have 

B 2 = 2(cos 2 + / TWB sin a ) , (6.22) 

and, for this specific setting, the maximum of B 2 is 

B 2 = 2y/l + /twb ■ (6.23) 

It turns out that the TWB state violates the Bell's inequality i6.2\ for every r ^ 0. The violation increases 
monotonically to the maximum value of 2%/2 as the function / TWB — * l,i.e., as the squeezing parameter r increases. 
This indicates that the EPR state maximally violate Bell's inequality. Furthermore, may be regarded as a 
quantitative measure of quantum nonlocality. 

In Ref. 11321 different representations of SU(2) àlgebra have been considered to exploit nonlocality of the 
TWB. In particular using the operators given in Eq. ( I6.16> one can show that the correlation £ PS is still given by 
Eq. (I6.2H . where now the function / TWB is substituted by 

2 

/twb — — arctan[sinh(2r)] . (6.24) 
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Therefore, the behavior of the Bell combination B2 is the same as above, but for any squeezing parameter r it 
gives a lower violation of local realism if compared to Eq. (I6.23> . In general, it is possible to demonstrate that 
the configurational parameterization given by Eq. (I6.14> leads to maximal violation for all vàlues of r. Finally, we 
mention that besides the representations of SU(2) given by Eq. J6.141 and Eq. dó.lót . different ones may be found 
for which the Bell combination B2 is not even a monotonic function of r, i.e. is not a monotonic function of the 
entanglement. 



6.2.2 Non-Gaussian states 

As already pointed out in Section l5~4l it is expected that non-Gaussian states are characterize by a larger nonlocality. 
Let us now exploit this possibility using the non-Gaussian states introduced in Section l5~4l Since the Wigner 
function of IPS and TWBA states is non-positive, all the nonlocality tests introduced will be considered. 



Displaced parity test 



In addressing nonlocality of IPS state, we will consider both the parameterizations \6.\1\ and (I6.20Í . We denote 
B(J) = Bidf for parameter Eq. (16. 171 . and C{J) = B2 D p for parameter Eq. (I6.20Í . As for a TWB, the violation 
of the Bell's inequality is observed for small r 11221 . For the rest of the section, we will refer to B(J) as £? |TWB) (j7) 
when it is evaluated for a TWB J5.28Í . and as B' WS) (J) when we consider the IPS state (15.481 . 



2.25 
2.2 

2.15 
B 2.1 

2.05 
2 

1.95 





0.02 0.04 0.06 0.08 0.1 



Figure 6.2: Plot of B( J) for J = 1(T 2 . The dashed line is B'™ 3 '^), while the sòlid lines are B (IPS, ( J) for different vàlues 
of T cfi (see the text): from top to bottom r cff = 0.999, 0.99 and 0.9. When r cff = 0.999, the maximum of B {WS \J) is 2.23. 
The plot on the right is a magnification of the region < r < 0.11 of the upper one. Notice that for small r there is always a 
region whereB (TWBI (.7) < B aes \J). 



We plot B (mli) (J) and B <IPS) ( t 7) in the Figs. l6.2l and l6.3l for different vàlues of the effective transmissivity t c s and 
of the parameter J: for not too big vàlues of the squeezing parameter r, one has that 2 < B'™ B) (J') < B (WS) (J). 
Moreover, when r e g approaches unit, i.e. when at most one photon is subtracted from each mode, the maximum of 
Bd?s) is always greater than the one obtained using a TWB. A numerical analysis shows that in the limit r c ff — > 1 
the maximum is 2.27, that is greater than the value 2.19 obtained for a TWB 1 122 1. The limit r c g — > 1 corresponds 
to the case of one single photon subtracted from each mode 1931 1941 . Notice that increasing J reduces the interval 
of the vàlues of r for which one has the violation. For large r the best result is thus obtained with the TWB since, as 
the energy grows, more photons are subtracted from the initial state 1 92 1 . Since the relevant parameter for violation 
of Bell inequalities is r e ff, we have, from Eq. (15.441 . that the IPS state is nonlocal also for low quantum efficiency 
of the IPS detector. 
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Figure 6.3: Plots of B(J) as a function of the squeezing parameter r for two different vàlues of J: (a) J — 5 ■ 10 2 and 
(b) J = 10~\ In all the plots the dashed line is B [ ™ B) (J), while the sòlid lines are B aeS) (J) for different vàlues of r cff 
(see the text): from top to bottom T e s = 0.999, 0.9, 0.8, 0.7 and 0.5. Notice that there is always a region for small r where 
B <TWB \J) < B i,es, (J). When t cB = 0.999 the maximum of B m \J) is always greater than the one f B irxY "(J). 
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The same conclusions holds when we consider the parametrization of Eq. ( 16.201 . In Fig. l6.4l we plot C {Tmí) {J) 
and C tles '(J'), i.e. C(J) evaluated for the TWB and the IPS state, respectively. The behavior is similar to that of 
B(J), the maximum violation being now C" rs '(J) = 2.43 for r cff = 0.9999 and J = 1.6 • 10~ 3 . Finally, notice 
that the maximum violation using IPS states is achieved (for both parameterizations) when r c ff approaches unit 
and for vàlues of r smaller than for TWB. 
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Figure 6.4: Plots of C(J) as a function of the squeezing parameter r for J = 1.6 ■ 10 -3 . In all the plots the dashed 
line is C ( ™ B> (J'), while the sòlid lines are C VPS) (J') for different vàlues of r e fi (see the text): from top to bottom r e fl = 
0.9999, 0.999, 0.99 and 0.9. When r cft = 0.9999 the maximum of C iWS '{J) is 2.43. 

Concerning the TWBA J5.54Í . let us consider the case of large ÍV2 and small N$, say 1V3 = 10 — ^JVa) -1 . As 
in the analysis of the entanglement properties of the tripartite state \T), the phase coefficients <j>2 and ^3 play no 
role in the characterization of nonlocality. Using the parametrization ai = \ot2 = h&'i = i\ÍJ and a' 2 = 0, an 
enhancement of the violation of Bell's inequality can be observed with respect to the TWB. Indeed the asymptotic 
violation turns out to be of B-2 D p = 2.41. It can be found, for large N2, when JN2 — 0.042 (see Fig. 16.51 3 . 




Figure 6.5: Bell combination obtained choosing optimized displacement parameters for TWBA state 15.541 (see text for 
details). Only vàlues violating inequality 16. 2t are shown. 

Although the IPS and the TWBA states allow for an enhancement of nonlocality with respect to the usual 
TWB state, they never reach the maximum violation admitted by quantum mechanics. As already pointed out the 
reason for this can be attributed to the fact that the displaced parity operator does not completely flip the parity 
of the entangled quantities characterizing the three states above (i.e. the number states). However, the maximum 
violation of the Bell's inequality in the contest of a DP test could be achieved if the following state |ECS) (entangled 
coherent state) 1 134| could be produced experimentally 

|ECS)=JV(|7>|-7>-|-7>l7», (6-25) 
where M is a normalization factor and I7) is a coherent state with 7 7^ 0. Its Wigner function read as follows: 

W ECS (a, (3) = AN 2 [ cxp{-2|a - 7 | 2 - 2\0 + 7 | 2 } + cxp{-2|a + 7 | 2 - 2|/3 - 7 | 2 } 
- exp{-2(a - 7 )(a* + 7) - 2(/3 + 7) 09* - 7) - 4 7 2 } 

- exp{-2(a* - 7 )(a + 7) - 2{f3* + 7) 09 - 7) - 4 7 2 }} , (6.26) 

where 7 is assumed to be real for simplicity. Notice that the ECS state may be represented in the 2 x 2-Hilbert 
space as 

|ECS) = ^=(\e)\d)-\d)\e)), (6.27) 




3 The same analysis holds if we reverse the role of the two modes, provided that the conditional measurement to obtain the TWBA is 
performed on mode 02 of the original tripartite state, rather then on mode 03 . 
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0.5 1 1.5 2 2.5 3 

r 



Figure 6.6: Plots of Bps' the dashed line refers to the TWB, whereas the sòlid lines refer to the IPS with, from top to bottom, 
r cff = 0.9999, 0.99, 0.9, and 0.8. 



where |e) = A/+(|7) + | — 7)) and \d) = Aí-(\j) — \ — 7)) are even and odd states with normalization factors 
Af+ and A/1 . Note that these states form an orthogonal basis, regardless of the value of 7, which span the two- 
dimensional Hilbert space. For this state the maximum violation can be achieved due to the fact the displaced parity 
operators act like an ideal rotation on the even and odd microscopic states |e) and \d) in which the ECS state may 
be decomposed [see Eq. (I6.27H . As a consequence the parity of |e) and \d), which are the orthogonal entangled 
elements in the entangled coherent state, can be perfectly fiipped by the displacement operator (for 7 — * 00), 
allowing for the maximum violation of Bell's inequality 1 129|. 



Pseudospin test 

Now we investigate the nonlocality of the IPS state by means of the pseudospin test considering the pseudospin 
operators given in Eqs. ( I6.16l >. If we set to zero the azimuthal angle, the correlation function (16. 151 reads 

l 4 

EÇ«Wi A) = 7 r V C k [ cos cos i? 2 + /„ sin i? a sin 6 ] , (6.28) 

Pn{r, Tes) ^ 



where we defined 



fe=i 



arctan 



2Br r 



■xA k V v A 

with 

A k = (b- f k ){b - g k ) - {2Bt cS + h k ) 2 , 

and all the quantities appearing in Eq. J6.28I have been defined in Section l5.4.1l 

In Fig.Élwe plot B PS for the TWB and IPS; we set #1 = 0, tí[ = tt/2, and i9 2 = -0' 2 = 7r/4. As usual 
the IPS leads to better results for small vàlues of r. Whereas B ( ™ B> — > 2^/2 as r — > 00, B£s ! has a maximum and, 
then, falls below the threshold 2 as r increases. It is interesting to note that there is a region of small vàlues of r 
for which £>™' < 2 < B^, i.e. the IPS process can increases the nonlocal properties of a TWB which does not 
violates the Bell's inequality for the pseudospin test, in such a way that the resulting state violates it. Note that the 
maximal of violation for the IPS occur for a range of vàlues r experimentally achievable. 

Concerning the TWBA state, a straightforward calculation shows that an expression identical in form to 
Eq. (I6.2H . where the following function / TWBA can be identified: 



^2 (l + iV 3 r?) 



1 + JVi N 3 (1 + Ni) rj 



^ (2/ri 1 /) 1 V 2A- + I ' '"'Vi X\ J V 1 + A'i 

k,p—0 



(i-(i-vT) — T^r ■ (6-29) 



2 k 



In order to compare the violations in the case of the TWB and the TWBA, let us fix as in the previous subsection 
a small value for N3, A plot of the functions / TWB and /twba versus the total number of photons of the TWB, and 
the total number of photons of the initial three-partite state, is given in Fig. J6.7b . It can be seen that the TWBA 
reaches large violations for smaller energies with respect to the TWB. As for the TWBA state, also in the ECS 
case an expression identical in form to Eq. J6.2 li may be found, where now the following function / ECS can be 
identified: 

( 00 4n+l \ ~ 2 

^ V(2n)!(2n + l)!; 

A remarkable feature of this case is that it allows for a maximum violation not only when 7 — > 00, but also when 
7^0. 
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Figure 6.8: Plots of B2 as a function of tanhr: (a) for different vàlues of r e a and for ideal homodyne detection (i.e. with 
quantum efficiency r]n = 1): from top to bottom r e fi = 0.99, 0.95, 0.90, 0.80 and 0.70; (b) with r e s ~ 0.99 and for different 
vàlues of the homodyne detection efficiency r^: from top to bottom r/u — 1,0.95,0.90,0.85 and 0.80. The maximum of 
the violation decreases and shifts toward higher vàlues of r as 77H decreases. For smaller vàlues of r e fi the violation is further 
reduced. 



Homodyne test 

The negativity of the Wigner function that may occur for non-Gaussian states suggests to perform a nonlocality test 
based upon a homodyne detection scheme. Indeed, homodyne nonlocality test has attracted much attention in the 
recent years 11271 11351 1971 1981 |95l . in view of the high quantum efficiency achievable with homodyne detection, 
which offers the possibility of a loop-hole free test of local realistic theories [127|. As we seen, the positivity of 
a Wigner function prevents the violation of homodyne Bell inequality (16. 2\ . On the other hand, its negativity is 
not sufficient, in general, to ensure a violation. Quantum states with negative Wigner function that doesn't violate 
local realism with homodyne test are given for example in Refs. 1135111361 . As shown in Ref. 1 127| also the ECS 
state does not allow for any violation unless it is subjected to an additional squeezing. The same situation occur if 
the TWBA is considered IT571 . 

Concerning the IPS state if one dichotomizes the measured quadratures as described in Section (16. H the Bell 
parameter reads B?. H = Efói, ipx) + J3(i?i, (pï) + E(§2i fi) ~ E('&2, <pï) where dh and iph are the phases of the 
two homodyne measurements at the modes a and b, respectively. Eq. (16. 8> can be rewritten as 

E(<& h ,<Pk)= dx# h dx Vk sign[x# h x tpk \P(x# h ,x lph ) , (6.30) 
Js. 2 

P(x$ h , x Vk ) being the joint probability of obtaining the two outcomes x$ h and x Vk , 

In Fig. 16.81 fa) we plot B2H for i?i = 0, $2 = 7r/2, fi — — tt/4 and ip2 = tt/4: as pointed out in Ref. |98|, 
the Bell's inequality is violated for a suitable choice of the squeezing parameter r. Moreover, when r e g decreases 
the maximum of violation shifts toward higher vàlues of r. As one expects, nonunit quantum efficiency rj H of the 
homodyne detection further reduces the violation [see Fig. l6.8l (b)l. Notice that, when f]n < 1, violation occurs for 
higher vàlues of r, although its maximum is actually reduced: in order to have a significant violation one needs a 
homodyne efficiency greater than 80% (when Tes = 0.99). On the other hand, the high efficiencies of this kind of 
detectors allow a loophole-free test of hidden variable theories 1 127 128 1, though the violations obtained are qui te 
small. This is due to the intrinsic information loss of the binning process, which is used to convert the continuous 
homodyne data in dichotomic results 1 135 1. Better results, even if the violation is always small, can be achieved 
using a circle coherent state 1127111281 or a superposition of photon number states [ 135 1, while maximal violation, 
i.e. B2H = 2\/2, is obtained by means of a different binning process, called root binning, and choosing a particular 
family of quantum states II138II1391 . 
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6.3 Three-mode nonlocality 

As we have seen in Section lóTTI nonlocality of multipartite systems may be analyzed by means of Bell-Klyshko 
inequalities. In particular, we have explicitly considered the constraints (16. 5> that every tripartite system must 
respect in order to be described by a local realistic theory. The aim of this section is to analyze the violation 
of Ineq. (16. 5> in tripartite systems. Both the states introduced in Section 1241 will be considered, as well as the 
parity-entangled GHZ state introduced in Ref . II 1261 . 



6.3.1 Displaced parity test 

Let us start our study of nonlocality for tripartite systems using the displaced parity test. Considering the correlation 
function E DP (a) given by Eq. ( 16.1 1> . the state V3 (i.e., the state whose covariance matrix V3 is given by 
Eq. (12.46» was found in 1 125 1 to give a maximal violation of Bj, DP ~ 2.32 in the limit of large squeezing and 
small displacement. The study in 1 125 1 however was performed for a particular choice of displacement parameters: 
«i = "2 = = and a[ — a' 2 = a' 3 = i\fj . One can identify a number of parameterizations that allow a 
significantly higher violation of Bell's inequality 111371 . Consider the one given by a\ = 012 = 013 = i\fj and 
a\ = a' 2 = a 3 = —2i^fj from which follows that 

B 3DP = 3exp{-12e~ 2, \7} - exp {24e 2r l 7} . (6.31) 

The asymptotic value B30P = 3 is found for large r and J ^ [see Fig. 16.91 (a)l. The importance of a suitable 
choice of the displacement parameters is apparent if this asymptotic value is compared to the violations obtained in 
Ref. 11251 . In that work in fact generalizations to more than three modes of state V3 were also considered, giving 
an increasing violation of Bell inequality as the number of modes increases, but never founding a violation greater 
than 2.8. 




Let us now consider the tripartite state \T), the correlation function is now given by Eq. ( 16. lli with the co- 
variance matrix Vt in Eq. ( 12.511 1. The symmetry of the state suggests a maximum violation of Bell inequality for 
N2 = N3 = jN [recali Eq. ( I2.50H . while the fact that the separability of the state doesn't depend on the phases <f>2 
and 03 suggests that they are not crucial for the nonlocality test. If we consider the same parametrization leading 
to Eq. (I6.3U and fix 4>2 — 4>3 — 7r ' we And: 

_j , g 6 J (l+W+2 V2y/N(2+N)} 2 e f J (4+7 AT+6V2 -0V (2+JV)) 

&3DP = 7 , =v , (6.32) 

from which follows an asymptotic violation of Bell's inequalities of B^dp = 2.89, for large N and small J . A 
slightly better result is found if a parametrization, more suitable and numerically optimized for the state |T), is 
considered: a.\ = \\f~J , 012 — 013 = a\ = 0, a' 2 = —\fj, a'3 = \fj, 4>2 = and <f) 2 = n. The Bell combination 
B30P for this choice of parameters is depicted in Fig.|^Jb). We note that in this case a larger choice of angles 
allows the violation of Bell inequality if compared with Fig. l6.9l (a). The asymptotic violation of Bell's inequality 
is now Bzbp = 2.99. Comparing the results obtained for the two states V3 and |T) it is possible to show that, 
even if the two states have qui te the same asymptotic violation, state V3 reaches it for lower energies 11371 . 



6.3.2 Pseudospin test 

Consider now the pseudospin nonlocality test. Let us calculate the expectation value of the correlation function 
(16. 151 for the state \T) (for simplicity we consider <f> 2 — <j>3 — 0). The only non vanishing contributes are given 
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by: 



ci = (4 ® 4 ® 4) = (4 ® 4 ® 4) 



^ $>*JV? (2S + 2< + 1)! , (6.33a) 



2(1 + Ni) * (2s)!(2í)! A /(2a + l)(2í+l) 

c 2 = (4 ® ® 4) = -{s\ ®s 2 z ® s 3 y ) 

- ^ ^^yffg^/E + g + I, (6.33b) 



2(1+ JVi)^ 2 à (2s)!(2í)! V 2Í+1 



c 3 = (4 ® 4 ® 4) = -<4 ® 4 ® 4) 



VJ^ ^ ., 2s ., 2t (2s + 2í)! / 2s + 2í+l 

Z.^ AA 3 r - vr -. v W - . 1 , (6.33c) 



2(l+JVi)^ z J (2s)!(2í)! V 2s + l 

with JVk = Nk/ (1 + Ni), and by (s), (g> 4 <S> 4) = 1- The correlation function then, according to Eqs. J6. 14b and 
J6. 15i . reads as follows: 



Ep S (d) = cos$i cos #2 cos $3 + ci cos#i sin$ 2 sini? 3 cos(v?2 — ^3) 

+ C2 cos i?2 sin #1 sin # 3 cos(y 1 + ^3) + C3 cos #3 sin $1 sin 



1Í2 cos(<^i - ip 2 ) ■ (6.34) 



Hence, without loss of generality, we can fix for example ipi = and (^2 = <^3 = tt and look for the maximum 
violation of Bell inequality A6.5Í constructed from Eq. ( 16.34b . Notice that if the coefficients Cfe, k = 1,2,3, 
were equal to 1 then the maximum violation admitted, Bj, PS — 4, should be reached. Considering Eqs. ( 16.331 1 
two limiting cases can be studied. First, for large N2 and small N$ (or viceversa) a numerical evaluation of the 
coefficients Ck shows that C3 — > 1, while the other two vanish. Hence, considering $3 = 0, the correlation function 
J6.34Í reduces to that of a TWB subjected to a pseudospin nonlocality test [see Eq. ( K5.21H . allowing an asymptotic 
violation of £>3 PS = 2 y/2. This result should be expected, since in this case the state ( 12.491 1 reduces to a TWB for 
modes ai and 02, while mode 0,3 remains in the vacuum state and factors out. Consider now the case in which 
N2 = N3 = jN. A numerical evaluation shows that the coefficients C4 — * ~ for large N, hence also in this case 
the maximum violation cannot be attained. The asymptotic violation turns out to be B% PS = 2.63. 

As already mentioned in Section 16.11 other representations for the pseudospin operators can be considered. 
Using Eqs. dó. 16l > and the Wigner function associated to state \T) it is possible to calculate the correlation function 
£ PS (d). Setting again the azimuthal angles ipk = 0, the latter shows the same structure as E PS (d) where now the 
coefficient Ck are replaced by 



2 arctan 



N 



2y/l + N ) j _ j _ 2 arctan y/N 



1 n(l + N) ' 20 tt(1 + ÍJV) ' 

An appropriate choice of angles leads to a violation of Bell's inequality given by 03 PS = 2.22, which is now 
reached for N ~ 1, value for which the coefficients c' k are approximately near their màxima. As already pointed 
out, we see that different representations of the pseudospin operators give rise to different expectation vàlues for 
the Bell operator. 

Applying now the same procedure to state V3 we find the same structure for the correlation function E' ps , 
where the coefficients are now given by 



—6 arctan 

Cn = Cq = 



4 coshr sinhí 



n y/5 + 4 cosh(4r) 



(6.36) 



After an optimization of the angles we obtain a maximal violation of £>3 PS = 2.09, for r ~ 0.42 (N ~ 0.56) 
that maximizes the coefficients c^. 

Finally, one may consider the nonlocality issue in the general case of an n-party system. We recali that in the 
case of discrete variable systems Mermin 1 73 1 showed that the multipartite GHZ state, defined as 

|GHZ)„ = ^(l+>i • ■ ■ l+>« - h)i • ■ • h>«) , (6-37) 

where \+)k is the eigenvector with eigenvalue +1 of the Pauli matrix <r z relative to the fc-qubit, admits a violation 
of local realism that exponentially grows with the number of party. The first attempt to compare this behavior with 
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continuous variables case was performed in Ref. 11251 . There, the violation of local realism by the states V n , 
a straightforward generalization to n modes of state V3, has been analyzed. Considering the Bell combination 
B n given in Eq. (I6.4Í in a DP setting, it was found that the degree of nonlocality of states V n indeed grows 
with increasing number of parties. This growth, however continuously decrease for large n, as opposite to the 
qubit case. Nevertheless, as already pointed out in Ref. 1 125 1, this analysis was performed for a particular choice 
of displacement parameters a, a' which unfortunately seems not to be optimal. In fact, the maximum violation 
attained with that choice neverreach the asymptotic value of B^ P3 = 3 obtained with Eq. J6.3H (e.g., B$5 = 2.8 in 
1 125 1). Another approach has been pursued in Ref. HI 261 . where eigenstates of the pseudospin operator s z (parity- 
entangled states) has been considered, in direct analogy to the n-party GHZ states J6.37b . Due to this analogy it 
is straightforward to show that these states lead to an exponential increase of the violation of local realism. For 
example in the tripartite case a maximum violation 63 = 4 can be found. Hence a behavior identical to that of the 
discrete variable systems is recovered. However, recali that to our knowledge there is no proposal concerning a 
possible experimental realization of the parity-entangled (non-Gaussian) states. 



Chapter 7 

Teleportation and telecloning 



In this Chapter we deal with the transfer and the distribution of quantum information, i.e. of the information 
contained in a quantum state. At first we address teleportation, i.e. the transmission of an unknown quantum 
state from a sender to receiver that are spatially separated. Teleportation is achieved by means of a classical and 
a distributed quantum communication channel, realized by a suitably chosen nonlocal entangled state. Indeed, 
quantum teleportation has no classical analog: the use of entanglement permits to transmit an unknown signal 
without classically broadcasting the whole information about its quantum state. On the other hand, we have that 
quantum information cannot be perfectly copied, even in principle. The no-cloning theorem follows from the 
linearity of quantum mechanics 11401 114Í1 . and forbids the existence of any device producing perfect copies of 
genèric unknown quantum states. Only approximate clones can be realized, that can be subsequently distributed 
in a quantum network by means of teleportation. Alternatively, the entire process can be realized nonlocally 
exploiting multipartite entanglement which is shared among all the involved parties. The latter process is known 
as telecloning, and will analyzed in the second part of this Chapter. 



7.1 Continuous variable quantum teleportation 



In this Section we address continuous variable quantum teleportation (CVQT), where the goal is teleporting an 
unknown state <j\ of a given mode 1, from Alice, the sender, to Bob, the receiver, i.e. reconstructing the quantum 
state onto another mode, on which Alice has no access. In the following we refer to the CVQT protocol sketched 
in Fig. 17.11 Alice and Bob share an entangled two-mode state of radiation described by the density matrix p23, 
where the subscripts refer to modes 2 and 3, respectively: mode 2 is own by Alice, the other by Bob. In order to 
implement the teleportation, Alice performs a joint measurement, i.e. the measurement of the normal operator Z 
on modes 1 and 2, getting as outcome a complex number z (see Section l5.6.1> : then, she sends her result to Bob 
via a classical communication channel. Once received this classical information, Bob applies a displacement D(z) 
to his mode 3 and obtains a quantum state g out which, in the ideal case, is identical to the initial state a\ 11421 . 

The original proposal for teleportation concerned states 
in a bidimensional Hilbert space |143|. The correspond- 
ing experiments have been performed in the optical domain, 
using polarization qubit 1 144j 11451 . and for the state of a 
trapped ion 11461 . Also CVQT can be realized by opti- 
cal means 11421 . and successful teleportation of coherent 
states has been realized [27 1. In optical CVQT, entangle- 
ment is provided by the twin-beam state (TWB) of radiation 
|A)) 23 = VT^En=o Xn \ n h\n) 3 , 2 and 3 being two 
modes of the field and A the TWB parameter. We assume 
A as real |A| < 1. TWBs |A))23 are produced by optical 
amplifiers (see Section [T.4.4> and, being a pure state, their 
entanglement can be quantified by the excess von Neumann 
entropy. We refer to Section lXTl for details and just remind 
that the degree of entanglement is a monotonically increasing function of A (or equivalently of the average number 
of photons). As we will see, the larger is the entanglement the higher (closer to unit) is the fidelity of telepor- 
tation based. There are different ways to describe CVQT 1 1421 11471 11481 . but, in general, two of them are the 
most common: the first makes use of photon number-state basis, the other is in terms of Wigner functions. This 
Section addresses the description of CVQT protocol following these two approaches and, in particular, we derive 
the completely positive (CP) map C describing the teleportation process also in the presence of noise. 



Bob 




Figure 7.1: Continuous variable quantum teleportation: 
scheme of the optical realization. 
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7.1.1 Photon number-state basis representation 



Here we describe the teleportation protocol in the photon number basis. The state Alice wishes to teleport to Bob 
is described by the density matrix 



(7.1) 



while the pure two-mode state they share is (in general) 

023 = |C)) 2 332 



(7.2) 



h.k 



where we used the matrix notation introduced in Section [L2l 

The first step of the protocol consists in Alice's joint measurement on modes 1 and 2, which corresponds to 
the measure of the complex photocurrent Z = a% + a\ (see Section l5.6.1l· The whole measurement process is 
described by the POVM 



n 12 (z) = -D í {z)\l))m l {{l\D\(z), 



(7.3) 



where |l))i2 = X™ |^)2, and Z)i(z) is a displacement operator on mode 1 (see Section l5.6.4> . The conditional 
state of mode 3 is then 



03 0) 



1 
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1 



Tr i2 [cri^te ni 2 (z)(g)I 3 ] 



(7.4) 
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(7.5) 



x |fc) 33 (m| ® |/i) 22 (w| ® |p)u<w|£>í(2) 



(7.6) 



^3 Y YY a P1 Chk C «™ l^PlWI")! 1< / ^I- C> í(^)lí')l |fc>33H 

-4-rC T D t («)CTD(2)C* ) 
np{z) 



(7.7) 



where (• • • ) T denotes transposition, the subscripts have been suppressed, andp(z) is the double-homodyne prob- 
ability density, given by 

p(z) = Tr 123 [di ® p 23 n 12 (z) <g> I 3 ] • (7.8) 

After the measurement, Alice sends her result to Bob through a classical channel and, then, he applies a displace- 
ment D(z) to his mode, in formula: 

g 3 (z) -> q' 3 (z) = D{z) g 3 (z) D\z) . (7.9) 
If the entangled channel is provided by the TWB, then C = (1 — A 2 ) 1 / 2 A a a and Eq. Í7.9Í rewrites as 



= ^ ,V A ata £> + (*) t D(z) X aia DUz) . 

TTp{z) 

Now, using the operatorial identity (II .491 , Eq. J7. 101 can be reduced to 



(7.10) 



(1-À 2 ) f d 2 wd 2 v 



TT P (Z) y C 2 7T 2 (l-A)' 



exp 



|^(M* + M*) 



x D{z) D{w) D\z) a D{z) D\v) D\z) , (7.11) 
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which, thanks to (ll.43>. becomes 



(1 + A) 
irp(z) 



d w d v 
tt 2 (1 - A) 



exp 



1 1 + A 

2 1 - A 



(h 



x exp{(u> — v)*z — (w — v)z*} D(w) a D\v). (7.12) 



The final output of CVQT is obtained integrating g' 3 over all the possible outcomes z of the double-homodyne 
detection, i.e. 



(7.13) 



ÉW = / d 2 zp(z) g' 3 (z) , 
Jc 

and, remembering the definition Jl .441 of the delta function S^ 2 \Ç), we obtain 



f?out 



W «Pj"£r^ D(w)aDÏ(w), 



(7.14) 



with 



a 2 _ = 



1 1- A _ 1 

41 + A = 4 6 



-2r 



r = tanh - A being the squeezing parameter of the TWB. Eq. il.141 corresponds to an overall Gaussian noise with 
parameter a 2 _ (see Section l4~2t : in this way the CVQT protocol can be seen as a thermalizing quantum channel 
1 149 1. Notice that g out approaches the input state a only for A — > 1 (or r — » do), i.e. for a TWB with infinite 
energy. 



7.1.2 The completely positive map of CVQT 



CVQT corresponds to a Gaussian completely positive (CP) map and, as we will see in Section l7.1.6l this result 
holds also in the presence of noise. If g- m is the state at Alice's side, the state at Bob's side will be 



gout — £-Q'u 



d 2 w FX , . = — D{w) g in D\w) , 



(7.15) 



w denoting the vector (5fte[w] , 3m[w]) T , and S is a 2 x 2 matrix describing a Gaussian noise, as we have addressed 
in Section l4~2l Notice that Eq. (17. 15i provides already the Kraus diagonal form of the teleportation map. In Section 
I7.1.6l we will explicitly derive the map (17 . 1 5 i for teleportation in the presence of noise. In the case of Eq. J7. 14i 
one has 

' o 2 _ 
a 2 _ 



S = Sn = 4 



(7.16) 



7.1.3 CVQT as conditional measurement on the TWB 

As we have seen in Section l5.6.4l when the modes in the measurement of Z are initially excited in a factorized 
state, then we can write the POVM as a single-mode POVM depending on the state of the other mode. This is the 
case of CVQT, which can be seen as the measurement of the POVM [see Eq. ( 15.1 17H 



U 2 (z) = -D^(z)<7 T D* 2 (z). 
ir 



(7.17) 



acting on the mode 2: in this way CVQT is reduced to a conditional measurement on the TWB followed by a 
displacement. Let a be again the state we wish to teleport and let the TWB be the entangled state shared between 
Alice and Bob. The conditioned state of mode 2 is then (we put p™, = 023) 



1 



p(z) 
(1-A 2 ) 



7175(2:) 

(1-A 2 ) 
irp(z) 



Tr 2 [e TWB n 2 (z)®I] 

Tr, ^ \hh\h) 33 (k\ 2 (k\DÏ(z) a T D* 2 {z) 



h,k 



A h+fe 2{k\D T 2 {z) a T D* 2 (z)\h) 2 \h) 33 (k\ , 



(7.18) 



h,k 
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with the double-homodyne density probability distribution p(z) given by 

p{z) -Tr 23 [ ftW Bn 2 (z)®I] . (7.19) 
Since [Z3 T (z)]* = D^z), in Eq. ( 17 .181 we can write 

2 {k\D T 2 (z) a T D* 2 {z)\h) 2 = 2{k\ [Dl(z)aD 2 (z)] T \h) 2 

= 2 {h\Dl(z) a D 2 (z)\k) 2 , (7.20) 

and, suppressing all the subscripts, we have 

Q 3 (z) = V \ h+k (h\D\z)oD{z)\k) \h)(k\ . (7.21) 

In order to have a full quantum teleportation, we must displace the state (I7.2H applying D(z), obtaining 

q' 3 (z) = D(z) g 3 (z) D\z) 
(1-A 2 ) 



irp(z) 

(1-A 2 ) 
irp(z) 



\ h+k {h\D\z) a D{z)\k) D{z)\h){k\D\z) (7.22) 



h,k 

2\ 



]T A' l+fc (Mz)W\Mz)) \M*)KM*)\ > ( 7 - 23 ) 



where we defined the new s.o.n.c. {\iph[z))}> with IV'íiÍ- 2 )) = Notice that Eq. M .221 can be written in 

operational form as 

= £ ^«A^l^)^») <fc|A at ^t (z) 

= ^—Ar D(z) X a>a DUz) a DU) X afa DUz) , (7.24) 
■np{z) 

which is the same as in Eq. (17. 10i . Finally, g out is obtained by means of Eq. (17. 131 . 

7.1.4 Wigner functions representation 

This Section addresses CVQT described in terms of Wigner functions. We first derive the teleported state Wigner 
function when the shared state has the general form given in Eq. (17.21 . then we specialize the results to the case of 
aTWB. 

Let W[ct](q;i) and W[g 23 ](a 2 , ot 3 ) be the Wigner functions associated to the states ( 17. U and (17.21 . respectively, 
where ah — k 2 (xh + i-Vh) ( see Chapter[0. Since the Wigner function corresponding to the POVM of ideal double- 
homodyne detection on mode 1 and 2 is 

W[U 12 (z)](ax, a 2 ) = \ 5 {2) {{a x - a*) - z) , (7.25) 
with z = k 2 (x + iy), using Eq. dl . lOli the double-homodyne density probability distribution ( I7.8l l reads 

p(z) = n 3 / d 2 a 1 d 2 a 2 d 2 a 3 W[a]{a 1 )W[g 23 ]{a 2 ,a 3 )W[U 12 {z)]{a ll a 2 )W[I 3 ]{a 3 ), (7.26) 
Jc 3 

while the conditioned state of mode 3 is 

W[g 3 (z)}(a 3 ) = ^- / d 2 a 1 d 2 a 2 W[a}(a 1 )W[g 23 }{a 2l a 3 )W[U 12 (z)}(a 1 ,a 2 )W[I 3 }(a 3 )^ (7.27) 
l>{~) Jc 2 

where W[l3](«3) = tt^ 1 . Thanks to Eq. (I7.25Í and after the integration over otx, Eq. J7.27I becomes 

W[g 3 {z)]{a 3 ) = f d 2 a 2 W{a}(a* 2 + z) W[g 23 ](a 2 , a 3 ) 



irp(z) 

^— í d 2 a 2 W[a](a 2 )W[g 23 }(a* 2 -z*,a 3 ). (7.28) 
np{z) Jc 
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Now we perform the displacement D(z), on mode 3, obtaining 

W[g' 3 (z)]{a 3 ) = — ^ ( d 2 a 2 W[a]{a 2 )W[g 23 ]{a* 2 ~z*,a 3 -z). (7.29) 



Trp(z) 

with £3(2) = D{z) Qs(z) D'(z), and where we used the property J 1 . 1 1 1 bl . The output state of CVQT is obtained 
integrating Eq. ( 17.291 1 over all the possible outcomes of the double homodyne detection, namely 

W[g out ](a 3 ) = í d 2 zp{z)W[g' 3 (z)](a 3 ). (7.30) 
Jc 

If the shared state is the TWB, the Wigner function reads as follows 

^m ^Í1^7 ] <731) 

witha, lfc = (fòe[a h ],Qm[a h ],$te[a k ),%m[ak}) T and 

I2 being the 2x2 identity matrix, a\_ — je ±2r , cr 3 — Diag(l, —1) is a Pauli matrix, and r — tanh -1 A is the 
squeezing parameter of the TWB. Substituting Eq. (I7.3H into Eq. ( 17.291 1 and integrating over z one has 

W[g out ](a 3 ) = f exp |-^_^| W [a](a 2 ) (7.33) 

exp {— w T S _1 tu} 



;2, 



W[D(w) a D* (w)](a 3 ) , (7.34) 



with w = a>2, w = Ç8le[w], ^sm[w]) T , and S = So is given in Eq. ( I7.16> . Finally, using Eq. (11.1 021 . we obtain 
the same density matrix as in Eq. 17. 151 . 



7.1.5 Teleportation fidelity 

Teleportation has occurred when the output signal g out is in the same quantum state of the unknown input a. 
Therefore, we need to define a quantity which gauges the similarity between a and g out . This task is achieved 
using the so called "fidelity" or "average fidelity" between the input and output state. When the input signal is a 
pure state a = \ip) the fidelity is defined in the following way 11521 1 

F = Tr[a g out ] = (ip\Qont\ip) , (7.35a) 
F = ix í d 2 a W[a](a) W[g out ]{a) , (7.35b) 

in terms of density matrix and Wigner function representation, respectively. F has the property that it equals 1 if 
and only if a is a pure state and g out = u; on the other hand, it equals if and only if the input and output states 
can be distinguished with certainty by some quantum measurement. In particular the average fidelity evaluates 
the extent at which all possible measurement statistics produceable by the output state match the corresponding 
statistics produceable by the input state. In order to explain this last consideration, let us consider the genèric 
POVM {IT a }, describing a certain observable, with measurement outcomes a. If the observable were performed 
on the input system, it would give a probability density for the outcomes a given by 

P ia (a)=Tt[aTl a }; (7.36) 

if the same observable were performed on the output system, it would give, instead, the probability density 

Pout(a)=Tr[É> out II a ]. (7.37) 

Now, a natural way to gauge the similarity of these two probability densities is by their overlap Q, defined as 
follows 

Q= í d 2 a y/P in (a)P mt (a) . (7.38) 

J£ 

1 When cr is not a pure state, a good measure for the fidelity is given by F = Tr y y/a gout y/5 . 
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It turns out that regardless of which observable is being considered Q> F and, moreover, one can show 1 150 1 that 
there exists an observable that gives precise equality in this expression. 

When the shared state is the TWB, substituting Eq. J7. 14b into the Eq. J7.35al >. one straightforward obtains 

F™ B (A) = F(X) = TT i- r = I±* . (7.39) 

The maximum average fidelity achievable by means of some classical (local) procedure to teleport a state 
is known as classical límit. This procedure should be characterized by a local measurement on the state to be 
teleported, a classical communication of the result, say 7Z, and, finally, a preparation stage at the receiver, according 
to a rule that associates a certain output state to 7Z such that fidelity is maximized. 

Let us suppose that Alice wishes to transmit to Bob an unknown coherent state without the resource of entan- 
glement, i.e. by no means of a shared entangled state. In such a case, we are interested in evaluating the maximum 
average fidelity achievable. First of all we assume that the coherent state is drawn from the set S constituted by the 
coherent states |/3), where the complex parameter f3 is distributed according to the Gaussian distribution 

p(0) = H e -n\l>\' t (7.40) 

SI being a real, positive parameter. Ultimately, of course, we would like to consider the case where Alice and Bob 
have no information about the drown coherent state: this is simply described by taking the limit íi — > 0. 

Alice's measurement for estimating the unknown parameter j3 when it is distributed according to a Gaussian 
distribution 11511 is the POVM {II Z } constructed from the coherent state projectors according to Ii z — 7T _1 |z)(z|: 
this kind of measurement is equivalent to optical heterodyning described in Section 15.6.41 where we send the 
vacuum in the other detector input port [see Eq. ( 15 . 1 151 with r = |0)(0|]. As in the case of the teleportation 
protocol, Alice's measurement outcome z is classically sent to Bob, that generates a new quantum state according 
to the rule z — > \f z ). Let us make no a priori restrictions on the states \f z ). 

Now, we find the maximum average fidelity F max (fi) Bob can achieve for a given íl. For a given strategy 
z ~~ * \.fz), the achievable average fidelity is lTÍ52l 

F(Q)= f d^p(P) f d 2 zp(z\(3)\(f z \(3)\ 2 = ^ í d 2 zd 2 (3e- n ^e-^ 2 \(m\ 2 

= ^L f d 2 ze-^ 2 (f z \O z (n)\fz) , (7.41) 

TT JC 

where p(z\j3) — Tr[|/3) (j3\ Ii z ] is the heterodyne probability density distribution and we defined the positive semi- 
definite Hermitian operator 

O z (ü)= J^d 2 {3exp{-(l + ïl)\{3\ 2 + 2ïïe[z*p}}\í3)(f3\ , (7.42) 

that depends only on the real parameter íl and the complex parameter z. It follows that 

(/,|0,(n)|/,)<max[0,(íï)] , (7.43) 

where max[AT] denotes the largest eigenvalue of the operator X. 

Now, for each z, Bob adjusts the state \ f z ) to be the eigenvector of O z (VÏ) with the largest eigenvalue. Then 
equality is achieved in Eq. ( 17.431 1 and it is just a question of being able to perform the integral in Eq. M A\\ . The 
first step in carrying this out is to find the eigenvector and the eigenvalue achieving equality in Eq. d7.43i . This is 
most easily evaluated by unitarily transforming O z (fi) into a new operator diagonal in the number basis, picking 
off the largest eigenvalue and transforming back to get the optimal \f z ) (we remember that eigenvalues are invariant 
under unitary transformations). 

In order to find the largest eigenvalue of O z (Cl), we consider the positive operator 



P= f d 2 (3e^ 1+n ^ 2 \P)((3\ 
Jc 



(7.44) 
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Since the number basis expansion of P has the following diagonal form 

H/31 2 



P 



d z l3e 



"(l+f2)l/3| 2 



E 



yrï! m! 



/T (/3*) m |n)(m| 



V n! m! Jo 



2tt 



d<j> e 



i(f>(n—m) 



dpp 



n+m+l -(2+£l)p 2 



|n)(m| 



27T(5„.., 



9 /*°° °° 

E^T / ^P 2n+1 e- (2+í1) " 2 |n)(n| =7 r]T(2 + fi) 



>+l) 



n)(n\ 



(7.45) 



in! (2 + n)- (n+1 *> 



its eigenvalues are {7r (2 + fi) and, then, (i(P) — 7r/(2 + fi), /.<?. the vacuum state's eigenvalue. 

Now consider the displaced operator 



,(fi) = D 



l + fi 



l + fi 



(7.46) 



where £>(y) is the Standard displacement operator. Using Eq. HA41 . one finds 



Q 2 (fi) = / d 2 /3exp{-(l + fi)|/3| 2 } 



y d 2 Ç exp j-(l + fi) 



1 + fi 

2 ï 



i + fi 



IO<ei 



exp 



exp 



l + fi 

\z\ 2 



l + fi 

and, substituting this into Eq. HAl\ . we have 



l + fi 

d 2 Ç exp { - (1 + m\ 2 + 2He[z*£]}|0tèl 
0,(íï) , 



(7.47) 



F(fi) 



< 



fi 

7T 2 
1 fi 



d z exp 



1 



7T 2 + fi 



d z exp 



l + fi 
fi 
l + fi 



4 



z\ 2 \{f z \Q z {ü)\f z 



l + fi 



fi 



Equality is obviously achieved in the previous expression by taking 



1 + fi 



\h) = D 

therefore the maximum average fidelity is given by 



|o) = 



i + fi 



1 + fi 

2 + fi 



(7.48) 



(7.49) 



(7.50) 



For fi — > oo we have F max (fi) — * 1 since this situation corresponds to the teleportation of a single known coherent 
state, a task that can be achieved classically by transmitting the value of the amplitude. On the other hand, in the 
limitfi — * 0, i.e. when the coherent state to be sent is drawn from a uniform distribution, we have F max (fi) — * 1/2. 

It should be noted that nothing in this argument depended upon the mean of the Gaussian distribution being 
j3 = 0: Bob would need to minimally modify his strategy to take into account Gaussians with a non-vacuum state 
mean, but the optimal fidelity would remain the same. 

7.1.6 Effect of noise 

In this Section we study CVQT assisted by a TWB propagating through a squeezed-thermal environment. Taking 
into account the results obtained in Section l4~4l the teleported state is now obtained from Eq. J7.34Í with 



S = S(r,7V th ,7V s )=4 



£ 2 o 



(7.51) 



80 



Chapter 7: Teleportation and telecloning 



E§, E 2 being given in Eqs. (14.451 . 

Finally, non-unit quantum efficiency i] in the joint measurement modifies the POVM, which becomes a Gaus- 
sian convolution of the ideal one, as pointed out in Chapter|5] In this case the output state is given by Eq. J7.34Í 
where 

E = S(T,N th ,N s ,r l )=[ 3 + •» 4E 2 +jD 2], (7.52) 



with D 2 = íl - r/Vr? Í96l. 



7.1.7 Optimized teleportation in the presence of noise 

In order to use CVQT as a resource for quantum information processing, we look for a class of squeezed states 
which achieves an average teleportation fidelity greater than the one obtained teleporting coherent states in the 
same conditions. If the input squeezed state to be teleported is u = \a, £) (a, £|, |a, £) ~ D(a)S(£,)\0), then the 
teleported state given by Eq. J7.34Í has average teleportation fidelity [see Eq. J7.35bH 

FtW = (V(e 2e + 4El)(e- 2 €+4S|))" 1 , (7.53) 

which attains its maximum 

F{\) = (1+4E2E3)- 1 , (7.54) 

when 




(7.55) 



For non squeezed environment N s —> we have E2 = E3, and thus then £ max — > 0, ï'.e. the input state that 
maximizes the average fidelity J7.53l > reduces to a coherent state. In other words, in a non squeezed environment 
the teleportation of coherent states is more effective than that of squeezed states. Moreover, Eq. 17.541 shows that 
meanwhile the TWB becomes separable, i.e. E 2 , E§ > 16 _1 [see Eqs. (I4.47H . one has F < 0.5. Finally, the 
asymptotic value of F for Tt — ► 00 is 

F (oo) = [2(l + N th )}-\ (7.56) 

which does not depend on the number of squeezed photons and is equal to 0.5 only if N t h = 0. This last result 
is equivalent to say that in the presence of a zero-temperature environment, no matter if it is squeezed or not, the 
TWB is non-separable at every time. In Fig. 17.21 we plot -Ftele as a function of Tt for different vàlues of A, -/V t h 
and N s . As N s increases, the nonclassicality of the thermal bath starts to affect the teleportation fidelity and we 
observe that the best results are obtained when the state to be teleported is the squeezed state that maximizes ( 17.53b . 
Furthermore the difference between the two fidelities increases as N s increases. Notice that there is an interval 
of vàlues for Tt such that the coherent state teleportation fidelity is less than the classical limit 0.5, although the 
shared state is still entangled. 



7.1.8 Teleportation improvement 

TWBs are produced either by degenerate (with additional beam splitters) or nondegenerate optical parametric 
amplifiers. The TWB parameter A = tanhr depends on the physical parameters as r oc X L> being the 
nonlinear susceptibility of the crystal used as amplifying medium and L the effective interaction length. For a given 
amplifier, the TWB parameter and thus the amount of entanglement are fixed. Therefore, since nonlinearities are 
small, and the crystal length cannot be increase at will, it is of interest to devise suitable quantum operations to 
increase entanglement and in turn to improve teleportation fidelity. 

In Section l6.2.1l we have seen that the nonlocal correlations of TWB are enhanced for small energies by means 
of the IPS process described in Section l5.4.1l motivated by this result, we will use the IPS state J5.48ft or, equiv- 
alently, ( 15. 511 as shared entangled state between Alice and Bob. In this case, Eqs. ( I7.35l > lead to the following 
expression for the average teleportation fidelity of coherent states 1 92 1 

p ÍA , 1 (1 + A)(l + Ar cff )(l - A 2 r cff )[2 - 2Àr cff + A 2 r cff ] 
lA ' Tcffj 2(l + A 2 r cff )[l + (l-r cff )A]{2-[2+(l-T eff )A]Ar cff } ' 1 ° j 

where r e g = 1 — 7/(1 — r) (see Section l5.4.H . In Fig. 17.31 we plot the average fidelity for different vàlues of 
r e ff: the IPS state improves the average fidelity of quantum teleportation when the energy of the incoming TWB 
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Figure 7.2: Plots of the average teleportation fidelity. The sòlid and the dashed lines represent squeezed and coherent state 
fidelity, respectively, for different vàlues of the number of squeezed photons N s : (a) N s = 0, (b) 0.1, (c) 0.3, (d) 0.7. In all the 
plots we put the TWB parameter A = 1.5 and number of thermal photons N t h = 0.5. The dot-dashed vertical line indicates 
the threshold Tt s for the separability of the shared state: when Tt > rí s the state is no more entangled. Notice that, in the case 
of squeezed state teleportation, the threshold for the separability corresponds to F — 0.5. 



is below a certain threshold, which depends on r e g and, in turn, on r and r/. When t c s approaches unit (when 
r) — > 1 and r — > 1), Eq. (I7.57Í reduces to the result obtained by Milburn et al. in Ref. |94| and the IPS average 
fidelity (line labeled with "a" in Fig. 17.21 is always greater than the fidelity í^twb(A) obtained with the TWB state 
[see Eq. J7.39H . However, a threshold value, Ath(r e ff), for the TWB parameter À appears when r c fr < 1: only if A 
is below this threshold the teleportation is actually improved [F(X, t c b) > í 1 TWB (A)], as shown in Fig. 17.31 Notice 
that, for r ff < 0.5, F(X, T c g) is always below F TWB (X). 

Ralph et al. demonstrated that entanglement is needed to achieve a fidelity greater than 1/2 1 153] and, using 
both the TWB and the IPS state ( I5.5H . this limit is always reached (see Fig. 17.31 . Nevertheless, we remember 
that in teleportation protocol the state to be teleported is destroyed during the measurement process performed by 
Alice, so that the only remaining copy is that obtained by Bob. When the initial state carries reserved information, 
it is important that the only existing copy will be the Bob's one. On the other hand, using the usual teleportation 
scheme, Bob cannot avoid the presence of an eavesdropper, which can clone the state, obviously introducing some 
error 1 154|, but he is able to to verify if his state was duplicated. This is possible by the analysis of the average 




Figure 7.3: On the left: IPS average fidelity F(\, r e s) as a function of the TWB parameter A for different vàlues of r e a = 
1— »7(1 — t): (a) T e s = 1. (b) 0.9, (c) 0.8, and (d) 0.5; the dashed line is the average fidelity F T wb(A) for teleportation with TWB. 
On the right: Threshold value Ath(ieff) on the TWB parameter x (sòlid line): when A < A t h we have F(X, r e ff) > Ftí/ B (X) 
and teleportation is improved. The dot-dashed line is A = 1/3, which corresponds to F T wb = 2/3: when fidelity is greater 
than 2/3 Bob is sure that his teleported state is the best existing copy of the initial state 1 155 |. The dashed line represents the 
vàlues A2/3(T e ff) giving an average fidelity F(X, t c b) = 2/3. When A2/3 < A < A t h both the teleportation is improved and 
the fidelity is greater than 2/3. 
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teleportation fidelity: when fidelity is greater than 2/3, Bob is sure that his state was not cloned [ 154-1 11551 . The 
dashed line in Fig. l7.3l (right) shows the vàlues A 2 /3(r c ff) which give an average fidelity J7.57Í equal to 2/3: notice 
that when A 2 /3 < A < A t h both the teleportation is improved and the the fidelity is greater than 2/3. Moreover, 
while the condition -F T wb(A) > 2/3 is satisfied only if A > 1/3, for the IPS state there exists a r e g -dependent 
interval of A vàlues (A 2 /3 < A < 1/3) for which teleportation can be considered secure [F(X, r c ft) > 2/3]. 



7.2 Quantum cloning 

A fundamental difference between classical and quantum information is that the latter cannot be perfectly copied, 
even in principle. This means that there exist no physical process that can produce perfect copies of genèric 
unknown quantum states. This so called no-cloning theorem emerges as an immediate consequence of the linearity 
of quantum dynamics 11401 114Í1 . Remarkably, if cloning was permitted, the Heisenberg uncertainty principle 
would be violated by measuring conjugate observables on many copies of a single quantum state. Nevertheless, 
even if perfect cloning is not possible, one may attempt to attain imperfect copies of genèric unknown quantum 
states. With an abuse of language, this is what is generally referred to as cloning process. With an n to m cloning 
process it is thus meant that m imperfect copies are produced from n identical original states (m > n). In this 
section we address the cloning issue for Gaussian states, first recalling the bounds that no-cloning theorem imposes 
in this case, then investigating some local and nonlocal (telecloning) cloning protocols. 



7.2.1 Optimal universal cloning 

The first concept to introduce is the universality of a cloning machine 1 156|. By universality we mean that the 
quality of the clones should be independent on the original states. As we have already seen in Section l7.1.5l we 
may use the fidelity as a measure of the similarity between two states, hence in a universal cloning machine every 
clone has the same fidelity with respect to the original state, independently of the original state itself. Furthermore, 
when the clones are equal one each other then we deal with symmetric cloning, while if we admit differences in 
the copies we have asymmetric cloning. Consider for the moment the first scenario. Universal cloning machines 
have been extensively studied for the case of discrete variables, for which it has been shown that the optimal n to 
m universal cloning machine of d-dimensional systems yields the fidelity 11571 : 

n(d- 1) +777(77+ 1) 

*= / —j-s • (/.58) 

m(n + d) 

In the limit of infinite dimensional systems one can show that the optimal universal cloner reduces to a classical 
probability distributor, attaining F = n/m 1 158 1, consistent with the d — > oo limit of Eq. J7.581 . This means that 
a universal continuous variable cloner behaves like a simple classical device that distributes the n original input 
states into n output states, chosen by chance between the possible m outputs and disregarding the remaining 
states. Nevertheless, a more interesting situation occur if we restrict the input states to the class of Gaussian states. 
Consider for the moment n identical arbitrary coherent states. The imperfection of the m copies may be regarded 
as an excess noise variance <r^ m in the quadratures, due to the nio m cloning process. Then, a procedure similar 
to what was done for qubits 1 159 1 allows one to estimate a lower bound cr^ m on the noise variance 1 160|. In fact, 
make use of the property that cascading two cloning processes results in a single cloning process whose excess 
noise variance is simply the sum of the variances of the two cloning. Then, the variance cr^ j of an optimal n to 
l cloning must satisfy < m + <rf n , (n < m < /)• In particular, if the m to l cloner is itself optimal and 
l — > oo, we have 

°n,oo < <T n ,m + ^m,oo • ( 7 - 59 ) 

As a matter of fact, a cloner that allows to build infinitely many copies corresponds to an optimal measurement 
of the original states, hence, with the aid of quantum estimation theory, one may identify '& 2 noo — X/n (we put 
fil = 2 -1 / 2 ). As a consequence, the lower bound we were looking for is given by 

o 777 — n 

< m = ■ (7-60) 

mn 

This result implies that the optimal cloning fidelity for coherent states is bounded by II 1601 

(7.61) 



77777 + 777 — 77 



Notice that this result does not depend on the amplitude of the input coherent states. If general squeezed states are 
considered, optimality can then be achieved only if the excess noise variance is squeezed by the same amount as 
the initial state, thus making the cloner state-dependent. 
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Figure 7.4: Scheme of local 1 to 2 cloning. 
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Figure 7.5: Schematic diagram of the 1 — > 2 telecloning scheme described in the text. 



Remarkably, optimal cloners achieving the fidelity given in Eq. Í7.6H may be implemented in an optical 
framework with the aid of only a phase-insensitive two-mode squeezer and a sequence of beam splitters 1 16l1 ll62l . 
As an example, consider the 1 to 2 cloner depicted in Fig. 17.41 Mode ao, excited in the state to be cloned, is sent 
to a two-mode squeezer with an ancillary mode «i. Using the notation introduced in Section n'.4.4l we have that 
the output mode b is given by 6 = / ia o + va \- Then a linear mixing of modes b and bi in a phase insensitive 
balanced beam splitter give rise to 

c = ^(^o + va\ + h) , ci = -^=(nao + va{- b{). (7.62) 

Now, considering vacuum inputs for modes a\ and b\, and squeezing parameters \i = v2, v = 1, it follows that 
(co) = (ci) = (ao). As a consequence, the scheme considered allows to copy the ampli tude of the original mode 
ao. Thus, if the latter is excited in a coherent state two clones are produced at the output modes cq and c\. The 
optimality of the cloner follows from the fact that the two-mode squeezing chosen yields an excess noise variance 
er 2 = 1/2 HI 631 - Finally, the generalization of this method allows to realize an optimal n to m local cloning 
machine. 



7.2.2 Telecloning 

The cloning process described above, even if local, may be applied in order to distribute quantum information 
among many distant parties, in what is called a quantum information network. Suppose that one wants to distribute 
the information stored into n states to m receivers. This may be achieved by two steps. One may first produce 
locally to copies of the original states by means of the cloning protocol presented above. Then, the teleportation 
of each copy, following the scheme described in Section 177X1 allows to attain the transfer of information 1 164|. 
This strategy has the obvious advantage to use only bipartite entangled sources. However, even in the absence of 
losses, it does not leave the receivers with m optimum clones of the original states, due to the non-unitary fidelity 
of the teleportation protocol in case of finite energy. This problem may be circumvented by pursuing a one-step 
strategy consisted of a nonlocal cloning. By this we mean that the cloning process is supported by a multipartite 
(m + n) entangled state which is distributed among all the parties involved. This so called telecloning process is 
thus nonlocal in the sense that it proceed along the lines of a natural generalization of the teleportation protocol to 
the many-recipient case 1 164|. To clarify this second scenario, let us describe now in details a 1 to 2 telecloning 
process based on the tripartite state \T) introduced in Eq. J2.49I II47I . 

A schematic diagram of the telecloning process is depicted in Fig. 17.51 After the preparation of the state |T), 
a joint measurement is made on the mode a\ and the mode b to be telecloned, which corresponds to the measure 



84 



Chapter 7: Teleportation and telecloning 



of the complex photocurrent Z = b + a\, as in the case of the teleportation protocol. The whole measurement is 
described by the POVM ( 15. 1 151 . acting on the mode oi, namely H(z) = tt^ 1 D(z) er T 'D'(z), where a is the the 
state to be teleported and cloned. The probability distribution of the outcomes is given by 

P(z) = Tr 123 [\T)(T\Il(z)®h®h] 

The conditional state of the mode a 2 and 03 after the outcome z is given by 
1 



P[z) 



Tri[|T)(T|n(z)(8)I 2 ®l3] 



1 1 yy^ f \ N? +k N§ +l l(p + q y.(k + iy. 



P(z) tt(1 + Ni) ^ V 0- + iVi)P+9+ fe+í V p! <?! fe! í! 

x (k + l\D{z)a T D^{z)\p + q) \p,q)(k,l\ . (7.64) 

After the measurement, the conditional state should be transformed by a further unitary operation, depending on 
the outcome of the measurement. In our case, this is a two-mode product displacement U z = D^z) ® D^(z). 
This is a local transformation which generalizes to two modes the procedure already used in the original CVQT 
protocol described in Section lTTl The overall state of the two modes is obtained by averaging over the possible 
outcomes 



023 = / d 2 Z P(z) T z 

Jc 



where r z = U z g z U\. 

If b is excited in a coherent state a = \a)(a\, then the probability of the outcomes is given by 



p ^ = ^IT^) exp rTTÏvri- (7 ' 65) 

Moreover, since the POVM is pure also the conditional state is pure. Is this way we have that g z — \ip z ))((ip z \ is 
the product of two states, namely 

\4>z)) = \{a + z*) e 2 ) <g> \(a + z*) e 3 ) , (7.66) 

where 

e *=Vï+k (/l=2 ' 3) - (7 - 67) 

Correspondingly, we have t z = U z \ipz)){{ipz\ U\ with 

U z |V>*>) = \ae 2 + z* (e 2 - 1)) ® \as 3 + z* (e 3 - 1)) . (7.68) 
The partial traces p 2 = T-3^3] and 03 = Tr 2 [g> 2 3] read as follows 

Qh = d 2 z P a (z) \a.Eh + z* (eh - l)){oc£h + z* (e h - 1)| . (7.69) 
J c 

From the teleported states in J7.69I we see that, depending on the vàlues of the coupling constants of the Hamil- 
tonian (I2.48Í the two clones can either be equal one to each other or be different. In other words, a remarkable 
feature of this scheme is that it is suitable to realize both symmetric, when 7Y 2 = N3 = N, and asymmetric cloning, 
N-2 ^ N 3 . This arise as a consequence of the possible asymmetry of the state that supports the teleportation. 
Let us first consider the symmetric cloning. According to Eq. J2.50Í the condition N2 = N 3 = N holds when 



■on 71 AT 4 7 i 72 
cos(ÍM) = - — 12 ; — ,2 , N = — 12 ; 12T2 • ( 7 - 7 °) 

2 72 - 7ir ( 2 72 - 7irr 



Since \(z \z )\ = exp{— \z —z | }, the fidelity of the clones is given by (we put e 2 = ^3 = e) 

f = L - {- wíf } -* h« + ^ - i)2 í 

-1 



= ( 2 + 3iV - 2y/N(2N +1)1 . (7.71) 
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As we expect from a proper cloning machine, the fidelity is independent of the amplitude of the initial signal and 
for < N < 4 it is larger than the classical limit F = 1/2. Notice that the transformation U z performed after the 
conditional measurement, is the only one assuring that the output fidelity is independent of the amplitude of the 
initial state. Exploiting Eq. J7.7H we can see that the fidelity reaches its maximum F = 2/2, for N = 1/2 which 
means, according to Eq. J7.70Í . that the physical system allows an optimal cloning when its coupling constants are 
chosen in such a way that I71/72I = (6 - V32) 1 / 2 ~ 0.586 . The total mean photon number required to reach the 
optimal telecloning is thus Ai + A2 + A3 = 2, hence, as we claimed above, it can be achieved without the need of 
infinite energy. The scheme presented is analog to that of Ref. 1 125 1 in the absence of an amplification process for 
the signal. There, the telecloning is supported by a state similar to the one given in Eq. 12.461 . where at the input 
ports of the tritter only two squeezed states are involved, the third mode being a vacuum state. Both the protocols 
described here and in Ref. 11251 achieve the optimality relying on minimal energètic resources, i.e. the total mean 
photon number is 2 in both cases. Notice that a generalization to realize a 1 — * m telecloning machine can be 
realized upon the implementation of SU(p, 1) Hamiltonian introduced in Section lT.4.51 In fact, having at disposal 
a 1 + m multipartite entangled state of the form ( 11.941 . it is straightforward to show that a measurement of Z 
on the mode to be telecloned and the sum-mode of 11.941 . followed by a local multimode displacement operation 
provides optimal clones in the remaining m modes. 

Let us now consider the asymmetric case. For A2 7^ A3 the fidelities of the two clones J7.69I are given by 



where h, k = 2, 3 (h 7^ fc). A question arises whether it is possible to tune the coupling constants so as to obtain 
a fidelity larger than the bound F = 2/3 for one of the clones, say Q2, while accepting a decreased fidelity for 
the other clone. In particular if we impose F3 = 1/2, i.e. the minimum value to assure the genuine quantum 
nature of the telecloning protocol, we can maximize F2 by varying the value of the coupling constants 71 and 
72. The maximum value turns out to be i^.max = 4/5 and it corresponds to the choice A3 = 1/4 and A2 = 1. 
More generally one can fix F3, then the maximum value of F2 is obtained choosing A2 = (I/Í3 — 1) and 
-W3 = (4/.F3 — 4) _1 . The relation between the fidelities is then 




(7.72) 



F 2 = 4 



(1 - £g) 
(4-3F 3 ) ' 



(7.73) 



which shows that F2 is a decreasing function of F3 and that 2/3 < F2 < 4/5 when 1/2 < F3 < 2/3 . The 
sum of the two fidelities F% + F$ = 1 + 3F2F3/4 is maximized in the symmetric case in which optimal fidelity 
Fi = Fs = 2/3 can be reached. The role of g 2 and Q3 can be exchanged, and the above considerations still hold. 
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Chapter 8 

State engineering 



In this Chapter we analyze the use of conditional measurements on entangled twin-beam state (TWB) of radiation 
to engineer quantum states, i.e. to produce, manipulate, and transmit nonclassical light. In particular, we will 
focus our attention on realistic measurement schemes, feasible with current technology, and will take into account 
imperfections of the apparata such as quantum efficiency and finite resolution. 

The reason to choose TWB as entangled resource for conditional measurements is twofold. On one hand, 
TWBs are the natural generalization to continuous variable (CV) systems of Bell states, i.e. maximally entangled 
states for qubit systems. On the other hand TWBs are CV entangled states that can be reliably produced with 
current technology, either by parametric downconversion of the vacuum in a nondegenerate parametric amplifier 
1 165 1, or by mixing two squeezed vacua from a couple of degenerate parametric amplifiers in a balanced beam 
splitter l27ll26l . 

The first kind of measurement we analyze is on/off photodetection, which provides the generation of condi- 
tional nonclassical mixtures, which are not destroyed by decoherence induced by noise and permits a robust test 
of the quantum nature of light. 

The second apparatus is homodyne detection, which represents a tunable source of squeezed light, with high 
conditional probability and robustness to experimental imperfections, such non-unit quantum efficiency and finite 
resolution. 

The third kind of measurement is that of the normal operator Z = b + c\ b and c being two modes of the field, 
as described in Section l5~6l In our case one of the two modes is a beam of the TWB, whereas the second one, 
usually referred to as the probe of the measurement, is excited in a given reference state. This approach allows to 
describe CV quantum teleportation as a conditional measurement, and to easily evaluate the degrading effects of 
finite amount of entanglement, decoherence due to losses, and imperfect detection \ 96\. 



8.1 Conditional quantum state engineering 



The general measurement scheme we are going to consider is schematically depicted in Fig. 18. II The entangled 
state subjected to the conditional measurement is the TWB |A)), A = \/l — A 2 A a1< \ with A = tanhr assumed 
as real. A measurement, performed on one of the two modes, reduces the other one accordingly to the projection 
postulate. Each possible outcome x of such a measurement occurs with probability P x , and corresponds to a 
conditional state a x on the other subsystem (Fig. l8.1> . Upon denoting by H x the POVM of the measurement 1 we 
have 



Px = Tr at 



OO 

|A»«A|I®n a (1-A 2 )V A 2 "(g|n :c |ç) = (l-A 2 )Tr fc [A 2bTb n-] 



26+6 - 



q=Q 



(8.1) 



and 



Qx 



P 



■ Tn, 



|A))((A|i®n 3 



1-A 2 

Px 



£ X P+CI (p\K\q) \p)(q\ 



Tr fc [A 2btí TÇ] 



(8.2) 



Notice that in the second line of Eq. J8.2t H x should be meant as an operator acting on the Hilbert space 7í a of 
the mode a. Our scheme is general enough to include the possibility of performing any unitary operation on the 
beam subjected to the measurement. In fact, if E x is the original POVM and V the unitary, the overall measurement 

'in this Chapter, in order to simplify notation, we denote the dependence of the element of the POVM {H X } X £ X on the outcome x as a 
subscript rather than on parenthesis as in Section l5"Tl 
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n, = v%v 

Figure 8.1: Scheme for quantum state engineering assisted by entanglement. 

process is described by H. x — V'E X V, which is again a POVM. In the following we always consider V = I, i.e. no 
transformation before the measurement. A further generalization consists in sending the result of the measurement 
(by classical communication) to the reduced state location and then performing a conditional unitary operation U x 
on the conditional state, eventually leading to the state a x = XJ x q x XJ\. This degree of freedom will be used in 
Section l8.1.3l where we re-analyze CV quantum teleportation as a conditional measurement. 



8.1.1 On/off photodetection 



By looking at the expression of TWB in the Fock basis, |A)) = \/\ — A 2 V \ q \a) \a). or at Eq. (I8.2Í it is apparent 
that ideal photocounting on one of the two beams, described by the POVM 11/. = (k\, is a conditional source 
of Fock number state \k), which would be produced with a conditional probability Pk = (1 — A 2 )A fc . However, 
realistic photocounting can be very challenging experimentally, therefore we consider the situation in which one 
of the two beams, say mode b, is revealed by an avalanche on/off photodetector (see Section l5.3.2> . The action of 
an on/off detector is described by the two-value POVM {ITo (rj), 111(77)} given in Eq. J5.20Í . The outcome "1" (i.e. 
registering a "click" corresponding to one or more incoming photons) occurs with probability 

Pl = ((A|^n l(ï7 )|A)) = T-^y = ^ , (8.3) 
with N\ = 2A 2 /(1 — A 2 ), and correspondingly, the conditional output state for the mode a is given by 1 166 1 

1 _ \2 00 

ei = -p- £A afc [i-(i -»,)*] l*X*l- (8-4) 
1 fe=i 

The density matrix in Eq. ( 18.41 describes a mixture: a pseMc/o-thermal state where the vacuum component has been 
removed by the conditional measurement. Such a state is highly nonclassical, as also discussed in Ref. 1 8 1. Notice 
that the nonclassicality is present only when the state exiting the amplifier is entangled. In the limit of low TWB 
energy the conditional state Q\ approaches the number state |1)(1| with one photon. 

The Wigner function W[f?i](a) of g\ exhibits negative vàlues for any value of A and 77. In particular, in the 
origin of the phase space we have 

-[*) = -ïï^fSRï- (8 - 5 > 

One can see that also the generalized Wigner function for s-ordering 

W s [g 1 ](a) = -— í í 7 ^[ft](T) exp(--|a- 7 | 2 l , 
shows negative vàlues for s E (—1, 0). In particular one has 

Ws[0l]{O) = -n(l + N x - S )l2(l + N x - S )- V N x (l + s)] ' (8 - 6) 

A good measure of nonclassicality is given by the lowest index s* for which W s is a well-behaved probability, 
i.e. regular and positive definite |83|. Eq. ( 18. 6> says that for gi we have s* = —1, that is gi describes a state as 
nonclassical as a Fock number state. 
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Since the Fano factor of Qi is given by 



F 



_ ({tfb-(bïb)} 2 ) _ (2 + iV A ) 
(òtò) 2 



4 (2 + N x ) 



2 + T]N x 4 + Nx (4 + í7ÍV à ) 



(8.7) 



we have that the beam ò is always subPossonian for (at least) N\ < 2. The verification of nonclassicality can be 
performed, for any value of the gain, by checking the negativity of the Wigner function through quantum homodyne 
tomography 1 166 1, and in the low gain regime, also by verifying the subPoissonian character by measuring the Fano 
factor via direct noise detection I 



I167II168I . 

Note that besides quantum efficiency, i.e. lost photons, the performance of a realistic photodetector may be 
degraded by the presence of dark-counts, i.e. by "clicks" that do not correspond to any incoming photon. In order 
to take into account both these effects we should describe the detector by the POVM ( 15.221 rather than ( 15.20b . 
However, at optical freqüències the number of dark counts is small and we are not going here to take into account 
this effect, which have been analyzed in details in Ref. 1 166). 



8.1.2 Homodyne detection 

In this Section we consider the kind of conditional state that can be obtained by homodyne detection on one of the 
two beams of the TWB . We will show that they are squeezed states. We first consider ideal homodyne detection 
described by the POVM Ií x = \x){x\ where \x) denotes the eigenstate Jl . 1 5I > of the quadratures x = i(a + a)) 
(throughout the Section we use k,\ — n-i — 1) and where, without loss of generality, we have chosen a zero 
reference phase (see Section l5~5l for details about homodyne detection). Then, in the second part of the Section we 
will consider two kinds of imperfections: non-unit quantum efficiency and finite resolution. As we will see, the 
main effect of the conditional measurement, i.e the generation of squeezing, holds also for these realistic scenarios. 

The probability of obtaining the outcome x from a homodyne detection on the mode ò is obtained from 
Eq. <|8J). We have 



OC - 



P x = {\-\f^ q \{x\q)f = -—^ 9 \- T ^ \ , (8.8) 

g=0 V 27R7 A L 



where 



„ 1 1 + A 2 1, 
°l = -,— = 4 il + Nx) - (8 ' 9) 

P x is Gaussian with variance that increases as À is approaching unit. In the (unphysical) limit A — * 1, i.e. infinite 
gain of the amplifier, the distribution for x is uniform over the real axis. The conditional output state is given by 
Eq. (I8.2Í . and, since H x is a pure POVM, it is a pure state g x = \ip x ) (ijj x | where 



/ 1 _ V2 00 

l^> = v-p— AB+B !*> = !>* !*>■ (8 - 10) 



Px 

k=0 



The coefficients of \tp x ) in the Fock basis are given by 



, n .4,1/4 ÍX 2 \ k/2 Hk{V2x) f 2A 2 o; 2 



(8.11) 



which means that is a squeezed state of the form 

|V*> = D(a x )S(()\0) , (8.12) 

where 



2xA _ Xy/N x {N x + 2) 



1 + A 2 1 + N\ 



(8.13a) 



C = tanbT^A 2 ) = tanh" 1 ( — ) , (8.13b) 



N x + 2 



and the quadrature fluctuations are given by 



Ax 2 = - — - — , Aï/ 2 = -(1 + N X ) . (8.14) 

4 1 + N x 4 
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Notice that (i) the amount of squeezing is independent on the outcome of the measurement, which only influences 
the coherent amplitude; (ii) according to Eq. J8.8Í the most probable conditional state is a squeezed vacuum. The 
average number of photon of the conditional state is given by 

2 N X (2 + N X ) 1 



N x = (ip x \a* a\ip x ) = x 2 



Ni 



(1 + Nx) 2 4 1 + N x ' 
The conservation of energy may be explicitly checked by averaging over the possible outcomes, namely 



dx P-r Nr,. 



1 



N 2 



Nx(2 + N x ) 



1 



41 



Nx 



(8.15) 



(8.16) 



Nx (1 + ÍVa) 2 

which correctly reproduces the number of photon pertaining each part of the TWB. 

We now take into account the effects of non-unit quantum efficiency 77 at the homodyne detector on the condi- 
tional state. We anticipate that g xrj will be no longer pure states, and in particular they will not be squeezed states 
of the form (18. 121 . Nevertheless, the conditional output states still exhibit squeezing, i.e. quadrature fluctuations 
below the coherent level, for any value of the outcome x, and for 77 > 1/2. The POVM of a homodyne detector 
with quantum efficiency 77 is given in Eq. J5.67b . Since the nonideal POVM is a Gaussian convolution of the ideal 
POVM, the main effect is that ïl XT] is no longer a pure orthogonal POVM. The probability P xr) of obtaining the 
outcome x is still a Gaussian with variance 



(8.17) 



where S 2 is given in Eq. J5.68I . The conditional output state is again given by Eq. 18.21 . After some àlgebra we 
get the mateix element in the Fock basis 



\/2 n+m n\m\ 



(1-A 2 ) A"+" / qp-TKl-AS)] 
(n\g XfJ \m) = — , , — \j ^ exp 



1-A 2 



4?7 2 A 2 a; 2 



1 - A 2 (l - 2r/) 



mm \rti.n 



. 1 1 1. i i 1 . ; 1 \ , v / \ 

£ 2fefc (T) U) Vv m+n - 2k H m+ n-2k (^x) , (8.18) 



where H n (x) is the n-th Hermite polynomials. The quadrature fluctuations are now given by 

2 _ ll+JV A (l-77) . o 1 



Axí 



2 -^(1 + Nx). 



(8.19) 



4 1 + 77 Nx 

As a matter of fact, Ay 2 is independent on 77, whereas Ax 2 increases for decreasing 77. Therefore, the conditional 
output g xrj is no longer a minimum uncertainty state. However, for 77 large enough we still observe squeezing in 
the direction individuated by the measured quadrature. We have that the conditional state is a general Gaussian 
state of the form ( 12. 151 with an average number of thermal photons given by 



^Vth = 




(1 + N x )[l + N x (l - v)] 



l + r]N x 



- 1 



and with amplitude and squeezing parameters 



Vx/NxiNx 



x-XT] 



6, 



1 



In 



(1 + Aí a )(1 + 77ÍV a ) 
1 + Nx(l~v) 



(8.20) 



(8.21) 



1 + vNx 

From Eqs. d8 . 1 91 and d8 .2 H we notice that g xv shows squeezing if 77 > 1/2, independently on the actual value x 
of the homodyne outcome. 

The outcome of homodyne detection is, in principle, continuously distributed over the real axis. However, in 
practice, one has always to discretize data, mostly because of finite experimental resolution. The POVM describing 
homodyne detection with binned data is given by 



1 



dt n 



x-S/2 



where Ií tv is given in Eq. J5.67Í . and 5 is the width of the bins. The probability distribution is now given by 




(8.22) 

i 

(8.23) 
(8.24) 



8.1 Conditional quantum state engineering 



91 



where A 2 ?J is given in Eq. J8.17t and 

Erf(x) = —= \ dt e~* 
V 71 " Jo 

denotes the error function. The conditional state is modified accordingly. Concerning the quadrature fluctuations 
of the conditional state we have, up to second order in S, 

which is below the coherent level for rj > 1 /2 and for 



_ 1 / 3(1 + 7 ? JV a )(2t 7 -1) 

W <XS= 8)I rf(N x + 2) ■ (8 ' 26) 

Therefore, the effect of finite resolution is that the conditional output is squeezed only for the subset \x\ < x$ of 
the possible outcomes which, however, represents the range where the probability is higher |96|. 



8.1.3 Joint measurement of two-mode quadratures 

In this Section we assume that mode b is subjected to the measurement of the the real and the imaginary part of the 
complex operator Z = b + c\ where c is an additional mode excited in a reference state S. As we have seen in 
Section l5"6l this kind of measurement is described by the POVM 



n, 



-D(a) S T D\a). 



(8.27) 



The present scheme is equivalent to that of CV teleportation, which, as pointed out in Section l7.1.1l can be viewed 
as a conditional measurement, with the state to be teleported playing the role of the reference state S of the 
apparatus. In order to complete the analogy we assume that the result of the measurement is classically transmitted 
to the receiver's location, and that a displacement operation D^(a) is performed on the conditional state g a . 
Eqs. diS. ll i and J8.2I > are rewritten as follows 



Pa 



Qa 



(1 - A 2 )Tr 2 



2a T a 



Tr 2 A 2ot X 
a a = D'(oe) g a D(a) 

while the teleported state is the average over all the possible outcomes, i.e. 





) X Ja D(a) 


Tr 2 







(8.28) 
(8.29) 

(8.30) 



d ap a a a 



d 2 a D\a) «A|I ® n Q |A)) D(a) 



After performing the partial trace, and some àlgebra, one has 



d a 
c 7T(J - 



exp 



D(a)SD\a) 



(8.31) 



(8.32) 



where a 2 _ = 1 + N x - \fN x {N x + 2), i.e. the result of SectionlTTl 
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